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FOREWORD 


This final report presents work which was conducted for Langley Research 
Center (LaRC) in response to requirements of Contract. NA31-16983. The work 
presented was performed by REMIECB. Inc* Huntsville* Alabama and is entitled 
"MINIVER Upgrade For The AVID System". The final report consists of three 
volumes. 

Volume I: LANM1N User's Manual 

Volume II: LANMIN Input Guide 

Volume III: EXITS User's and Input Guide 

The NASA technical coordination for this study was provided by Ms. Kathryn 
E. Wurster of the Vehicle Analysis Branch of the Space Systems Division. 
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Fig. 1.2 LANMIN Macro Flow Diagram 
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Section 2.0 
CARO INPUT 

r 


Tho card inpot to LANMIN has been broken into sixteen functional aresa. 
These areas are as follows: 

(1) Title 

(2) Timing Paraaeters and Print Control 

(3) Trajectory Data 

(4) Atmosphere Data 

(5) Flowfield Options 

(6) Flow Angle Data 

(7) Angle of Attack Option 

(8) Crossflow Data 

(9) Transition Criteria 

(10) Heat Transfer Options 

(11) Heating Multipliers 

(12) Transforaations 

(13) Surface Condition 

(14) Geoaetry Data 

(15) Initial Conditions 
(14) Control Paraaeters 

Most of the input for these areas involves specifying the values for the *W* 
A definition of all the eleaents of the *1* array is given in Appendix 
A. 
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A 


1 
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(1) TITLE 
Format (18A4) 


Cat* Description 

(2) IlMIUfi PARAMETERS AND print control 
F ormat (3F20.6) 


CARD 

NO. 

COLUMN 


INPUT VALUE 

1 

1-20 

tl 

Initial time, sec (V(l)) 


21-40 

Atl 

Printout interval 1, sec (»(2)) 


41-60 

t2 

Second time, sec (W(3)) 

2 

1-20 

At2 

Printout interval 2, sec (W(4)) 


21-40 

t3 

Third time, sec (f(5)) 


41-60 

At3 

Printout interval 3, sec (W(6)) 

3 

1-20 

t4 

Fourth time, sec (W( 7 ) ) 


21-40 

*CALC 

Calculation interval factor (W(8)) 

Note: If less than three printout intervals are desired, 

the unused parameter must be input as zero* 

(3) TRAJECTORY DATA 



Format 

(3F20.6, 2F10.6) 


CARD 

NO. 

COLUMN 


INPUT VALUE 

1 

1-20 

Number 

of time dependent table entries (W(50) limit * 50) 

2 

1-20 

Time (t), sec (f (51)— W(100) ) 


21-40 

41-60 

61-70 

71-80 


Altitude (a), ft (f (101)— f(150)) 

Velocity (v), ft/aec (W(151)— 1(200) ) 

Angle of Attack (a), degree* (W(211)— W(260) ) 
Yaw Angle (0), degrees (W( 651) — 1(700)) 
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(4) A TMOSPHERE DATA 
Format 5(13, F10.6) 


LOCATION 
FORMAT (13) 


OPTION OR DATA 
FORMAT (F10.6) 


010 


Ataosphere/Freestrean Properties Option Flag (ATFLAG) 

°* prop ® rt1 ®® defined nalng the U& P.S. 

RtlBdftr . fl AlMRlfilUtXR* Trajectory data are required input. 

1. find tnnnel input option - 

Freeatream static teaperature and pressure are input in 

ct°ui°r .. n^f unction 

or tine. The tine data and freeatreaa velocity data are 

input in the TRAJECTORY DATA card act with the altitud^ 

'*00 ) wt equal to zero* 

2. Freeatreaa static teaperature and pressure are input in 

if ‘.1 5 w 7 550 ' — 

aoi iii* d ^* j Altitod ® ▼•lues are input in locations 

iniii *‘ J :; t0ry dgU (tlae ' ® ltitBd ®* velocity and 
angle of attack) are required input* 

4. Freeatreaa properties are defined using Ifee 1363 Patrick 

iSSS *: '££.*** 

5 ‘ SfwMc^A^rr 1 ** *£* def ined n,ln * IfeE 1221 tt utdenbert 
Reference AtF°fPftcrc . Trajectory data are required input* 

*' fr ? rti " ard deflaed n » lB « 21s 1221 aatoatm 

iRi £lZ AiS£l£&£££. Trajectory data are required input. 

7 ’ !!f!!** rot, ?. prop ? rt1 ®* lre defined using m> 13 7J Vandenbe*. 
£alA £lX Almistel*. Trajectory data ereTe quite d input! 

H^ e n!nt* Pr T rtle V ,e deflned Min * Sfi 1221 Konnedv 
“ AtMPtPhf fg * Trajectory data are required input. 

9 * t!trn re *T propertl ®» ®*® defined using IJg 1221 Kennedy 
i?9 Id Day Atnoschera. a.s. - . . 


400 

Freeatreaa Properties Data Input as Follows: 
Nuaber of table entries (Linit 50) Enalts 

401-450 

Table of altitude data, ft 

451-500 

Table of freestrean temperatures, *R 

501-550 

Table of ireeetrean pressures, lb/sq ft 



(5) JELQHEJJBLB .OPTIONS 
Format 5(13, F10.6) 


LOCATION 
FORMAT (13) 

03 1-03 6 

046-048 


OPTION OR DATA 
FORMAT (F10.6) 


Flowfield and Local Preaanra Option Flags (FF) 

Flowfiald and preaanra option flags are inpnt beginning in 

Location 031. A total of 9 locations are available 

(031-036 and 046-048). The firat f Wield optical., 

?? 10 Locatloa 031 » followed by additional * 

nn^hf pte “ ttr * °P tlon f • • reeded. Location 

oreasure^tfl ^ *? #d *“ • ,ceadln * ord ®*‘ Flowfield and 
pressure flags always occur in pairs except for FF ■ 39 

(S..,t Cylinder) e.d FP - 29 (Pr.edu 
Flowfield Options: 

35. Sharp Wedge Shock Angle. 

The shock angle for a sharp wedge is deternined from a 
bniit-in table of sharp wedge shock angle as a function 

iLli it* A h "“ b ? r * nd the wed * e or flow Reflection 
•ogle. The wedge angle is inpnt. 

36. Sharp Cone Shock Angle. 

The shock angle for a sharp cone is deternined from an 
internal table as a function of the cone half-angle and 
the upstream Mach number. Ihe cone half-angle is input. 

38. Oblique and Normal Shock. 

The actual shock angle is input, 
a normal shock. 


An angle of 90 represents 


39. Parallel Shock. 

U«ed alone or as a last option to Indicate that a swept 

flM 8 Wed) Par Tk 101 , J oc V ol,ltlon *• Ro*ireu (no pressure 
i a * aeedod) * The angle-of-attsck or shock angle is incut 

cyliad#r th0 ln P nt * a fl® i* the (90-A) where * 
k is the sweep angle in degrees. 




LOCATION OPTION OR DATA 

FORMAT (13) FORMAT (F10.6) 


031-031!! Local Pressure Options: 

046-048 14. INPUT Cp as a function of M a , for the location V(649) • 1.0, 

Data Stored in TCPM< ) - Cp, TMCP( ) - H m , NCPlfT - number of 
Mach Numbers. 

13. Tangent Wedge Pressure Coefficient. 

The local pressure coefficient for a tangent wedge is deter- 
mined from a built-in table as a function of the wedge angle 
and the upstream Mscb number. The wedge angle is input, 

16. Tangent Cone Preasure Coefficient. 

The local pressure coefficient for a tangent cone is determined 
from a built-in table as a function of the cone half-angle and 
the upstream Maoh number. The cone half-angle is input. 

17. Oblique Surface Pressure. 

Provides pressure solution for a surface whose shock angle is 
slightly greater than the surface angle-of-attack or flow de- 
flection angle. The shock angle must be known and input under 
the flowfield shock option (flowfield option 38). The input 
pressure angle ia the surface angle-of-attack or flow deflection 
angle. 

18. Modified Newtonian Pressure - Oblique Shock. 

Local pressure is calculated as a function of the local slope. 
Local total pressure is used with this option. 

The local body angle is input. 

Prandtl-Meyer Expansion: 

29. This option provides a Prandtl-Meyer expansion and is used fol- 
lowing any one of the local pressure options. The input angle 
is the angular difference between the two slopes (flow turning 
angle). 


(6) FLOW ANGLE DATA 
Format 5(13, F10.6) 


LOCATION OPTION OR DATA 

FORMAT (13) FORMAT (F10.6) 


037-045 Flow Angle Data Input. 

The angles required as input by the flowfield options are input 
beginning in location 037. These data must be input in the same 
order as the flowfield option location requiring their input, 
(i.e., the angle input in location 037 must correspond to flow- 
field option in location 031, etc.). All angles are input in 
degrees. 
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(7) ANGLE OF ATTAC K OPTION 
Format 5 <13, F30.6) 


LOCATION 
FORMAT (13) 


261 

262 

263 

264 

265 

266 

267 

268 


OPTION OR OATA 
FORMAT <1106 ) 


Anglcof-Attack Option 

Angles input in locations 03 1 and 038 are adjusted within 
tba program to account for angle-of-at tack. If it is 
desirable for otlu r angles input in locations 039 through 
045 to be adjusted for angle-of-attack, the following 
control parameters must be input: 

1> Initialized in Main 

2. Initialized in Main 


3. 

4. 

5. 

6 . 

7. 

8 . 

9. 


Varies angle in Location 039 
Varies angle in Location 040 
Variea angle in Location 041 
Varies angle in Location 042 
Varies angle ia Location 043 
Varies >ngle in Location 044 
Varies angle in Location 045 


(input data) 
(input data) 
(input dat> ,* 
(input data) 
(input data) 
(input data) 
(input data) 


269 


(8) CROSSFLOW DATA 
>'\.»nut 3(J '•*, 5) 


LOCATION OPTION OR DATA 

FORMAT (13) FORMAT (F10.6) 


201 Crossflow Option Flsg (CFFLG) 

1. Constant Width Roctsngls (Idssl Gas). 

This method assumes a eonstant width rectangle and ideal 
gas ehordwise velocity gradient. Required inpnta are 
corner radius. Ro # and width. D o « 

2. Conatant Width Rectangle (Real Oas). 

This method assumes a constant width rectangle and 
real gas ehordwise velocity gradient. Required inpnta 
are width. Do. and the velocity gradient correction 
factor. V. The term U can he varied from 0.31 for a 
flat surface to 1.0 for a swept cylinder. 

3. Sharp Edged Delta Configurations (Ideal Gas). 

This method assumes a sharp edged delta configuration 
and ideal gas ehordwise velocity gradient. Required 
input is the delta sweep angle X. 

4. Delta Configuration (Real Gas). 

This method assumes a delta configuration and real gaa 
ehordwise velocity gradient. Required Inputs are the 
delta sweep angle. X. and the velocity gradient 
correction factor, ,0. A sharp edged delta configuration 
is represented by U a 0.31« An input of M.O produces 
crossflow on a pointed cone whose half angle is (90-X). 


Crossflow Data Input 

202 Rectangle Width. D Q (ft) 

203 Delta wing sweep angle. X (deg) 

204 Reel Gas velocity gradient correction factor. U 

205 Rectangle corner radius. R Q (ft) 


(9) TRANSITION rare 
Porn at 5(13, F10.6) 




1EM 


LOCATION 
FORMAT (13) 


OPTION OR DATA 
FORMAT (F10.6) 


027 


014 


020 


Transition Criteria Option Flag (TRFLAG) 

1. This option allows transition fron laminar to turbulent flow to 
oocnr at a particular time in the trajectory. The tiie (.H) 
when transition begins is input in location 014. The time (sec) 
when fully turbulent flow is reached is input in location 020. * 

*• 1 l i2Lzr “ ° p,1M *• •* mp ' u 

3. Transition based on the Reynolds number R . . Input requirements 

R fy nold * (PARA1, location 014) and the fully 

turbulent Reynolds number (PARA2, location 020). 7 

4. Transition based on R q . Input requirements are onset 

Reynolds number (PARAS, location 014) and the fully turbulent 
Reynolds number (PARA2, location 020). 7 

5 * T'!!! itl0a ?!* ed ♦ (MDAC “ E transition parameter). Required 

inputs are the onset value of * (PARA1, location 014) and the 

02 »>- *«— « 


♦ 


R e 0 IM 0 (p o O # /p e ) 0 * 2 ) 


6. Transition based on \p (see option 5). Transition onset 

determined from built-in curve of V as a function of angle of 

P^cmeters PARA1 and PARA2 are not required input. 
Transition rune length is determined from the ratio, 

T "Fully turb^^trans onset 

1# in P ut in location 028, if desired. If tu> value is 
input, the program will select a value from a built-in table of 
l0 *i* P # 6 vcrcns ratio. 

? * on bnllt in cur >® of onset R _ versus 

Me (Rockwell International criteria). Transition rone length 
aetermined in the same manner as option 6. 

8. Tr *nsition based on R # $ /M . Input requirements arc onset 

value (PARA1, location 014^ and the fully turbulent value (PARA2 
location 020). Nominal values are 150. and /T z 150 respectively. 


Transition Data 
PARA1 


PARA2 


028 


Transition Zone Length Ratio 


(10) ffiAI TRANSFER OPTIONS 
Font at 5(13, F10.6) 


location option or data 

FORMAT (13) FORMAT (F10.6) 


Heat Tra&sfar Method Option Flag (HIFLAG) 


1. Fay-Riddell Stagnation Point Method. 

Required inpot data are nose radiua R,. (location 012) and 

optio,u and 1® with shock and pressure angles 
of SO degrees. An option available (location 313) will 
allow dissociation to be considered (Lewis no. ■ 1.4) or 
not considered (Lewis no. -1.0). 


2. Cato/Johnson Swept Cylinder Method. 

■Ujiirad Inpot dlti Ira local cylinder radiua l, (location 




tion 38 


012), sweep angle A (location 017), flowfield < 

with a shock angle of 90 degrees, and flowfield local 
pressure option 18 with a pressure angle of (90-A) degrees. 


3. Eckert's Reference Enthalpy Flat Plate Method. 

Required input data are the running length (location 013) 
and Mangier transformation factors for both laainar and 
turbulent flow (locations 015 and 016, respectively). 


4. Eckert/ Spa 1 ding-Chi Flat Plate Method. 

Required input data are running length, L (location 013) 
and Mangier transforaation factors for both laainar and 
turbulent flow (location 015 and 016, respectively). If 
the Von Raman Reynold's Analogy factor is desired, 

input a 1.0 in location 319 otherwise it uses Colburn's 
equation. 


5. Boeing Flat Plate Method. 

Required input is the running length L (location 013). 


Modified Beckwith/Gallagher Swept Cylinder Method. 

Required input data are local cylinder radius R* (location 
012) and flowfield option 39 with a shock anglTof (90-A). 


7. Boeing p r g r Swept Cylinder Mpthod. 

Required input data are looal cylinder radius R*, 

012) * nd * lowfi#ld option 39 with a sSock angle 

Of F(rX)e 


8. Lees, Detra and Bidalgo Hemisphere Method. 

Required input data hemisphere nose radius R*. (location 
012), running length L (looation 013), loca?V>dy slope 8 
(location 017), flowfield option 38 with a shock angle of 

90 degrees and flowfield option 18 with a local slope angle 
of o* 
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9. *J«tin and Goodrich, Lee side Orbitcr Ho. ting. 

• F *^ •®J 1 Rldde11 *“P«t (8« - 1 foot for full .c.lo ) 

•nd wiudward wall enthalpy for t* t •caia) 

«»P»«. f(J» . 4,0 BID/lbm moJ.riai) (21> ■ °- # 

I0 ‘ S2r}‘ II* M *P Boating. 

**• *“• •• {0T l»*t transfer option 4 plot Flop 
Ull h j! lnput Jf W(22, ‘ Running length it to the hinge 

4i:«oii ::rnT * a4 re#tucwt -• - oiUon 

11. Fivol, Fin-Plate Peak Interference Heating. Input ia the 

*hock*de t a che ° 4, - Whlch h the d ®*«»lt option if the fin 
TI? pltt * •P eclal ▼•lttea for thie option 

IF W(30) * 0 the fin angle is used, IF W(30) » 1 0 the 

thJ 1 J.J f .iIu Ck B 4 *i* dd ^ J° the fl “ #ngle ' * nd IP f < 30 > - 2.0 

ne yaw angle, 0, is added to the fin ansle The Im.+i. 

»•*»» 5 

loJ% Tl in f( ! 5) \ Th * fln #n * le *«l»tion to the 
local flow at a • fi - o is input in f( 26 ). 


This option is used with the Fay-Riddell ataanatiee ***** 
equation (HIFLAG - 1.0, locational). * 

0. Dissociation considered (Lewis nuaber - 1 . 4 ) 

1. Dissociation not considered (Lewi, nuaber - 1.0) 


(id smm mxmjm 

Format 5(13, F10.6) 


LOCATION 
FORMAT (13) 


018 

019 

320 

321-330 

331-340 

341-350 

360 

361-370 

371-380 

381-390 


OPTION OR DATA 
FORMAT (F10.6) 


Multiplication Factor Input 

tranafcr multiplication factors can ba input aa a 
constant Taint, a function of trajectory tine and/or a 
function of upatroan Mach number. All three input options 
may be used, if desired. All multiplier values are 
multiplied together to obtain an overall value. 

Constant Multiplier (Laminar) 

Constant Multiplier (Turbulent) 

Number of Time Dependent Table Entries (Limit 10) 

Time Table, see 

Multiplier Table (Laminar) 

Multiplier Table (Turbulent) 

Number of Mach Number Dependent Table Entries (Limit 10) 
Mach Number Table (Log 19 M # ) 

Multiplier Table (Laminar) (Log-, h./h ), 

1 U LAN 

Multiplier Table (Turbulent) (Log a9 VVturb 
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(12) nANSFORHATlPHS 

Foraat 5(13, F10.6) 


LOCATION OPTION OR DATA 

FORMAT (13) FORMAT (F10.6) 

319 Reynold's Analogy Factor 

0. As taffies Colburn's Reynold's Analogy factor defined as 

8 “ (P^) ^ (Generally used for design) 

1. Assumes Von Raman Reynold's Analogy factor defined by 
S - 1 ♦ 5(C f )*' t(P r - 1) ♦ In ([5P r + 13/6)3 

(Generally nsed for prediction) 

025 Virtual Origin Option (VRFLG) 

0. Gaea the geonetrie running length in the heat transfer equations. 

1. Corrects the geoaietrie running length to account for the onset of 
transition. Uses the corrected running length in the turbulent 
heat transfer equations. 



Mangier Transforaation Factors 

015 

Laminar Factor “ 3 Cone 

016 

Turbulent Factor ■ 2 Cone 

(13) SURFACE CONDITION 

Foraat 5(13, 

F10.6) 

LOCATION 
FORMAT (13) 

OPTION OR DATA 
FORMAT (F10.6) 

021 

Windward Bottom Centerline fall Enthalpy (Btu/lba). 

Used with heat transfer option (9) for orbiter leeside heating 


Material Eaissiwity 

023 

Eaissiwity (must be > 0.0) 
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(14) GEOMETRY DATA 
Foratt 5(13, F10.6) 




LOCATION 
FORMAT (13) 

OPTION OR DATA 
FORMAT (F10.6) 


Geoaetry Data Input (Constant) 

012 

Note or Leading Edge Radisa, ft 

013 

Gee 'ric Running Length, ft 

017 

Local Slope or Sweep Angle, deg 


Geoaetry Data Input (Funetion of Time) 

550 

Number of Table Entriea (Limit 10) 

561-570 

Time, aec 

571-580 

Note Radii, ft 

581-590 

Running Length, ft 

591-600 

Local Slope or Sweep Angle, deg 

601-610 

View Factor (not uaed) 

354 

Surface Distance to start of turbulent boundary layer. The 
value input in location 354 ia subtracted from the total 
running length input in location 013 to provide the new 
running length for turbulent heating. 

(15) INITIAL CONDITIONS 

Foratt 5(13, 

F10.6) 

LOCATION 
FORMAT (13) 

OPTION OR DATA 
FORMAT (F10.6) 

024 

Initial Temperature, *F 

316 

Integrated Convective Beat Flux, Btu/sq ft 
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<ltf> fQNTROL PARAMETERS 

Format 5(13. F10.6) 


LOCATION 
FORMAT (13) 


209 


314 


<11 

641 


642 


OPTION OR DATA 
FORMAT (F10.6) 


Continuity Option 

0. Convective Heating Load initialled, 

1. Convective Heating Load not initialed. 


Rarefied Geometry Option (NFCS) 

0. Sharp cone rarefied equation! used 

1, Flat plate rarefied equation! uied 


Body Point Number 


Program Input Control Parameter (JFI) 




• d,t * u#in « title, timing parameter! 

•nd trajectory data from previou. caae. 

<D all unchanged caae data from previous case.) 


2. End of input data (use with last case). 


3 ’ 22 f roa/previoua “T 

(B«e! all unchanged case data from previous caae.) 


4. title, timing, trajectory and caae data 

(Initially zeros full f array) 


3. Same a. l (except zero, all c ... data from previou, caae). 


6. Sa« a. 3 (ezc.pt zero. .11 timing parameter, and c... 


Print Option (JTFQQ) _ 

0. Normal Miniver printout 

1. Normal Miniver printout plus VANODT (aummary) printout 

2. VAN<>DT printout only 


18 


. .... • 


' j 


i 


■ I 

i 

• ! 


! 


i 


0. No output fil. fro* VANOUT created 

1. OUTPUT FILE fro* VANOUT created 
Rarefied Flow Option (NONCON) 

0. Conti nou* flow equations are used throughout flight 

1» Continuum and rarefied flow equations are used* 
depending on flow regime encountered during flight. 

Output Units Option 

0. English 

1. Metric 

Pressure Coefficient Input Option (JFCP) 

These cards follow 541 specifications 
Card 1 NCPMP (13) 

Card 2 M^, Cp (2F10.6) 

Card 3 M«o» Cp (2F10.6) 

0. Pressure coefficients are not input 

1. Pressure coeffi cients are input (Flowfield option 14). 
Cone Flow Option (NSB) 

This option is used only when rarefied flow conditions 
are considered (Control Parameter 646.) 

0. Uses sharp oone correlations 

1. Uses blunt cone correlations 
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Section 3.0 

PREPROCESSOR INTERACTIVE INPUT 


HI. .action coot. in. de. crip, ion. of ..hr..,!... bjr 

...... ..d. premin. premin i. .. pr.pro.....r .... *„ lanmin. 

"•*••• * ‘ lU th,t “ *• «• •— «rd 1.... for... .. di..„...d i. s..- 
tl..J.0.nd i. »..d .. th. input til. for LANMIN. Th. pri.ery 10 fo,«„ lo „ u _ 
.loded I, thi. fil. !. ... A definition of ot 

'W* array ia given in Appendix A. 


T». PREMIN program .on.!.,. 0 f 18 rontln.. ,nd on. ..in progra. n.d 1. ... 

pnbl. of .renting . new ontpnt fil. or .odifpin, .id one. », INI>0T 

i. n..d to Input th. fil. ,o „. ..difi.d „d ,h. rontin. OUTPUT 1. ...d to 

.rente th. ontpnt file. A .aero flow dl.gru of PREMIN i. included in Fig... 
3-1. 


3.1 MAIN 

MAIN performs the task of managing the rest of PREMIN. The dr ceri-.inatioi 
of whether PREMIN is to crests a new ontpnt file or modify an rid ontpnt fil. i, 
«*de within MAIN. After the method is determined, MAIN calls the required 
rontinea in their proper seqnence a case at a time. The variables transferred 
in Common or nsed as control flags by MAIN are: 

!£* ‘ J Bter#ctl ™ Ia P«t - FORTRAN Unit Nnmber 

MW Aft-n " I 0ttt P wt - FORTRAN Unit Nnmber 

MFLAG-0 - English units 

EFLAG«1 - Metric units 

FLAG a O - Create new output file 

H*Afl“l - Modify old output file 

JFLAG*1 - Previous case deleted 

IN SEEM - Current esse inserted 

WI^6"1 - TITL1 currently defined 

FNAM( - Name of input 

JFKS - Program Input Control Parameter (W(641)) for the 
previous esse 
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NC - Casa nimbus 

ND - Number ciki deleted 

NI " Number ceeee inserted 

THtl - Title of i set of ceeee 

NCPHT - Number of Mech numbers in CP vs. Meoh number table 

NCPMT8 - Saved value of NCPMT from previous esse 

TMCP - Mach number array for CP vs. Mach number table 

TMCPS - Saved values of TMCP array from previous case 

TCPM - CP array for CP vs, Mach number table 

TCPMS - Saved values of TCPM array from previous case 

V - Array of values created and modified by PR EM IN 

(always in English units) (see Appendix A for complete 
description) 

WW - Find array copied to output file. Copy of 'f' array 

except changed to desired output units. 


3.2 SUBROUTINE INPUT 

Subroutine INPUT is called by MAIN and performs two functions. If IB which 
is the parameter passed in the subroutine call equals 2 then INPUT xeros out the 
values in the f array according to the Program Input Control Parameter of the 
last ease and returns to MAIN. 

IF IB-1 then INPUT reads an existing LANMIN input file for the purpose of 
PREMIN cresting s new LANMIN input file. INPUT first xeros out the 1 array ac- 
cording to the Program Input Control Parameter of the last case and then reads 
in data for a single case from the input file. INPUT keeps track of the current 
case number and inserts and deletes cases. 


3.3 SUBROUTINE MODIFY 

Subroutine MODIFY is called by MAIN to modify the V array for the current 
case. MODIFY calls a requested imroutlne in order to modify or redefine a 
specific section of the W array. MODIFY is also capable of changing the value 
of a specific V number. 

3.4 SUBROUTINE TIMING 

Subroutine TIMING is called by MAIN or MODIFY to define the Timing Parame- 


tar* and Print Control values. TIMING assigns valnea for W(l) through W(8) and 
checks to aaa If the number of print times exceeds 100. SITM equals the total 
number of print timas. 

3.5 SUBROUTINE THAT 

Subroutine TRAJ is called by MAIN or MODIFY to define the trajectory data 
to be used. TRAJ obtains this data from a user supplied trajectory file or in- 
teractively via the terminal. 

If the trajectory is read from a trajectory file, TRAJ first determines 
whether Beta data is included and the number of time dependent table entries. 
TRAJ then reads ALT, VEL, Angle of Attack, and YAW Angle (if included) for each 
TIME. After obtaining the trajectory data from the trajectory file, TRAJ will 
display the trajectory data on the terminal 20 lines at a time waiting for a 
carriage return between pages. After display TRAJ gives the user the opportuni- 
ty to modify th© trajectory data. 

The format for the trajectory file is an A80 for the title on line one. 
Line two contains the Beta flag using an 12 and the number of trajectory times 
using an F10.5. The rest of the file contains the trajectory data using either 
a 4E15.4 or a 5E15.4, depending on whether or not the yaw angle Jat* is includ- 
ed. 

The other method of creating the trajectory data is via the terminal by in- 
putting a line at a time. After defining the trajectory in this manner, TRAJ 
proeeeds to display the trajectory data on tt terminal as described above. 

After the trajectory data table is displayed on the terminal, the tuisr la 
given the opportunity to modify part or all of the trajectory. TRAJ enables the 
user to insert, delete, or change a line of the trajectory. After modification, 
the trajectory is displayed on the terminal as before. 
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Once the user le satisfied that the trajectory is complete TRAJ will arable 

the uaer to create a new trajectory file with a new user supplied file name. 

The wariablaa need in TRAJ are: 

FNAM1 - The filenaae of the trajector file. 

LINE - The number of a line at the trajectory table. 

TITLE — Title of the trajectory file. 

BFLAG*1 - YAW Angle incloded with trajectory. 

3.6 ‘SUBROUTINE ATMS 

Subroutine AIMS ia called by MAIN or MODIFY to aelect the Atmosphere option 
to be need. 

3.7 SU r IOUTINE WNDTUN 

Subroutine WNDTUN ia called by MAIN or MODIF” if the Ataoaphere Option 2 

was selected. WNDTUN defines the Atmospheric data for a Wind Tunnel case. 

WNDTUN sets values for 1^, and for each tine point of the trajectory NPTS is 
the number of time points. 

3.8 SUBROUTINE ATMDTA 

Subroutine ATMDTA is called by MAIN or MODIFY if the Atmosphere Option 3 
was selected. ATMDTA is the routine the user uses to define his own set of At- 
mospheric data. ATMDTA asks for the number of Altitude entries (NALT) and then 
asks for the values of Altitude. T w .and P„ for each entry. 

3.9 SUBROUTINE HEAT IN 

Subroutine HEAT IN is called by MAIN or MODIFY if a Beating Indicator case 

is desired. HEAT IN sets the array values for a 1 ft. radius sphere with a 

wall temperature of 0 F and a Lewis number of 1.0. BEATIN uses the Heat 
transfer option 1 which is the Fay and Riddell Stagnation Point Method. BEATIN 
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the ODUanfi Shod Plowf j«, opt 


ion with a fchocik angle of 90 . The modified 


" «‘'uu j eu 

Newtonian Praasnre option ia oaed with ♦»,. , * 

uaed with th. body angle of 90 . HEATJN al.o 

tk. PM.t "Ml «n to 2 which ,h. ontpnt for Slu „„ rv ^ 

3.10 SUBROUTINE HXRNTI) 

S -kr..U.. HTKMTII I. celled by 

u, '« «... mm, .... a. r. r p.,... t . r . b> u . 

.... Traveler Opt loo. After th 1 ... .re HIUHID .... .... 

p Tttnt, w.ll cmi e e iv i ty. end continuation option, HTOMTD Jot. tbo oter «p. c i- 

O .bother b. Rarefied flow ood virtool or,,,, optio... H1HMTD .boob. lt 

the rarefied flow option waa eelected. Rarefied r, 

0# Rarefied flow cannot be need by Beat 

Transfer Options 2, 3, 5. 7, 9, 10, or 11. 

3.11 SUBROUTINE HTMULT 

Sobrontiao BTMULT i. C.U.d bp MAIM or MODIFT i f . 

n Factor it desired. There tre 3 option type. end BTMULT allow, the neer 

“ a ’ Cld ' " “ ,lrtM *° * « »>»«■ 1 i. cboaen. . Le.i.tr and 

Tnrbnlont factor .net b. .applied. If 2 „ p( 

«.« b. aapplied. For eaeb ti.. . Laminar and tnrbnl.n, .... 

Plied. If option - i. eelected. ,». nnmber nn. b „. „„ 

For each Ma.b nnmber. a Laminar and t.rbnl.n, ...tip,!., u 

■TMULT convert a ,b... 1 .., multiplier. t . u , p „ ^ |# 

lb.ro net be a mlnimnm of 2 and a maxima, of l 0 entriea. 

3.12 SUBROUTINE TRANS 

Subroutine TEAMS i. called b, MAIM or MODIFY to , h . pp . 

•leu to be need. HUNS the. r.,.ira. ,h. bi,b.e, .„d ,o..., 




th» transition region depending on the option seleoted. 


3.13 SUBROUTINE CROSS 

Subroutine CROSS ie celled by MAIN or MODIFY if a Cross Flow adjnstaent is 
needed. (BOSS selects the Cross Flow adjnstaent option to be used and then asks 
for data depending on which option is selected. 

3.14 SUBROUTINE FLOW 

Snbrontine FLOf is called by MAIN or MODIFY to set the flowfield and local 
pressure options. The Paraaeter ID is passed through the snbrontine call. If 
ID-1 then FLOW will create a new set of flowfield data. If ID-2 then FLOW aodi- 
fies an old set of flowfield data. If ID-2 then the flowfield data anst be 
transferred froa the f array to the flowfield and local pressure arrays. FF is 
the flowfield option array. FFA is the array containing the angles correspond- 
ing to the flowfield option. P is the local pressure option array. PA is the 
array containing the angles corresponding to the local pressure option. Once 
converted, FLOW echos the psirs of flowfield and local pressnre options with op- 
tion abrevlations and corresponding angles. If no changes are desired, then 
FLOW is exited. Otherwise, FLOW provides a table of flowfield and local pres- 
sure options with option nuabers. The user aust select options by pairs with 
corresponding angles. Choosing option -1. for the flowfield option and the 
local pressure option signifies that no aore options are to be included. If the 
CP vs. Mach nuaber table pressure option is chosen, FLOW requires a filenaae 
containing the table, or the nuaber of table entrees, and M w and CP for each 
•ntry. The foraat for the CP vs. Maoh nuaber file is an 13 on line one for the 
nuaber of Mach nuabers. The reaaining lines contain the Mach nuaber and the 
pressure coefficient using a ZF10.6 foraat. After all options are selected. 
nm echos the pairs of options. If the options are correot, the flowfield and 
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Joc.l pressure err.ys ere converted to the W array. FNAM5 is the variable for 
the filename of the pressure coefficient vs Mach number table. 

3.15 SUBROUTINE TDGEOM 

Subroutine TDOEOM i. called by MAIN or MODIFY if time dependent Geometry is 
needed. First, the number of time dependent entries is set. Then the time, ra- 
dius, length, and slope or sweep angle are required for each entry. 

3.16 SUBROUTINE 00N1RL 

Subroutine CONTRL is called by MAIN or MODIFY to set the LANMIN control 
flags. CONTRL sets the body point number for the case, the print option, output 
nnits option, and the program input control parameters for the case. CONTRL 
•«ks if LANMIN will create an output file and if the current case is the first 

of a set of streamline cases. If the case i. the first of a set of streamline 
cases subroutine STREAM is called. 

3.17 SUBROUTINE UNITS 

Subroutine UNITS 1. called by MAIN to copy tb. .orUn, « array 

output WW array. If HFLAG»1. than metric unite arc required and UNITS convert, 
the WW array into Metric units. 

3.18 SUBROUTINE STREAM 

Subroutine STREAM is called by CONTRL if a aet of streamline case, i. to be 
created. If a streamline file is created by STREAM, it will be used by OUTPUT. 

STREAM first checks if the time dependent geometry option was selected. If 
it was, then streamline cases can not be used and STREAM returns control to 
CONTRL. If time dependent geometry was not used, then STREAM opens a temporary 
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file to store the date for the streamline cases. Neat, STREAM determines the 
shock and local prestnre options to be nsed by the streamline cases from the 
previous case. STREAM then requires X distance, shock angle, body angle, and 
body point number for each streamline case. STREAM writes these values and the 
’W’ number that these values are to be stored in to the streamline file. The 
parallel shock and Prandt-Meyer pressure options don't require both the shock 
and body angles to be input. A negative value of X indicates the end of the 
streamline cases. 

The variables used by STREAM are: 

X - X distance 

SA - Shock angle 

BA - Body angle 

BP - Body point number 

CASE - Case number 

TPRT — Program input control parameter 

• * - 'W* index for CASE 

J2 - 'W' index for X 

J3 - index for SA 

J4 ~ 'W' index for BA 

J5 - *W' index for BP 

J€ - *W* index for XPRT 

“ Flag used to mark parallel shock or Prandt— Meyer cases. 
3.19 SUBROUTINE OUTPUT 

Subroutine OUTPUT is called by MAIN to write the *f‘ array values for a 
case to the output file. FNAM4 is the variable for the filename of the output 
file. OUTPUT opens this file unless previously opened. Next OUTPUT writes the 
data for the eurrent case to the output file according to the program input con- 
trol parameter of the previous case. OUTPUT writes to the file the title, tim- 
ing parameters, and print control, the nomber of time dependent trajectory table 
entries, the trajectory data, and the case data. The last I of the case data to 
be written to the file is the program input oontrol parameter. OUTPUT checks to 
•ee if the case Just written is the first of a set of streamline cases. ISIRM»1 
indicates streamline cases. If streamline cases are indicated, then OUTPUT co- 
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Piet the streamline ceaee to the OUTPUT file from the temporary streamline file. 
After the last streamline caae is written to the OUTPUT file, OUTPUT sets the 
program input control parameter for the last streamline ease. OUTPUT then 
closes and deletes the temporary streamline file and returns to MAIN. 


NIVW 



PREMIN Macro Flow Diagram 




















Section 4.0 


OUTPUT DESCRIPTION 


The printed ontpnt fro* LANMIN comet in two forns. The detailed printout 
for each tlae point la given if JFQQ - 0 or 1. Thie detailed printout containa 
-,h<§ vaiiables Hated in Table 4.1* The definitioni and units are given in Table 
4.1 for each of the variablea printed. 

If the suaaary printout it eelected (i.e. JFQQ - 1 or 2) the output con- 
taina the variablea shown in Table 4.2. This table containa a listing of twelve 
ti*e dependent variables for each body point. The variable and units are clear- 
ly labeled on the output. The euaaary page ean be obtained for two sets of out- 
put units. If 1(448) - 0 the output is in Englieh units and if 1(448) » 1 the 
output ia in Metric units as illustrated in Table 4.2. 

An output, tape or file, ean be created fro* LANMIN if 1(443) » 1.0. The 
output file is a card inage of the printed suaaary page. 
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Table 4.1 


INTERMEDIATE PRINT OUTPUT PARAMETERS 

VARIABLE DEFINITION UNITS 


TIME 

TIME 

ZINF 

ALTITUDE 

VINF 

FREESTREAM VELOCITY 

A INF 

FREESTREAM SPEED OF 

VE 

EDGE VELOCITY 

RE INF 

FREESTREAM REYNOLDS 

REL 

FREESTREAM REYNOLDS 

MINF 

FREESTREAM MACH NO. 

MU 

UPSTREAM MACH NO. 

MB 

EDGE MACH NO. 

L 

RUNNING LENGTH 

RN 

NOSE RADIUS 

PR 

PRANDTL NUMBER 

EMS*F 

EMISSIVITT * VIEW Fi 



(SBC) 


(FT) 


(FT/SEC) 

SOUND 

(FT/SEC) 


(FT/SBC) 

NO. /FT 

(1/FT) 

NO. 

( ~ ) 


( — ) 
( — ) 
( — ) 
(FT) 
(FT) 

( — ) 
( — ) 


PINF 

FREESTREAM PRESSURE 

PU 

UPSTREAM PRESSURE 

PB 

EDGE PRESSURE 

FT 

POST SHOCK TOTAL PRESSURE 

ALPHA 

LOCAL ANGLE OF ATTACK 

CF/2 

LOCAL SKIN FRICTION C0EF./2 

TAD V 

WALL SHEAR FORCE 


(LBF/SFT) 

(LBF/SFT) 

(LBF/SFT) 

(LBF/SFT) 

(DEG) 

( — ) 
(LBF/SFT) 


TINF FREES 1REAM TEMPERATURE 

TO UPSTREAM TEMPERATURE 

IE EDGE TEMPERATURE 

IT POST SHOCK TOTAL TEMPERATURE 

T* ECKERT REFERENCE TEMPERATURE 

TR p|t REFERENCE TEMPERATURE 

TV WALL TEMPERATURE 


(F) 

(F) 

(F) 

(F) 

(F) 

(F) 

(F) 


1 INF FREESTREAM STATIC ENTHALPY 

HU UPSTREAM STATIC ENTHALPY 

HE EDGE STATIC ENTHALPY 

IT TOTAL ENTHALPY 

H* ECKERT REFERENCE ENTHALPY 

HR pit REFERENCE ENTHALPY 

HW WALL ENTHALPY 

PHI SWEEP ANGLE 


(BTU/LBM) 

(BTU/LBM) 

(BTU/LBM) 

(BTU/LBM) 

(BTU/LBM) 

(BTU/LBM) 

(BTU/LBM) 

(DBG) 


RBOI 

FREESTREAM DENSITY 

RBOU 

UPSTREAM DENSITY 

RBOB 

EDGE DENSITY 

RIOT 

POST SHOCK TOTAL DENSITY 

RHO* 

ECKERT REFERENCE DENSITY 

RHOR 

pp REFERENCE DENSITY 

RHOV 

WALL DENSITY 


(SLUG/CFT) 

(SLUG/CFT) 

(SLUG/CFT) 

(SLUG/CFT) 

(SLUG/CFT) 

(SLUG/CFT) 

(SLUG/CFT) 
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PARAI 

MO INF 
MO 0 
MO B 
MO T 

MO • 

MO R 
MO V 

PARA2 

FF FLAG 
ANGLE 

NC L 
NC T 

HRECOV L 
HRECOV T 

QC L 
QC T 

KSOB 12 
KSUB 12 

KSOB L3 
KSOB T3 

KSOB L 
KSOB T 

PARA 

PCTT 

NSOB C 
1 IDEAL 

H RECOV 
T RECOV 

QCONV 

QCTOTAL 

QRAD 
OR TOT 

Q NET 
ON TOT 

TRAD EQ 

QC Cff 
QC CWT 


LAMINAR TRANSITION PARAMETER 

FREESTREAM VISCOSITT 
UPSTREAM VISCOSITT 
EDGE VISCOSITT 
POSTSBOCE TOTAL VISCOSITT 
ECKERT REFERENCE VISCOSITT 
p|t REFERENCE VISCOSITT 
WALL VISCOSITT 

TURBULENT TRANSITION PARAMETER 

FLOWFIELD FLAG NUMBER 
ANGLE FOR CORRESPONDING FLAG 

LAMINAR HEAT TRANSFER COEFFICIENT 
TURBULENT HEAT TRANSFER COEFFICIENT 

LAMINAR RECOVER! ENTHALPY 
TURBULENT RECOVERY ENTHALPY 

LAMINAR HEATING RATE 
TURBULENT HEATING RATH 

LAMINAR TIME DEPENDENT MULTIPLIER 
TURBULENT TIME DEPENDENT MULTIPLIER 

LAMINAR MACH NO. DEPENDENT MULTIPLIER 
TURBULENT MACH NO. DEPENDENT MULTIPLIER 

TOTAL LAMINAR MULTIPLIER 
TOTAL TURBULENT MULTIPLIER 

TRANSITION PARAMETER VALUE 
PERCENT TURBULENT HEATING/ 100 

HEAT TRANSFER COEFFICIENT 

IDEAL GAS AIR HEAT TRANSFER COEFFICIENT 

RECOVERY ENTHALPY 
RECOVERY TEMPERATURE 

HEATING RATE BASED ON HI 
HEATING LOAD FOR TV 

WALL RADIATION RATE AT TV 
WALL RADIATION LOAD AT TV 

NET HEATING RATE (QCONV-QR) 

NET HEATING LOAD 

RADIATION EQUILIBRIUM WALL TEMPERATURE 

COLD WALL HEATING RATE (TW - 0) 

COLD WALL HEATING LOAD 

3f 


( — ) 

(LBF-SEC/SFT) 

(LBF-SEC/SFT) 

(LBF-SEC/SFT) 

(LBF-SEC/SFT) 

(LBF-SEC/SFT) 

(LBF-SEC/SFT) 

(LBF-SEC/SFT) 

( — ) 

( ~ ) 

(DBG) 

(L8M/SFT-SEC) 

(LBM/SFT-SEC) 

(BTU/LBM) 

(BTU/LBM) 

(BTU/SFT-SEC) 

(BTU/SFT-SEC) 

( ~ ) 

( ~ ) 

( ~ ) 

( — ) 

( — ) 

( ~ ) 

( ~ ) 

( — ) 

(LEM/SFT-SEC) 

(BTU/HR-SFT-R) 

(BTU/LBM) 

(F) 

(BTC/SFi-»EC) 

(BTU/SFT) 

(BTU/SFT-SEC) 

(BTU/SFT) 

(BTU/SFT-SEC) 

(BTU/SFT) 

(F) 

(BTU/SFT-SEC) 

(BTU/SFT) 


ORIGINAL PACE Eu> 
OF POOR QUALITY 
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APPENDIX A 


V ARRAY DESCRIPTION 
ARRAY / SYMBOL/ VAR IABLE/DEFINITION 
1(1) » T1 ■ tl (first calculation tiae, tias 1) 

1(2) ■ DTI ■ Atl (delta tlao step between tiae tl end t2) 

1(3) » 12 - t2 (tiae 2) 

1(4) ■ DT2 ■ At2 (delts tiae step between t2 end t3) 

1(5) - T3 - t3 (tiae 3) 

1(d) ■ DT3 “ At3 (delts tiae step between t3 end t4) 

1(7) » T4 • t4 (tiae 4) 

1(8) ■ DTCALC (KCALC eslenlstion interval factor) 

1(9) ■ (Open Location) 

1(10) ■ ATFLAG (ataospberic freestreaa properties option flap) 

1(11) * HTFLAG (NHFLAG - beat transfer aetbod option flag) 

1(12) ■ RN ■ R N (nose radius or loesl cylinder radios, leading edge 
radios, ft,) 

1(13) ■ EL (ronning length, ft.) 

1(14) ■ PARA1 (transition data) 

1(15) • ENL (2nd power of Mangier transformation factor for laainar flow) 
1(16) ■ ENT (fifth power of Mangier transforaation factor for torbnlent flow) 
1(17) • Pf-r (sweep angle. A, local slope, 6, deg.) 

1(18) ■ AXL ■ (constant laainar aoltiplication factors) 

1(19) ■ ACT ■ (constant torbolent aoltiplication factors) 

1(20) ■ PARA2 (trsnsition dsta) 

1(21) ■ lindward Bottoa Centerline lall Enthalpy (Bto/lba) 

1(22) - Flap Length (ft) for heat transfer option 10 
1(23) - EMIS (eaissivity > 0): If (EMIS.LB.0)EMIS»0.8 

1(24) ■ TIN (initisl teaperatore, *F) 


A-l 


W(25) - 
W (26) « 
W(27) - 
V(28) 
1(29) 
1(30) 

1(31) 
1(32) 
1(33) 
1(34) 
1(35) 
1(36) 
1(37) 
1(38) 
1(39) 
1(40) 
1(41) ■ 
1(42) ■ 
1(43) - 
1(44) ■ 
1(45) - 
1(46) - 
1(47) ■ 
1(48) « 
1(49) • 


Pi i Length to Point of interest for Heet Transfer Option 11 (ft) 
Pin Angle a » fi - 0 for Beat Transfer Option 11 (Deg) 

TRFLAG (transition criteria option flag) 

ELPAC (tramitlon tone length ratio) 

VKLFLG (virtual origin option) 


Pin Option Flag 0 - Pin angle 

1 ■ a ♦ Pin angle 

2 ■ p ♦ Fin angle 


•> Flowfield and local pressure option flags 


GF(1) » FF(1) 

QF(2) - FF(2) 

GF(3) - FP(3) 

GF(4) - FF(4) 

GF(5) - FF(5) 

GF(6) - FF(6) 

BLFA(l) (ALFA(l) - BLPA(l) + ATAK) degree 
BLFA(2) (ALFA (2) - BLFA(2) + ATAK) degree 
BLFA(3) - ALFA (3) degree (add ATAK to ALFA if 1(263) - 3.) 

BLFA(4) - ALFA (4) degree (add ATAK to ALFA if 1(264) - 4.) 

BLFA(5) - ALFA(5) degree (add ATAK to ALFA if 1(265) - 5.) 

BLFA(6) - ALFA(6) degree (add ATAK to ALFA if 1(266) - 6.) 

BLFA(7) - ALFA (7) degree (add ATAK to ALFA if 1(267) - 7.) 

BLFA(8) - ALFA (8) degree (add ATAK to ALFA if 1(268) - 8.) 

BLFA(9) - ALFA (9) degree (add ATAK to ALFA if 1(269) - 9.) 

HH(1) - FF(7) 

I 

HB(2) - FF(8) > Flowf ield and local pressure option flags 

I 

HH(3) - FF(9) 


(Open Location) 

ENTB (nonber of trajectory points) 


1(50) 


1(51) - 12(1)- 


-> (Trajectory Tine, aeo) 


W(100). 

W(101). 

1(150)* 

1(151)* 


12 ( 30 )- 

22 ( 1 )- 


’> (Trajectory Altitude, feet) 


ZZ(50)° 

VZ(1)— 


*> (Trajectory Velocity, ft/sec) 


1(200)- 12(50) 
1 ( 201 ) 

1 ( 202 ) 

1(203) 

1(204) 

1(205) 

1(205) 


- CFFLG (crottflow optioa flag) 

DSOBO (croeaf low rectangular width, 0 , ft) 

c 

ELNBDA (delta wing aweep angle, X (deg)) 

(real gaa velocity gradient correction factor, O) 

CORNS (rectangle corner radiua, S (ft)) 

0 


1(208)' 

1(209)> 

1 ( 210 )' 

1 ( 211 )' 


I 


-> (Open Locationa) 


CONFLG (continuity option) 

(Open Location) 

Al(l) 

I 

(Trajectory poaitive a, degree) 


1(260)- Al(50) 
1(261) 


I 


1(262)' 


*• Built into o added to 1(37): 

If (1(261). LB.0)1(261)-1.0 

2. | program a added to 1(38) : 

If (1(262). LB.0)l(262)-2.0 


A- 3 


W(263)= 

f(264)* 

W(265)- 

W(266)» 

W(267)» 

W<268)« 

W(269)» 

W(270)« 


3. 

1 

Values 

“eff 

4. | 

1 

am st be 

“eff 

5. 1 

1 

input for 

°eff 

6. ~~> 

1 

t o be 

“eff 

7* 1 

1 

added to 

°eff 

8. | 

1 

body 

%ff 

9. 

angle 

“eff 


•dded to W(39) 
•dded to »(40) 
*<lded to W( 4 J) 
added to W(42) 
Added to W(43) 
added to W(44) 
•dded to W(45) 


W(313)« 

W(314)= 

W(315)- 

W(316)« 

*(31?). 
W(318) = 
W(3l9)-= 

W(320)« 

W(321)« 

• 

W(330)« 
W(331)« 


> (Open Locations) 


I 


NFCS (rarefied flaa 0 - c v a __ 

I *h*rp cone option need 

~ plate option osed) 

NLEW <<lia,o.i. tio „ fi, t o - Lewis rat „ . 1.4 

1 “ Lewis number « 1,0) 

OCSmiA coovective he., 8../., f„ 

(Open Location) 

(Open Location) 

' WLG (Flag trietl ° n J «»•» »..0) 

0 Colburn a equation used) 

‘SST,? "•« *» •• tantM of «i.„ 

1X1(1)- 


I 


IX1(10)~ 
AKL2 (1) — 


W(340)« AKL2(10)- 


I 

| > (Ttble of ®«ltiplier times, sec) 


| > <Table 0f ltoin ar multiplier factors] 


W(34l) 


AKT2(1) 


W(350). 

1(351)« 

• 

K3S3)- 

1(354)- 

W(335)- 

• 

1(359)- 

1(360)- 

1(361)- 

• 

1(370)- 

1(371)- 

• 

1(380)- 

K38D- 

• 

1(390)- 

1(391)- 

• 

1(399)- 
1(400)- 
1(401). 


I 

| > (Tabl * of turbulent multiplier factors) 

AKT2(10) 


9 

j > (Open Locations) 

Ximc (md.. dl.t.ac. «c of tnrbole&t bowdor, l.y.r, f„ 


I 


-> (Open Locations) 


ENM3 (number of Mach dependent multipliers) 
TMACE(l) 


IWACB(IO)- 
AKL3(1) — 


I 

. > »** 10 ) <*») 


AEL3(10)- 
AKT3 (1)— 


AET3 (10)* 


I 

| * *10 *i^n* laminar multiplier 

I 

. > (Lob, . h./h ) 

| 10 i u' turbulent multiplier 


I 


•> (Open Locations) 


ENAL1Z (Table of wind tunnel freestream properties, limit 50) 
FSALT(l)- 


1(450)- FSALT(50) 


"> (Table of altitude data, ft) 


5 


1(451)- FSTEMP(l)- 


1(500)- 

1(501)- 

• 

1(550)- 

1(531)- 

• 

1(559)- 

1(560)- 

1(561)- 

♦ 

• 

1(370)- 

1(571)- 

• 

1(580)- 

1(381)- 

• 

1(590)* 

1(591)- 

• 

1(600)- 

1(601)- 

1(610)- 


FSTEMP(SO)- 
FSPRES(l)— 

FSPRES(50) 


I 


I 


■> (fable of wind tunnel freestream temperature*. 


*> (Table of wind tunnel free stream pressure, 
lb/aq ft) 


•> (Open Locations) 


ENT3 (Number of Table Entries for Geometry as Function of Time) 

Tied) 


tie do- 
med)— 


mc(io)- 

ELZ(l)— 


ELZ(10)* 

FHIZ(l)- 


PHIZ(10)- 


•> (Body Geometry is Varying with Time, See) 


*> (Nose Radius, Ft) 


*> (Running Length, Ft) 


•> (Local Slope or Sweep Angle, Degree) 


*> (Open Locations) 


A- 6 


1(611)- Body Point Number 


1 ( 612 )- 

• 

1 ( 640 )- 

1 ( 641 )- 

1 ( 642 )- 

1 ( 643 )- 

1 ( 644 )- 

1 ( 645 )- 

1 ( 646 )- 

1 ( 647 )- 

1 ( 648 )- 

1 ( 649 )- 

1 ( 650 )- 

1 ( 651 )— 


1 ( 700 ) 


I 

, > (Open Locations) 


JFK (Input Control Flap) 

JFQQ (Print Option) 

Output Fila from VANODT 
(Open Location) 

JFOO (Initialisation Option) 

NON CON (Rarefied Flow Option) 

(Case Number) 

Output Unite Option, O.-Bnglish, 1. -Me trie 
JFCP (Pressure Coefficient Input Option) 
NSB (Cone Flow Option for Rarefied Flow) 

I 

. > (IAW Angle g, degree) 
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SAMPLE CASES 

Two sample oases are presented in this appendix whieh exercise both the 
PREMIN and LANMIN oodes. The components shown here follow those depicted in 
Fi|. 1.1. 

flmiPg IaAl£lt£X Sample SbIS. 

The inpnt for a heating indicator run is rather straight-forward using PRE- 
MIN. After the execution command, units selection, and printout selection, the 
trajectory is input. PREMIN then displays the trajectory data twenty lines at a 
time waiting for a carriage return before displaying the next twenty. After the 
trajectory is input, the atmospheric option is chosen. Hie preceeding steps are 
standard in most input. If the heating indicator option is selected, the re- 
maining input is much simpler than for most other options. During execution of 
PREMIN, if the answer to a question is an end of file character, then PREMIN 
will stop execution. If the answer to a question is incompatible with the read- 
ing format, then the question will be reasked until an acceptable answer is re- 
ceived. 

The PREMIN interactive input is followed by a listing of the input trajec- 
tory file. The first line of the trajectory file contains the title. The 
second line contains a Beta flag and the number of trajectory table entries. 
The third line tbroogh the end of the file contains the trajectory data. Each 
line of trajectory data contains a time, altitude, velocity, angle of attack, 
and, if the Beta flag equals one, the yaw angle, A listing of the PREMIN output 
file follows the trajectory file. The LANMIN output listing is given last. The 
susimery output is labeled with units and is self explanatory. 
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TRAJECTORY INPUT 



DO VOU HAVE A TRAJECTORY INPUT FILE ? 

vcc 

ORIGINAL PAGi; 53 

WHAT IS THE FILE NANE ? 


OF POOR QUALITY 

REENTRY. TRJ 





TIME 

ALTITUDE 

VELOCITY 

ANGLE ATTACK 


(SEC) 

(FT) 

(FT/SEC) 

(DEG) 

1 

0.3000E+00 

0.3963E+06 

O.2457E+05 

0.41 13E+02 

2 

0.4530E+02 

0.3738E+06 

0.2459E+05 

0.4126E+02 

3 

0.9030E+02 

0.351 5E+06 

0.2462E+05 

0.4121E+02 

4 

0.1353E+03 

0.3295E+06 

0.2464E+05 

0.4046E+02 

5 

0.1803E+03 

0.3080E+06 

0.2466E+05 

0.4065E+02 

6 

0.2253E+03 

0.2880E+06 

0.2466E+05 

0.41 90E+02 

7 

0.2703E+03 

0.2697E+06 

0.246 1E+05 

0.393 9E+02 

8 

0.3153E+03 

0.2563E+06 

0.2447E+05 

0.407 1E+02 

9 

0.3603E+03 

0.2502E+06 

0.2422E+05 

0.4174E+02 

10 

0.4053E+03 

0.2470E+06 

0.2392E+05 

0.3995E+02 

11 

0.4503E+03 

0.2446E+06 

0.2361 E+05 

0.3928E+02 

12 

0.4953E+03 

0.2425E+06 

0.2328E+05 

0.3966E+02 

13 

0.5453E+03 

0.2400E+06 

0.2288E+05 

0.3907E+02 

14 

0.5693E+03 

0.2388E+06 

0.2268E+05 

0.3932E+02 

15 

0.5933E+03 

0.2377E+06 

0.2247 E+05 

0.3902E+02 

16 

0.6173E+03 

0.2366E+06 

0.2225E+05 

0.3925E+02 

17 

0.641 3E+03 

0.2353E+06 

0.2202E+05 

0.3948E+02 

18 

0.6653E+03 

0.2337E+06 

0.2178E+05 

0.5997E+02 

19 

0.6893E+03 

0.2320E+06 

0.2153E+05 

0.4001 E+02 

20 

0.7133E+03 

0.2323E+06 

0.21 26E+05 

0.4025E+02 


TIME 

ALTITUDE 

VELOCITY 

ANGLE ATTACK 


(SEC) 

(FT) 

(FT/SEC) 

(DEG) 

21 

0.7373E+03 

0.2303E+06 

0.2098E+05 

0.4046E+02 

22 

0.7613E+03 

0.2279E+06 

0.2068E+05 

0.401 4E+02 

23 

0.7853E+03 

0.2252E+06 

0.2036E+05 

0.401 6E+02 

24 

0.8093E+03 

0.222 3E+06 

0.2002E+05 

0.401 6E+02 

25 

0.8333E+03 

0.21 90E+06 

0.1 965E+05 

0.4027E+02 

26 

0.8573E+03 

0.2153E+06 

0.1924E+05 

0.4040E+02 

27 

0.881 3E+03 

0.21 14E+06 

0.1879E+05 

0.401 2E+02 

28 

0.9053E+03 

0.2076E+06 

0.1830E+05 

0.4030E+02 

29 

0.9293E+03 

0.2056E+06 

0.1776E+05 

0.4001 E+02 

30 

0.9533E+03 

0.2025E+06 

0.1718E+05 

0.4200E+02 

31 

0.9773E+03 

0.1 972E+06 

0.1654E+05 

0.4090E+02 

32 

0.1001E+04 

0,1 922E+06 

0.1587E+05 

0.407 5E+02 

33 

0.1 025E+04 

0,1 872E+06 

0.15I3E+05 

0.3987E+02 

34 

0.1049E+04 

0.1 826E+06 

0.1435E+05 

0,3 943 E+02 

35 

0.1 073E+04 

0.1796E+06 

0.1357E+05 

0.3953E+02 

36 

0.1097E+04 

0.1766E+06 

0.1279E+05 

0.3975E+02 

37 

0,1 121E+04 

0.1725E+06 

0.1 201 E+05 

0.3896E+02 

38 

0,1 145E+04 

0.1675E+06 

0.11 22E+05 

0.381 2E+02 

39 

0,1 169E+04 

0.1619E+06 

0.1044E+05 

0.3697E+02 

40 

0.1 240E+04 

0.1500E+06 

0.8336E+04 

0.3 407 E+02 
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v. 


TIME 

(SEC) 


ALTITUDE VELOCITY ANGLE ATTACK 
(FT) (FT/SEC) (DEG) 


41 

42 

43 

44 

45 

46 

47 

48 

49 

50 


0.1302E+04 
0.1364E+04 
0.1426E+04 
0.1488E+04 
0.1550E+04 
0.1612E+04 
0.1674E+04 
0.1736E+04 
0.1860E+04 
0.1 925E+04 


0.1339E+06 
0.1 174E+06 
0.1062E+06 
0.891 0E+05 
0.7635E+05 
0.5744E+05 
0.4167E+05 
0.2760E+05 
0.3337E+04 
-0.3000E+01 


0.6757E+04 
0.5342E+04 
0.4064E+04 
0.2917E+04 
0.1915E+04 
0.1151E+04 
0.8000E+03 
0.681 0E+03 
0.51 20E+03 
0.2020E+03 


0.2820E+02 

0.2305E+02 

0.2027E+02 

0.1647E+02 

0.1086E+02 

0.7790E+01 

0.7720E+01 

0.7090E+01 

0.3780E+01 

0.1220E+01 


DO 

NO 

DO 

NO 


YOU WISH TO CHANGE ANY OF THE TRAJECTORY INPUT ? 
YOU WISH TO WRITE THIS INPUT TO A FILE ? 


TRAJECTORY INPUT IS COMPLETE 


ATMOSPHERE DATA 


OPTIONS 


1. 1962 U.S. STANDARD ATMOSPHERE 

2. WIND TUNNEL OPTION 

3. INPUT ATMOSPHERIC DATA (ALT, T- INF, PH NF) 

4. 1963 PATRICK AIR FORCE BASE ATMOSPHERE 

5. 1971 VANDENBERG REFERENCE ATMOSPHERE 


OPTION SELECTED ? 
4 


1963 PATRICK AIR FORCE BASE ATMOSPHERE 

IS THIS OPTION CORRECT ? 

YES 


DO YOU WANT TO RUN A 
YES 


HEATING INDICATOR ? 


HEATING INDICATOR 


FAY AND RIDDELL 
RADIUS » 1 FT SPHERE 
WALL TEMP » 0 DEG. F 
LEWIS NO. « 1.0 
SUMMARY PRINT ONLY 


CONTROL FLAGS 


YOU ARE COMPLETING INPUT FOR CASE 1 
WHAT IS THE BODY POINT NUMBER FOR CASE 1 ? 

SHOULD LANMIN CREATE AN OUTPUT FILE FOR CASE 1 ? 

OUTPUT UNITS OPTIONS 

0. ENGLISH 

1. METRIC 

OPTION SELECTED ? 

0 . 


AS CASE 1 THE SAME STREAMJNE 

NO 



INPUT OONTROL FLAG 

' • Ksa Tin£ ' Ti "'* « 

2. END OF INPUT (LAST CASE) 

* 85 8 Si SSSf 

OPTION SELECTED ? 

2 . 


***** INPUT OOMPLETE FOR CASE 1 »#»#* 

00 YOU WANT TO MAKE ANY MODIFICATIONS TO CASE I 1 
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TRAJECTORY FILE LISTING 
HEATING INDICATOR CASE 


ORIGINAL PAGE la 
OF POOR QUALITY 


TRAJ. (ORB ITER) 


STS-1 REENTRY 
0 50.00000 
0.3000E+00 
0.4530E+02 
0.9030E+Q2 
0.1353E+03 
0 . 1 803E+03 
0.2253E+03 
0.2703E+Q3 
0.3153E+03 
0 .3603E+03 
0.4053E+03 
0.4503E+03 
0.4953E+03 
0.5453E+03 
0.5693E+03 
0.5933E+03 
0.6173E+03 
0.641 3E+03 
0.6653E+03 
0.6893E+03 
0.7133E+03 
0.7373E+03 
0.7613E+03 
0.7853E+03 
0.8093E+03 
0.8333E+03 
0.8573E+03 
0.881 3E+03 
0.9053E+03 
0.9293E+03 
0.9533E+03 
0.9773E+03 
0.1001E+04 
0.1025E+04 
0.1049E+04 
0.1073E+04 
0.1097E+04 
0.1 121E+04 
0.11 45E+04 
0.1 169E+04 
0.1240E+04 
0.1302E+04 
0.1364E+04 
0.1426E+04 
0.1488E+04 
0.1 550E+04 
0.1612E+04 
0.1674E+04 
0.1736E+04 
0.1860E+04 
0.1 925E+04 


0.3963E+06 
0.3738E+06 
0.351 5E+06 
0.3295E+06 
0.3080E+06 
0.2880E+06 
0.2697E+06 
0.2563E+06 
0.2502E+06 
0.2470E+06 
0.2446E+06 
0.2425E+06 
0.2400E+06 
0.2388E+06 
0.2377E+06 
0.2366E+06 
0.2353E+06 
0.2337E+06 
0.2320E+06 
0.2323E+06 
0.2303E+06 
0.2279E+06 
0.2252E+06 
0.2223E+06 
0.2190E+06 
0.2153E+06 
0.2114E+06 
0.2076E+06 
0.2056E+06 
0.2025E+06 
0.1 972E+06 
0.1 922E+06 
0.1872E+06 
0.1826E+06 
0.1796E+06 
0.1766E+06 
0.1725E+06 
0.1675E+06 
0.1619E+06 
0.1500E+06 
0.1339E+06 
0.1 174E+06 
0.1062E+06 
0.891 0E+05 
0.7635E+05 
0.5744E+05 
0.4167E+05 
0.2760E+05 
0.3337E+04 
-0.3000E+01 


0.2457E+Q5 

0.2459E+05 

0.2462E+05 

0.2464E+05 

0.2466E+05 

0.2466E+05 

0.246 1E+05 

0.2447E+05 

0.2422E+05 

0.2392E+05 

0.2361E+03 

0.2328E+05 

0.2288E+05 

0.2268E+05 

0.2247E+05 

0.2225E+05 

0.2202E+05 

0.2178E+05 

0.2153E+05 

0.2126E+05 

0.2098E+05 

0.2068E+05 

0.2036E+05 

0.2002E+05 

0.1 965E+05 

0,1 924E+05 

0.1879E+05 

0.1830E+05 

0.1776E+05 

0.1718E+05 

0.1654E+05 

0.1587E+05 

0.1513E+05 

0.1435E+05 

0.1357E+05 

0.1279E+05 

0.1201E+05 

0.1 122E+05 

0.1044E+05 

0.8336E+04 

0.6757E+04 

0.5342E+04 

0.4064E+04 

0.291 7E+04 

0.191 5E+04 

0.1 151E+04 

0.8000E+03 

0,681 0E+03 

0.5120E+03 

0.2020E+03 
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0.41 13E+02 
0.41 26E+Q2 
0.41 21 E+02 
0.4046E+02 
Q.4065E+02 
0.41 90E+02 
0.3939E+02 
0.4071 E+02 
0.4174E+02 
0.3995E+02 
0.3928E+02 
0.3966E+02 
0.3907E+02 
0.3932E+02 
0.3902E+02 
0.3925E+02 
0.3948E+02 
0.3997E+02 
0.4001 E+02 
0.4025E+02 
0.4046E+02 
0.401 4E+02 
0.401 6E+02 
0.401 6E+02 
0.4027E+02 
0.4040E+02 
0.401 2E+02 
0.4030E+02 
0.4001 E+02 
0.4200E+02 
0.4090E+02 
0.4075E+02 
0.3 987 E+02 
0.3943E+02 
0.3953E+02 
0.3975E+02 
0.3896E+02 
0.381 2E+02 
0.3697E+02 
0.3407E+02 
0.2820E+02 
0.2305E+02 
0.2027E+02 
0.1647E+02 
0.1086E+02 
0.7790E+01 
0.7720E+01 
0.7090E+01 
0.3780E+01 
-0.1220E+01 


PREM IN OUTPUT FILE . 

HEATING INDICATOR CASE t'/v 'i . 

OF POOR QUAI T»> 


REENTRY TRAJ. (ORB ITER) 

PAFB REF. HTG. 

IND. CASE 



0.000 

25.000 

500.000 



14.000 

1200.000 

29.000 



1925.500 

1.000 




50.000 





-300 

396300.000 

24570.000 

41.130 

0.000 

4; .300 

373800.000 

24590.000 

41.260 

0.000 

90.30C 

351500.000 

24620.000 

41.210 

0.000 

135.300 

329500.000 

24640.000 

40.460 

0.000 

180.300 

308000.000 

24660.000 

40.650 

0.000 

225.300 

288000.000 

24660.000 

41 .900 

0.000 

270.300 

269700.000 

24610.000 

39.390 

0.000 

315.300 

256300.000 

24470.000 

40.710 

0.000 

360.300 

250200.000 

2 4220.000 

41 .740 

0.000 

405.300 

247000.000 

23920.000 

39.950 

0.000 

450.300 

244600.000 

23610.000 

39.280 

0.000 

495.300 

242500.000 

23280.000 

39.660 

0.000 

545.300 

240000.000 

22880.000 

39.070 

0.000 

569.300 

238800.000 

22680.000 

39.320 

0.000 

593.300 

237700.000 

22470.000 

39.020 

0.000 

617.300 

236600.000 

22250.000 

39.250 

0.000 

641 .300 

235300.000 

22020.000 

39.480 

0.000 

665.300 

233700.000 

21780.000 

39.970 

0.000 

689.300 

232000.000 

21530.000 

40.010 

0.000 

713.300 

232300.000 

21260.000 

40.250 

0.000 

737.300 

230300.000 

2 0980.000 

40.460 

0.000 

761.300 

227900.000 

20680.000 

40.140 

0.000 

785.300 

225200.000 

20360.000 

40.160 

0 000 

809.300 

222300.000 

20020.000 

40.160 

0.000 

833.300 

219000.000 

19650.000 

40.270 

0.000 

857.300 

215300.000 

19240.000 

40.400 

0.000 

881.300 

211400.000 

18790.000 

4C.170 

0.000 

905.300 

207600.000 

18300.000 

40.300 

0.000 

929.300 

205600.000 

17760.000 

40.010 

0.000 

953.300 

202500.000 

17180.000 

42.000 

0.000 

977.300 

197200.000 

16540.000 

40.900 

0.000 

1001.000 

192200.000 

15870.000 

40.750 

0.000 

1025.000 

187200.000 

15130.000 

39.870 

C.000 

1049.000 

182600.000 

14350.000 

39.430 

0.000 

1073.000 

179600.000 

13570.000 

39.530 

0.000 

1097.000 

176600.000 

12790.000 

39.750 

0.000 

1121.000 

172500.000 

12010.000 

38.960 

0.000 

1145.000 

167500.000 

11220.000 

38.120 

0.000 

1169.000 

161900.000 

10440.000 

36.970 

0.000 

1240.000 

150000.000 

8336.000 

34.070 

0.000 

1302.000 

133900.000 

6757 000 

28.200 

0.000 

1364.000 

117400.000 

5342.000 

23.050 

0.000 

1426.000 

106200.000 

4064.000 

20.270 

0.000 

1488.000 

89100.000 

2917.000 

16.470 

0.000 


B~8 


ORIGINAL PAG, F 
OF POOR QUALITY 


1550.000 

1612.000 

1674.000 

1736.000 

1860.000 

1925.000 

4.0000 11 1.0000 12 

90.0000 38 90.0000 315 

1001.0 

2.0000 643 1.0000 647 


76350.000 

1915.000 

10.860 

0.000 

57440.000 

1151.000 

7.790 

0.000 

41670.000 

800,000 

7.720 

0.000 

27600.000 

681 .000 

7.090 

0.000 

3337.000 

512.000 

3.780 

0.000 

- 3.000 

202.000 

- 1 .220 

0.000 

1.0000 31 
1 .0000 

38.0000 32 18.0000 


1.0000 641 2.0000 
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LANMIN OUTPUT 
HP AT I NG INDICATOR CASE I 

STS ! REENTRY TRW. (ORIITERI PAF1 REF. HTO. IW. CASE 


TIMING 


Tl 

- .000 SEC 

DTI 

« 25.000 SEC 

T2 

« 500.000 SEC 

DT2 

- 14.000 SEC 

T3 

> 1200.000 SEC 

DT3 

a 29.000 SEC 

T4 

« 1925.300 SEC 

DT CALC 

• 1.000 SEC 

DTEHP MAX 

' .000 SEC 


CONTROL PARAMETERS 

H(641)«3.000 M(642)«2.000 
WI643M.OOO H(644)* .000 
W(649)» .000 H(A46)« .000 
W(647)«1.000 H(646)e .000 
H<649)= .000 H(650)« .000 


HEAT TRANSFER 

HT METHOD » 1.000 

RN = 1.000 FT. 

I « .000 FT. 

N SUB L « 1.000 

N SUB T » 1.000 

PHI » .000 DEB. 

VIRT-L CJPT= .000 


MULTIPLICATION FACTORS 
K SIS L 1 « 2.000 

R9«Tla 1.000 


MATERIAL 

EMISSIVITY® .600 


FUMFIELD 


FLAG 

1 36. 

2 IS. 

ANGLE 

90.00 

90.00 



TRADITION 

OPTION 

TRAJECTORY 

TIME 

SEC. 

ALTITUDE 

FT 

VELOCITY 

FT/SEC 

ANGLE OF ATTACK 
DEG 

.30 

45.30 

90.30 

39*300. 

'373800. 

351500. 

24570. 

24590. 

24620. 

41.130 

41.260 

41.210 



B-1Q 



original ,, i. y 

OF POOR QUAUT\ 


135.30 

3295X. 

100.30 

306000. 

235.30 

288000. 

270.30 

2697X. 

315.30 

2563X. 

360.30 

250200. 

405.30 

2470X. 

450.30 

2446X. 

495.30 

242SX. 

545.30 

2400X. 

569.30 

238X0. 

593.30 

2377X. 

617.X 

2366X. 

641.X 

235300. 

665.X 

2337X. 

6S9.X 

2320X. 

713.X 

2323X. 

737.X 

2309X. 

761.X 

2279X. 

7BS.X 

22S2X. 

X9.X 

2223X. 

833.X 

2190X. 

857.X 

2153X. 

861.X 

2114X. 

905.X 

X76X. 

929.X 

20S6X. 

953.X 

202SX. 

977.X 

1972X. 

1X1 .00 

1922X. 

1025.X 

1872X. 

1049.X 

1826X. 

1073.X 

1796X. 

1097.X 

1766X. 

1121.X 

1725X. 

1145.X 

1675X. 

1169.X 

1619X. 

1240.X 

1500X. 

1302.X 

1339X. 

1364.X 

1174X. 

1426.X 

1062X. 

148B.X 

891X. 

1550. X 

76350. 

1612.X 

57440. 

1674.X 

41670. 

1736.X 

276X. 

1860.X 

3337. 

1925.X 

-3. 


1963 PAFB STD AfTOSPHERE 



C ‘ iv -W:u 

° F POu,* 

24640. 

40.460 

24660. 

40.6M 

24660. 

41 .9X 

24610. 

39.390 

24470. 

40.710 

24220. 

41.740 

23920. 

39.950 

23610. 

39.2X 

232X. 

39.660 

22880. 

39.070 

22680. 

39.320 

22470. 

39.0X 

222X. 

39.250 

22020. 

39.4X 

217X. 

39.970 

21530. 

40.010 

21260. 

40.250 

20980. 

40.460 

20680. 

40.140 

X360. 

40.160 

20020. 

40.160 

19650. 

40.270 

19240. 

40.4X 

18790. 

40. IX 

183X. 

40.3X 

17760. 

40.010 

17180. 

42.0X 

16540. 

40.9X 

15870. 

40.7X 

151X. 

39.870 

14350. 

39.4X 

13570. 

39.5X 

12790. 

39.7X 

1X10. 

X.960 

11220. 

X.1X 

10440. 

36.970 

8336. 

34.070 

6757. 

a. 2 x 

5342. 

23.0M 

4064. 

X.270 

2917. 

16.470 

1915. 

10.860 

US1. 

7.790 

BX. 

7.7X 

681. 

7.090 

512. 

3.7X 

X2. 

-1.2X 
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fid 


STS-I HEEKUt raw. (ORBITED POT REF. m,. ° F P °^ 


Tift 

SC 


ALT VEL 

aft ft/sc 


HACK ANGLE 
» ATTACH 


REYNOLDS 

NO./FT 


l 


.0 396.3 24970.0 19.29 

29.0 364.0 24961.0 20.99 

90.0 371.9 24993.1 21.69 

79.0 399.1 24609.6 23.30 

100.0 346.8 24624.3 24.21 

129.0 334.9 24639.4 29.23 

190.0 322.9 24646.9 26.14 

179.0 310.9 24697.6 26.88 

200.0 299.2 24660.0 27.62 

229.0 268.1 24660.0 27.90 

290.0 278.0 24632.6 27.87 

279.0 268.3 24999.4 27.99 

300.0 260.9 24917.6 27.01 

329.0 299.0 24416.1 26.96 

390.0 291.6 24277.2 26.22 

379.0 249.2 24122.0 29.92 

400.0 247.4 23995.3 29.65 

425.0 245.9 23784.3 25.39 

490.0 244.6 23612.1 25.14 

475.0 243.4 23428.9 24.88 

900.0 242.3 23242.4 24.63 

514.0 241.6 23130,4 24.47 

528.0 240.9 23018.4 24.32 

942.0 240.2 22906.4 24.16 

996.0 239.5 22790.8 24.01 

970.0 238.6 22673.9 23.85 

984.0 238.1 22551.4 23.69 

998.0 237.5 22426.9 23.93 

612.0 236.8 22298.6 23.36 

626.0 236.1 22166.6 23.19 

640.0 235.4 22032.9 23.01 

694.0 234.5 21893.0 22.83 

668.0 233.5 21751.9 22.63 

682.0 232,5 21606.0 22.44 

696.0 232.1 21494.6 22.26 

710.0 232.3 21297.1 22.10 

734.0 231.4 21139.2 21.90 

738.0 230.2 20971,3 21.66 

792.0 228.8 20796.3 21.43 

766.0 227.4 20617.3 21.19 

780.0 225.8 20430.7 20.93 

794.0 224.1 20236.8 20.66 

808.0 222.5 20038.4 20.38 

822.0 220.6 19824.2 20.09 

836.0 218.6 19603.9 19.78 

890.0 216.4 19364.7 19.46 

864.0 214.2 19114.4 19.12 

878.0 211.9 18851.9 18.77 

892.0 209.7 18571.5 18.41 

906.0 207.5 18284.3 18.09 


JtAT COEF REC ENTHALPY RAD EOUIL HEAT RATE 
UW8FT-5 BTU/JH deb f 1TU/9F7-S 


SAT LOAD 
ITU/BFT 


PRESSURE 

LB/SFT 


FLOW 

TYPE 


41.13 
41.20 
41.29 
41.23 
41.05 
40.63 

40.92 
40.63 
41.20 

41.89 
40.52 

39.93 

40.26 

40.93 

41.90 

41.16 

40.16 
39.66 
39.28 
39.49 
39.60 
39.44 

39.27 
39.11 
39.18 
39.31 

39.14 

39.07 
39.20 

39.33 
39.47 
39.74 
39.97 
40.00 

40.08 
40.22 

40.34 
40.45 
40.26 

40.14 
40.16 
40.16 
40.16 
40.22 

40.28 
40.36 
40.32 
40.16 
40.20 

40.29 


.239*001 
.426*001 
.793*001 
.157*002 
.296*002 
.986*002 
.119*003 
.240*003 
.486*003 
.936*003 
.165*004 
.273*004 
.382*004 
.495*004 
.970*004 
.630*004 
.675*004 
.712*004 
.748*004 
.780*004 
.813*004 
.833*004 
.853*004 
.874*004 
.895*004 
.916*004 
.939*004 
.995*004 
.974+004 
.997*004 
.102*005 
.109*005 
.109*005 
.112*005 
.113*005 
.112*005 
.115*005 
.119*009 
.129*009 
.131*005 
.138*005 
.146*005 
.154*005 
.164*005 
.174*005 
.187*005 
.200*009 
.214*005 
.229*009 
.244*005 


.146-003 
.189-003 
.246-003 
.331-003 
.436-003 
.597-003 
.827-003 
.119-002 
.199-002 
.220-002 
.293-002 
.383-002 
.458-002 
.523-002 
.562-002 
.991-002 
.613-002 
.638-002 
.661-002 
.681-002 
.702-002 
.712-002 
.720-002 
.728-002 
.736-002 
.742-002 
.746-002 
.749-002 
.752-002 
.756-002 
.761-002 
.768-002 
.779-002 
.782-002 
.779-002 
.767-002 
.771-002 
.789-002 
.802-002 
.820-002 
.840-002 
.862-002 
.884-002 
.909-002 
.939-002 
.978-002 
.102-001 
.106-001 
.109-001 
.112-001 


.122*009 

.122*009 

.122*009 

.122*005 

.122*009 

.122*005 

.122*009 

.122*009 

.122*009 

.122*009 

.122*005 

.122*009 

.121*005 

.120*005 

.118*005 

.117*005 

.119*009 

.114*005 

.112*005 

.110*005 

.109*009 

.108*009 

.107*005 

.106*005 

.105*005 

.104*005 

.102*005 

.101*005 

.100*005 

.990*004 

.978*004 

.966*004 

.954*004 

.941*004 

.928*004 

.919*004 

.901*004 

.887*004 

•873*004 

.858*004 

.843*004 

.827*004 

.811*004 

.794*004 

.777*004 

.798*004 

.739*004 

.719*004 

.699*004 

.677*004 


1001.0 

1097.1 

1204.6 
1327,9 

1453.8 

1609.3 

1782.4 

1972.1 

2179.5 

2391.8 

2998.5 

2803.9 
2944.3 

3090.0 

3100.7 

3132.1 

3151.8 

3173.0 

3190.7 

3202.8 

3213.8 

3217.1 

3217.5 

3217.8 

3217.5 

3214.8 

3208.8 

3202.2 

3194.9 

3188.2 
3181.7 

3177.1 

3172.3 

3167.3 

3191.1 

3123.1 

3113.4 

3113.6 

3116.2 

3119.0 

3122.4 

3126.7 

3128.6 

3132.6 

3138.3 

3190.2 

3199.0 

3164.6 

3164.3 

3150.7 


.424*000 

.688*000 

.116*001 

.209*001 

.371*001 

.689*001 

.100*002 

.139*002 

.193*002 

.266*002 

.353*002 

.461*002 

.548*002 

.621*002 

.699*002 

.684*002 

.701*002 

.719*002 

.734*002 

.749*002 

.755*002 

.759*002 

.760*002 

.761*002 

.761*002 

.760*002 

.799*002 

.751*002 

.745*002 

.740*002 

.736*002 

.733*002 

.730*002 

.727*002 

.715*002 

•693*002 

.666*002 

.688*002 

.69K002 

.699*002 

.699*002 

.703*002 

.707*002 

.712*002 

.719*002 

.731*002 

.741*002 

.749*002 

.752*002 

.743*002 
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.000 

.172*002 

.463*002 

.984*002 

.191*003 

.363*003 

.614*003 

.961*003 

.144*004 

.211*004 

.299*004 

.414*004 

.551*004 

.707*004 

.872*004 

.104*005 

.122*005 

.140*005 

.198*009 

.177*009 

.196*009 

.206*009 

.217*005 

.228*005 

.238*005 

.249*009 

.299*005 

.270*005 

.280*005 

.291*005 

.301*005 

.311*505 

.321*005 

.332*009 

.342*005 

.351*009 

.361*005 

.371*005 

.380*005 

.390*005 

.400*005 

.410*005 

.420*005 

.429*005 

.440*005 

.430*005 

.460*005 

.471*005 

.481*005 

.492*005 


.266-001 Ufl 
.445-001 LAN 
.758-001 LAN 
.133*000 LAN 
.229*000 LAN 
.429*000 LAN 
.818*000 LAN 
.157*001 LAN 
.299*001 LAN 
.965*001 LAN 
.997*001 LAN 
.149*002 LAN 
.241*002 LAN 
.314*002 LAN 
.362*002 LAN 
.400*002 LAN 
.431*002 LAN 
.466*002 LAN 
.900*002 LAN 
.931*002 LAN 
.564*002 LAN 
•580*002 LAN 
.593*002 LAN 
.606*002 LAN 
.619*002 IAN 
•630*002 LA" 
.636*002 LAN 
.643*002 LAN 
.648*002 LAN 
.Af<o+O02 LAN 
.664*002 LAN 
.676*002 LAN 
.688*002 LAN 
.701*002 LAN 
.697*003 LAN 
.676*002 UW 
.684*002 Ltf 

.vwj+ 002 urn 
•740*002 LAN 
.773*002 LAN 
.812*002 Utft 
.854*002 LAN 
.898*002 UW 
.991*002 LAN 
.101*003 um 
.110*003 LAN 
.119*003 LAN 
.129*003 LAN 
.138*003 LAN 
.143*003 LAN 



OR!C;ir,T 

OF POO.'V 



? 


8T8-I REENTRY TRAJ. CORBITER) PAFB REF. HTQ. IND. CASE I.P, NO. 100! 

(CONTINUED) 


2? in rr/*r *5? £££ * AT m Ki EHTH * J,Y ^ ^RUIL HEAT RATE (CAT LOAD PREB8URE 

SEC AFT FT/8EC NO ATTACK NO. /FT LIR/SFT-S ITU/LW BED F ITU/8F7-6 ITU/8FT LI/8FT 


W.0 206.4 17969.3 17.70 40.12 .250*009 .111-00) 

934.0 203.0 17646.4 17.33 40.40 .296*009 .111-00) 

946.0 203.2 17306.1 16.94 41.96 .270*009 .113-00) 

962.0 200.6 16946.0 16.51 41.60 .269*009 .116*00) 

976.0 197.5 16974.7 16.09 40.96 .319*003 ,120*00) 

990.0 194.9 16181.0 15.99 40.82 .311*009 .124-001 

1004.0 191.6 19777.9 15.12 40.64 .368*009 .127-001 

1018.0 186.7 19349.8 14.63 40.13 .396*003 .130-001 

1032.0 165.9 14902.5 14.15 39.74 .423*009 .132-001 

1046.0 183.2 14447.5 13.66 39.49 .490*009 ,134-001 

1060.0 181.2 13992.5 13.19 39.48 .466*009 .133*001 

1074.0 179.5 13937.5 12.73 39.54 .478*005 .134*001 

1088.0 177.7 13082.5 12.27 39.67 .491*009 .132-001 

1102.0 175.7 12627.5 11.81 39.59 .907*005 .132-001 

1116.0 173.4 12172.5 11.36 39.12 .530*009 .133-001 

1130.0 170.6 11713.8 10.89 38.64 .560*009 .134-001 

1144.0 167.7 11292.9 10.43 38.16 .995*009 .139-001 

1198.0 164.5 10797.5 9.98 37.90 .639*005 .138-001 

1172.0 161.4 10391.1 9.95 36.85 .682*005 .139*001 

1186.0 159.1 9936.2 9.17 36.28 .716*005 .140-001 

1200.0 156.7 9521.4 8.79 39.70 .793*005 .140-001 

1229.0 151.8 8662.0 8.03 34.52 .834*005 .140-001 

1298.0 145.3 7877.6 7.35 32.37 .997*005 .144-001 

1287.0 137.8 7139.0 6.73 29.62 .126*006 .152-001 

1316.0 130.2 6437.5 6.14 27.04 .160*006 .199-001 

1349.0 122.5 9775.6 5.57 24.63 .205*006 .168-001 

1374.0 115.6 5139.9 9,01 22.60 .293*006 .172-001 

1403.0 110.4 4936.1 4.46 21. X .288*006 .168-001 

1432.0 104.9 3953.0 3.92 19.90 .333*006 .164-001 

1461.0 96.5 3416.5 3.43 18.12 .428*006 .167-001 

1490.0 88.7 2884.7 2.92 16.29 .534*006 .167-001 

1519.0 82.7 2416.0 2.46 13.66 .603*006 .197-001 

1948.0 76.8 1947.3 2.00 11.04 .663*006 .142-001 

1977.0 68.1 1982.3 1.69 9.92 .861*006 .140-001 

1606.0 99.3 1224.9 1.30 8.09 .110*007 .134-001 

1635.0 51.6 1020.8 1.09 7.76 .138*007 .144-001 

1664.0 44.2 896.6 .90 7.73 .157*007 .152-001 

1693.0 37.4 763.5 .78 7.53 .179*007 .164-001 

1722.0 30.8 707.9 .70 7.23 .174*007 .180-001 

1791.0 24.7 660.6 .63 6.69 .212*007 .194-001 

1780.0 19.0 621.0 .98 5.92 .230*007 .208-001 

1809.0 13.3 981.5 .93 5.14 .290*007 .222-001 

1838.0 7.6 542.0 .49 4.37 .271*007 .239-001 

1867.0 3.0 478.6 .43 3.24 .269*007 .237-001 

1896.0 1.5 340.3 .X 1.01 .198*007 .202-001 

1929.0 .0 202.0 .18 -1.22 .122*007 ,198-001 

1994.0 .0 X2.0 .18 -1.22 .122*007 .198-001 


.699*004 

3110.2 

.713*002 

.502*009 

.141*003 

.632*004 

X77.2 

.689*002 

.511*005 

.142*003 

.608*004 

3051.5 

.673*002 

.521*005 

.147*003 

.964*004 

XX.8 

.663*002 

.530*009 

.155*003 

.599*004 

3021.0 

.698*002 

.539*009 

.168*003 

.933*004 

3000.0 

.647*002 

.548*005 

.179*003 

.908*004 

2974.1 

.632*002 

.597*009 

.189*003 

.481*004 

29X.9 

.609*002 

.566*009 

.197*003 

.494*004 

2899.6 

.987*002 

.974*005 

.206*003 

.428*004 

2849.0 

.560*002 

.982*009 

.213*003 

.401*004 

2781.5 

.520*002 

.989*005 

.212*003 

.377*004 

2723.2 

.489*002 

.596*005 

.214*003 

.390*004 

2644.3 

.447*002 

.602*009 

.212*003 

.330*004 

2990.2 

.421*002 

.608*009 

.214*003 

.X7*004 

2927.6 

.393*002 

.613*009 

.216*003 

.285*004 

2465.6 

.367*002 

.619*009 

.221*003 

.264*004 

2403.6 

.343*002 

.623*009 

.227*003 

.244*004 

2344.2 

.321*002 

.628*009 

.235*003 

.226*004 

22X.5 

.299*002 

.632*009 

.241*003 

.209*004 

2212.7 

.276*002 

.636*009 

.142*003 

.193*004 

2142.8 

.254*002 

.639*009 

.243*003 

.161*004 

1990.5 

.210*002 

.646*005 

.245*003 

.139*004 

1X5.0 

.179*002 

.691*005 

.263*003 

.113*004 

1726.8 

.155 '002 

.695*005 

.299*003 

.937*003 

1985.0 

.132*002 

.699*009 

.334*003 

.773*003 

1437.0 

.111*002 

.662*009 

.380*003 

.631*003 

1268.1 

.894*001 

.669*005 

.410*003 

.514*003 

1085.0 

.677*001 

.667*009 

.405*003 

.413*003 

895.3 

.497*001 

.668*009 

.405*003 

.332*003 

715.5 

.371*001 

.669*005 

.442*003 

.263*003 

929.9 

.259*001 

.670*009 

.467*003 

.212*003 

368.9 

.199*001 

.671*009 

.440*003 

.170*003 

2X.9 

.842*000 

.671*009 

.395*003 

.141*003 

116.0 

•433*000 

.671*009 

.424*003 

.118*003 

26.7 

.107*000 

.671*009 

.448*003 

.108*003 

-13.0 

-.297-001 

.671*005 

.511*003 

.109*003 

-25.4 

-.781-001 

.671*009 

.583*003 

.108*003 

-15.5 

-.458-001 

.671*005 

.716*003 

.112*003 

4.1 

.357-001 

.671*009 

.894*003 

.117*003 

23.7 

.131*000 

.671*009 

.110*004 

.121*003 

40.9 

.228*000 

.671*009 

.132*004 

.125*003 

59.5 

.323*000 

.671*009 

.199*004 

.128*003 

67.8 

.413*000 

.671*005 

.190*004 

.130*003 

79.5 

.463*000 

.671*009 

.217*004 

<129*003 

70.1 

.371*000 

.671*009 

.215*004 

.128*003 

67.5 

.286*000 

.672*005 

.217*004 

.128*003 

67.9 

•286*000 

.672*005 

.217*004 
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Msfl l&M£i Staple Cue 

The example ease given herein corresponds to the results presented in 
Volume I Fig. 7.15 and 7.16 for X/L - 0.1, 0.2 and 0.3. The problem solved is 
heating to the bottom centerline of an orbiter naing the effective running 
length concept. In the wind tnnnel case, time is used to denote different run 
conditions. Since the effective running length changes with angle of attack and 
angle of attack with run conditions, the time dependent geometry option is used. 

The PREMIN interactive listing is followed by the PREMXN output file. The 
output listing from LANMIN resulting from the PREMIN output file is given last. 
The output lisiting contains the intermediate print as well as the summary page. 
The description of the intermediate output symbols and units is given in Table 
4.1. The summary print follows each case. Cases 1, 2, and 3 are for body 
points 100, 200 and 300 corresponding to X/L « 0.10, 0.20, and 0.30 respective- 
ly. The heating load in the summary print has no meaning for a wind tunnel case 
where time is a run number indicator. 




INTERACTIVE PREMIN 
WIND TUNNEL RUN CASE 



INTERACTIVE INPUT TO LARC MINIVER - LANMIN 


OUTPUT FILE METHOD 

1. CREATE A NEW OUTPUT FILE 

2. MODIFY AN EXISTiNG OUTPUT FILE 

OPTION SELECTED ? 

1 

engush WANT T ° ,NPUT data ,n ENGL,SH metric 7 
INTERACTIVE INPUT FOR CASE 1 


SPECIFY PRINTOUT INTERVALS 

INITIAL TIME (SEC) 

1.0 

PRINTOUT INTERVAL 1 (SEC) DELTA TIME 

1.0 

SECOND TIME (SEC) 

4.0 

PRINTOUT INTERVAL 2 (SEC) DELTA Tll€ 

0.0 

THIRD TIME (SEC) 

0.0 

PRINTOUT INTERVAL 3 (SEC) DELTA TIME 

0.0 

FOURTH TIME (SEC) 

0.0 


ARE THE PRINTOUT TIMES CORRECT ? 
YES 


DO YOU HAVE A TRAJECTORY INPUT FILE ? 

WHAT IS THE NUI©ER OF TRAJECTORY POINTS ? (50 TRAJ.PTS. MAXIMUM) 

WILL BETA VALUES BE INPUT ? 

NO 

TYPE IN THE FOLLOWING TRAJECTORY VARIABLES SEPERATED BY COMMAS 
‘ " " “ “ " * 50 TIMES MAXIMUM ........... 

TIME(SEC),ALTITUDE(FT),VELOCITY(FT/SEC),ANGLE OF ATTACK(DEG) 

1.0. 0.0.3816.0.20.0 
2 

2.0. 0.0.3821.0.29.86 

3 

3.0. 0.0.3821.0.39.98 

4 

4.0. 0.0.3859.0.39.98 


TIME 

(SEC) 


ALTITUDE VELOCITY ANGLE ATTACK 
(FT) (FT/SEC) (DEG) 


1 0.1000E+01 

2 0.2000E+01 

3 0.3000E+01 

4 0.4000E+01 


O.OOOOE+OO 0.381 6E+04 0.2000E+02 
O.OOOOE+OO 0.3821 E+04 0.2986E+02 
O.OOOOE+OO 0.3821 E+04 0.3998E+02 
O.OOOOE+OO 0.3859E+04 0.3998E+02 


DO YOU WISH TO CHANGE ANY OF THE TRAJECTORY INPUT ? 
NO 

DO YOU WISH TO WRITE THIS INPUT TO A FILE ? 

YES 

WHAT IS THE NEW FILE NAME ? 

TRAJ.FIL 


TRAJECTORY INPUT IS COMPLETE 


ATMOSPHERE DATA 


OPTIONS 


1 . 

2 . 

3. 

4. 

5. 


1962 U.S. STANDARD ATMOSPHERE 
WIND TUNNEL OPTION 


INPUT ATMOSPHERIC DATA(ALT, T- I NF, P- I NF ) 
1963 PATRICK AIR FORCE BASE ATMOSPHERE 
1971 VANDENBERG REFERENCE ATMOSPHERE 


OPTION SELECTED ? 

2.0 


WIND TUNNEL OPTION 

IS THIS OPTION CORRECT ? 
YES 


WIND TUNNEL OPTION 

EKS 


------ 4 values required 

T-INF(R),P-lNF(lfl/SFT) 

94 . 7 , 8.064 

2 

95 . 0 , 8.208 

3 

95 . 0 , 8.208 

4 

96 . 9 , 12.672 

ARE ALL INPUTS CORRECT ? 



DO YOU WANT TO RUN A HEATING INDICATOR ? 


HEAT TRANSFER l€THOD 


OPTIONS 


1, 

2 . 

3 . 

4 . 

5 . 

6 . 

7 . 

8 , 
9 . 

10 . 

11 , 


Uft!i? PHERE stagnation point 

CATO/ JOHNSON SWEPT CYLINDER 
ECKERT REF § ENTHALPY FLAT PLATF lUFTwnn 
ECKERT/SPAULD INS-CHI FLAT PLATE ttTHOD 
BOEING RHO-MU FLAT PLATE METHOD 

B0^NG T ^O^M| l i A cwCBr S lJ! Pr C ' fLINDER WTHOO 

rcic/^S!?"^ cvlinder nethod 

UBmSMnft «A?fJH SPHERE DISTOIBLrri0N 

FLAP REATTACH^NT heating 
fin-plate peak interference heating 


OPTION SELECTED ? 

4.0 

SHOULD RAREFIED aOW HEATING BE INCLUDED ? 

YES™ E HEAT transfer option correct ? 


4. ECKERT / SPAULD I NG-CH I FLAT PLATE METHOD 

RUNNING LENGTH (FT) ? 

.2097 

SURFACE DISTANCE TO START OF TURBULENT B.L. 

THIS DISTANCE IS SUBTRACTED FROM THE RUNNING LENGTH 
FOR TURBULENT HEATING CALCULATIONS. 

DESIRED LENGTH ? 

0.0 

IS AN AUTOMATIC VIRTUAL ORIGIN CORRECTION DESIRED ? 
NO 

TURBULENT MANGLER FACTOR ? 

.6667 

LAMINAR MANGLER FACTOR ? 

1.0 

REYNOLDS-ANALOGY FACTOR 

0. COLBURN 

1. VON KARMAN 

DESIRED FACTOR ? 

1 . 

ANY CHANGES ? 

NO 


* * * WALL CONDITIONS * * * 

WALL TEMPERATURE (DEG F) ? 

0.0 

WALL EMISSIVITY ? 

0.8 


* * * CONTINUATION OPTION * * * 

DO YOU WANT TO PROVIDE AN INITIAL HEATING LOAD GT 0.0 (BTU/ SQ.FT) 7 


DO YOU WANT TO USE A HEAT TRANSFER MULTIPLICATION !€THOD 
NO 


TRANSITION OPTIONS 


OPTIONS 


1. TIME DEPENDENT: LAM TO TURB 

2. TIME DEPENDENT: TLB TO LAM 

3. REYNOLDS NO. DEPENDENT 

4. RE-THETA 

5. MDAC-E TRANSITION 

6. MDAC-E TABLE LOOK-UP 

7. NAR RE VS re TABLE LOOK-UP 

8. RE-THETA/ME 


OPTION SELECTED ? 

1.0 


1. TIME DEPENDENT: LAM TO TURB 

TRANSITION BEGINS AT TIME(SEC) ? 
3.2 

FULLY TURBULENT AT TIME(SEC) ? 
3.6 

ANY CHANGES ? 

NO 


DO YOU WANT CROSS FLOW ADJUSTMENT OPTION ? 
NO 


•O 


FLOWFIELD AND LOCAL PRESSURE OPTIONS 


FLOWFIELD - - ------ PRESSURE 

-1. FLOWFIELD TYPE NOT NEEDED -1. PRESSURE TYPE NOT NEEDED 

1. SHARP WEDGE SHOCK ANGLE 1. INPUT CP VS MACH NO. TABLE 

2. SHARP CONE SHOCK ANGLE 2. TANGENT WEDGE PRESSURE 

3. OBLIQUE AND NORMAL SHOCK (90 DEG) 3. TANGENT CONE PRESSURE 

4. PARALLEL SHOCK (PRES NOT NEEDED) 4. OBLIQUE SURFACE PRESSURE 

5. MODIFIED NEWTONIAN 

6. PRANDTL-MEYER EXP. (FF NOT NEEDED) 

* * * * INPUT OPTIONS IN PAIRS WITH ASSOCIATED DELTA ANGLES * » « * 

(NOTE: TO SIGNIFY END OF CASE USE -1.0 FOR BOTH FF AND PRESS. OPTIONS.) 
FLOWFIELD, DEL ANGLE, PRESSURE, DEL ANGLE 
1 

2 . 0 , 2 . 8 , 2 . 0 , 2.8 

2 

-1 .0,0. 0,-1 .0,0.0 


TOTAL EFFECTIVE ANGLE 

SET1 

FF SHRP-CONE ALPHA + 2.800 

P TAN-WEDGE ALPHA + 2.800 


ANY CHANGES ? 
NO 
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DO YOU WANT TO USE TIME DEPENDENT GEOMETRY ? 


TIME DEPENDENT GEOMETRY 


NIJffcER OF TIME DEPENDENT ENTRIES ? (MAX-10) 

KM!o^M7 D ;i.'r , ' LEN6TH<FT, ' SLOPE * SWEP(DE6) 

2. 0. 0.0.. 1239.0.0 

3. 0. 0.0.. 0858.0.0 

4. 0. 0.0. .0858.0.0 
ANY CHANGES ? 

NO 


CONTROL FLAGS 


YOU ARE COMPLETING INPUT FOR CASE 1 
toS T ,S ™ E B0DY P°fNT NUI*ER FOR CASE 1 ? 

SHOULD LANMIN CREATE AN OUTPUT FILE FOR CASE 1 ? 


0 . 

1 . 

2 . 


PRINT OONTROL OPTIONS 
DETAILED PRINTOUT 
DETAILED PLUS SUMMARY PRINTOUT 
SUMMARY PRINTOUT 


NOTE: IF 
OPTION 1. 


AN OUTPUT FILE IS TO BE CREATE 
OR 2. MUST BE SELECTED 


EITHER 


OPTION SELECTED ? 

1.0 


OUTPUT UNITS OPTIONS 

0. ENGLISH 

1. METRIC 

OPTION SELECTED ? 

0 . 


le ? AL0NG THE SAME STREAM I NE 

AS CASE 1 ? 

NO 


ORIGINAL PAGE 19 
0F POOR QUALITY 


INPUT CONTROL FLAG 

XmmZu'. s us cas E title ' t,h,ng » ™*. 

*• END OF INPUT (LAST CASE' 

i HIS SSS3SS- 

(INITIALLY ZERO W ARRAY) 1 itle *T|MING,TRAJ. AND CASE DATA. 

6. SAME AS (3) EXCEPT f^TIAL IZ^ZERO^U F ??r PREV,0US CASE 

wtr| ,n, TIALIZE ZERO ALL TIMING AND CASE DATA 

OPTION SELECTED ? 

1 .0 


***** INPUT OOMPLETE FOR CASE 1 ***** 

DO YOU WANT TO MAKE ANY MODIFICATIONS TO CASE 1 ? 


CREATE OUTPUT FILE 

PREMIn!oUT E NAME ° F ™ E ° UTPUT FILE T0 BE CREATED ? 

HHAT IS THE TITLE FOR CASE 1 ? (NOTE* 7 2 CHAR l imiti 

BOTTOM CENTERLINE HEATING TO GENERIC ORBITER AT SaS'b-SeL B 

INTERACTIVE INPUT FOR CASE 2 


IS CASE 2 ATMOSPHERE DATA SANE AS FOR CASE 1 ? 

T tv 


IS CASE 2 WIND TUNNEL DATA SAME AS FOR CASE 1 ? 

Ttw 


DO YOU WANT TO RUN A HEATING INDICATOR ? 
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J£ CASE 2 HEAT TRANSFER DATA SAME AS FOR CASF 1 ? 

HEAT TRANSFER METHOD 

OPTIONS t. 

2 . 

3. 

4. 

5. 

6 . 

7. 

8 . 

9. 

10 . 

11 . 

OPTION SELECTED ? 

4.0 

SHOULD RAREFIED FLOW HEATING BE INCLUDED ? 

NO 

IS THE HEAT TRANSFER OPTION CORRECT ? 


HEMISPHERE STAGNATION POINT 
CATO/ JOHNSON SWEPT CYLINDER 

ENTHALPY FLAT a ATE METHOD 
/ L P ULD,NG ^« FLAT FLATE METHOD 
BOEING RHO-MU FLAT a ATE METHOD 

iii!^ T ^ri E T«AT H ^ SB,ERED,STRIBUTI0N 

FLAP RE ATTACH f€ NT HEATING 
FIN-aATE PEAK INTERFERENCE HEATING 


4. ECKERT/SPAULDING-G • FLAT a ATE METHOD 


RUNNING LENGTH (FT) ? 

•333 2 

SURFACE DISTANCE TO START OF TURBULENT B.L. 

THIS DISTANCE IS SUBTRACTED FROM THE RUNNING LENGTH 
FOR TURBULENT HEATING CALCULATIONS. 

DESIRED LENGTH ? 

0.0 

IS AN AUTOMATIC VIRTUAL ORIGIN CORRECTION DESIRED ? 
NO 


TURBULENT MANGLER FACTOR ? 
.6667 

LAMINAR MANGLER FACTOR ? 

1.0 


REYNOLDS-ANALOGY FACTOR 

0. COLBURN 

1. VON KARMAN 

DESIRED FACTOR ? 

1 . 

ANY CHANGES ? 

NO 


* * * WALL conditions * # # 

WAM TFMPFFiATIJRf (DEG F) ? 

WALL EMI SS IV ITT ? 

0.8 


* * * CONTINUATION OPTION * « # 


DO YOU WANT TO PROVIDE 


AN INITIAL HEATING LOAD GT 


0.0 (BTU/5Q.FT) ? 



DO YOU WANT TO USE 
NO 


A HEAT TRANSFER FtJLTlPL I CAT I ON METHOD ? 



YES 0 ^ 2 TRANSITI0N DATA THE SAME FOP CASE 1 ? 


DO YOU WANT CROSS FLOW ADJUSTMENT OFT ION ? 



IS CASE 2 FLOWFIELD DATA THE SAME AS FOR CASE 1 ? 
FLOWFIELD AND LOCAL PRESSURE OPTIONS 


“ FLOWFIELD 

“1. FLOWFIELD TYPE NOT NEEDED -1 

1. SHARP WEDGE SHOCK ANGLE , 

2. SHARP 00 NE SHOCK ANGLE 2* 

J* 2?k[? U f AND N0RMA L SHOCK (90 DEG) 3 ! 
4. PARALLEL SHOCK (PRES NOT NEEDED) 4 , 

5. 

6 . 


- * - - - PRESSURE - - - 

PRESSURE TYPE NOT NEEDED 

INPUT CP VS MACH NO. TABLE 

TANGENT WEDGE PRESSURE 

TANGENT CONE PRESSURE 

OBLIQUE SURFACE PRESSURE 

MODIFIED NEWTONIAN 

PRANDTL-MEYER EXP. (FF NOT NEEDED) 


(note! TO l SIGN|^ r £Nrof N CASE R USE IT ? Ft, mH™ DaTA /mLeS * * * * 
fUWFIELD.DEL ANGLE, PRESSURE, DEL ^NGLE ? ” B0 ™ ” AN ° PRESS ' 0f TIONS.) 

2 . 0 , 0 . 3 , 2 . 0 , 0.3 
2 

“ 1 . 0 , 0 . 0 ,- 1 . 0 , 0.0 


SET1 

FF 

o 


SHRP-CONE 

TAN-WEDGE 


TOTAL EFFECTIVE ANGLE 


ALPHA + 
ALPHA + 


0.300 

0.300 


ANY CHANGES ? 
NO 



L^TOU WANT TO USE TIME DEPENDENT GEOMETRY ? 


IS CASE 
NO 


2 time dependent geometry data the same as for case 


1 ? 


time dependent geometry 

NUffcER OF TIME DEPENDENT ENTRIES ? (MAX*10> 

KM!o"5692 0 !i. , 0 T> ' LENGTH< ' :T1,SLOFt W S«EP<DEG> 

2.0. 0. 0..3152.0.0 

3. 0. 0.0.. 2277.0.0 

4. 0. 0.0. .2277.0.0 
ANY CHANGES ? 

NO 


CONTROL FLAGS 


YOU ARE COMPLETING INPUT FOR CASE 2 

WAT IS THE BOOT POINT NUHiER FOR CASE 2 ? 

SHOULD LANMIN CREATE AN OUTPUT FILE FOR CASE 2 ? 


0. 

1. 

2 . 


nrT.MPH B WNIR0L OPTIONS 

detailed printout 
detailed plus summary printout 

SUMMARY PRINTOUT UT 


«SJiN , f.«W? , S! SELECTED ® E * TED E ' THER 

OPTION SELECTED ? 

1.0 


OUTPUT UNITS OPTIONS 

0. ENGLISH 

1. METRIC 

OPTION SELECTED ? 

0 . 


IS CASE 
AS CASE 
NO 


' ALONG THE SAME STREAM I NE 


1 . 

2 . 

3. 

4. 

5. 

6 . 


INPUT CONTROL FLAG 

X'S T "U. T, «» and traj. 

END OF INPUT (LAST CASE) 

NEN T?TL E^AND^T^M^NG?** I N^T^ AL^CASE^DATA^nwrw js E MSE * 
S^W.'SlSrr W-CASE DATA. 

® - - SR 


OPTION SELECTED ? 

1 .0 


***** INPUT COMPLETE FOR CASE 2 «... 

OO^TOU NANT TO MAKE ANY MODIFICATIONS TO CASE 2 ? 
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modification sections 

1. TIMING PARAMETERS 

2. TRAJECTORY DATA 

3. ATMOSPHERE DATA 

J. FLOWFIELD AND PRESSURE DATA 

5. CROSSFLOW DATA 

6. TRANSITION CRITERIA 

7. HEAT TRANSFER OPTION 

8. HEATING MULTIPLIERS 

9. GEOMETRY DATA 

10. CONTROL PARAMETERS 

11. HEATING INDICATOR 

12. OR CHANGE A SPECIFIC VARIABLE IN W ARRAY 

SECTION TO BE MODIFIED ? 

9. 


is CASE 2 TIME DEPENDENT GEOfETRY DATA THE SAME AS FOR CASE 


time dependent GEO^TRY 


NUteER OF TIME DEPENDENT ENTRIES ? (MAX*10) 

1 l ^o S n C> ^?? l i ( r r),LENGTH<FT) » SL0PE 0R SWEEP (DEG) 

J #u,u.u, #3532,0.0 

2. 0. 0.0. .1943. 0.0 

3. 0. 0.0.. 1448.0.0 

4.0. 0.0.. 1448.0.0 
ANY CHANGES ? 

NO 

DO YOU WISH TO MODIFY ANY OTHER SECTIONS FOR CASE 2 ? 


***** INPUT OOMPuETE FOR CASE 2 ***** 

DO YOU WANT TO MAKE ANY MODIFICATIONS TO CASE 2 ? 


INTERACTIVE INPUT FOR CASE 3 





IS CASE 3 WIND TUNNEL DATA SAME AS FOR CASE 2 ? 


DO YOU WANT TO RUN A HEATING INDICATOR ? 
NO 


IS CASE 3 HEAT TRANSFER DATA SAME AS FOR CASE 2 ? 


HEAT TRANSFER KETHOD 


OPTIONS 


1 . 

2 , 

3. 

4. 

5. 

6 . 

7. 

8 . 

9. 

10 . 

11 . 


HEMISPHERE STAGNATION POINT 
CATO/ JOHNSON SWEPT CYLINDER 
ECKERT REF. ENTHALPY FLAT PLATE METHOD 
ECKERT/ SPAULO I NG-CH I FLAT PLATE IETHOD 
BOEING RHO-MU FLAT PLATE METHOD 
BECKWITH/ GALLAGHER SWEPT CYLINDER NETHOD 
BOEING RHO-MU SWEPT CYLINDER KETHOD 
LEES/DETRA-HIDALGO HEMISPHERE DISTRIBUTION 
LEES IDE ORB I TER HEATING 
FLAP REATTACH KENT HEATING 
FIN-PLATE PEAK INTERFERENCE HEATING 


OPTION SELECTFD ? 

4.0 

SHOULD RAREFIED FLOW HEATING BE INCLUDED ? 
NO 

IS THE HEAT TRANSFER OPTION CORRECT ? 


4. ECKERT/SPAULD I NG-CH I FLAT PLATE KETHOD 


RUNNING LENGTH (FT) ? 
.5254 


SURFACE DISTANCE TO START OF TURBULENT B.L. 

THIS DISTANCE IS SUBTRACTED FROM THE RUNNING LENGTH 
FOR TURBULENT HEATING CALCULATIONS. 

DESIRED LENGTH ? 

0.0 


IS AN AUTOMATIC VIRTUAL ORIGIN CORRECTION DESIRED ? 
NO 


TURBULENT MANGLER FACTOR ? 
.6667 

LAMINAR MANGLER FACTOR ? 

1.0 


REYNOLDS-ANALOGY FACTOR 

0. COLBURN 

1. VON KARMAN 

DESIRED FACTOR ? 

1 . 

ANY CHANGES ? 

NO 


* * * WALL CONDITIONS * * * 

WALL TEMPERATURE (DEG F) ? 

0.0 

WALL EMISSIV ITY ? 

0.8 


* * * CONTINUATION OPTION » »« 

DO YOU WANT TO PROVIDE AN INITIAL HEATING LOAD GT 0.0 (BTU/SQ.FT) ? 
NO 


DO YOU WANT TO USE A HEAT TRANSFER KJLTIPLICATION METHOD ? 
NO 


IS CASE 3 TRANSITION DATA THE SAME FOR CASE 2 ? 
YES 


DO YOU WANT CROSS FLOW ADJUSTMENT OPTION ? 
NO 


JS CASE 3 FLOWFIELD DATA THE SAME AS FOR CASE 2 ? 
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FLOWFIELD AND LOCAL PRESSURE OPTIONS 


! * FLOWFIELD 

* FLOWFIELD TYPE NOT NEEDED 7 “ “ ‘ PRESSURE 

1. SHARP WEDGE SHOCK ANGLE “!* ^P^SURE TYPE NOT NEEDED 

i: awxaSc « 0EG , : wasar 

4. PARALLEL SHOW (PRES WT NEEDED) ^ ?! g SaS^re 

5. MODIFIED NEWTONIAN 

# # # # 6# PRANDTL-MEYER EXP. (FF NOT NEEDED) 

(NOTE: TO S^GNI^ENDy^S^USE -? ^pS^I ED DELTA ^LES * # # # 
FLOWFIELD, DEL ANGLE, PRESSURE, DEL^ANGLE F0R B0TH FF AND PRESS * OPTIONS.) 


2 . 0 , 0 . 0 , 2 . 0 , 0.0 

2 

- 1 . 0 , 0 . 0 ,- 1 . 0 , 0.0 


SET1 

FF 

P 


SHRP-CONE 

TAN-WEDGE 


TOTAL EFFECTIVE ANGLE 

ALPHA + 0.000 

ALPHA + 0.000 


ANY CHANGES ? 
NO 


OO^YOU WANT TO USE TIME DEPENDENT GEOteTRY ? 


IS CASE 3 TIME DEPENDENT GEOtCTRT DATA THE SAME AS FOR CASC 


2 ? 


TIME DEPENDENT GEOI€TRY 


NUfrBER OF TIME DEPENDENT ENTRIES ? (MAX«10) 

J .V ’ ' LENGTN <FT ’ • SL 0PL OR SWEEP (DEG) 

2. 0. 0.0. .31 52)o!o, 

3.0. 0.0.. 2277. 0.0 

4. 0. 0.0. .2277.0.0 
ANY CHANGES ? 

NO 


CONTROL FLAGS 


YOU ARE COMPLETING INPUT FOR CASE 3 

300 T ,S TWE B0DY P0,NT NW ® ER F0R CASE 3 7 
SHOULD LANMIN CREATE AN OUTPUT FILE FOR CASE 3 ? 


0. 

1. 

2 . 


PRINT OONTROL OPTIONS 
DETAILED PRINTOUT 
DETAILED PLUS SUMMARY PRINTOUT 
SUMMARY PRINTOUT 


NOTE: IF AN 
OPTION 1. OR 


FILE IS TO BE CREATED EITHER 
2. MUST BE SELECTED 


OPTION SELECTED ? 

1.0 


OUTPUT UNITS OPTIONS 

0. ENGLISH 

1. METRIC 

OPTION SELECTED ? 

0. 


IS CASE 4 ALONG THE SAME STREAM I NE 
AS CASE 3 ? 

NO 


INPUT OONTROL FLAG 


1 . 

2 . 

3. 

4. 

5. 

6 . 


0 N A?A C K^OUS S iA N S E T ' TLE ' TIMING 

END OF INPUT (LAST CASE) 

NEW Trn E p F Shh°!?m?r N6 TRAJ * DATA FR0M PREVIOUS CASE. 

NEW CASE E F^?nIi M |Kmr !Ji T,AL CASE DATA UN °HANGED. 

TITLE - TIH,NS '™ J - **> «« DATA. 

SAMF AC !1! riSS * ER0 ALL Ms£ DATA FR0M PREVIOUS CASE 
SAME AS (3) EXCEPT INITIALIZE ZERO ALL TIMING AND CASE DATA 


OPTION SELECTED ? 

2.0 


4 

-■* i 

’! ***** INPUT COMPLETE FOR CASE 3 ***** 

DO YOU WANT TO MAKE ANY MODIFICATIONS TO CASE 3 ? 

NO 

„r : 

44 

*«««**»« OUTPUT FILE COMPLETE ******** 

1 > 



PREMI N OUTPUT FILE 
WIND TUNNEL CASE 


‘ i 


BOTTOM CENTERLINE HEATING 
1.000 
0.000 
0.000 

4.000 

1.000 
2.000 

3.000 

4.000 

4.0000 

1 .0000 


36.0000 

2.0000 

96.9000 

4.0000 


10 1.0000 

11 

130.2097000 


14 3.2000 

15 

160.6667000 


20 3.6000 


230.8000000 


27 1 .0000 

31 

261 1 .0000 

262 


11 


453 95.0000 454 
504 12.6720 560 
564 4.0000 

5810.2097000 
5820.1239000 

5830.0858000 

5840.0858000 
611 100.0000 642 

10 1.0000 
130.3532000 
14 3.2000 

160.6667000 
20 3.6000 

230.8000000 
27 1 .0000 

370.3000000 

380.3000000 
261 1.0000 262 
453 95.0000 454 
504 12.6720 560 
564 4.0000 
5810.3532000 
5820.1943000 

5830.1448000 

5840.1448000 
611 200.0000 642 


1 .0000 

4.0000 


15 1.0000 


TO GENERIC ORB I TER AT MACH 8-TUNNEL B 
1.000 4.000 

0.000 0.000 

1.000 


0.000 

0.000 

0.000 

0.000 


3816.000 

3821 .000 

3821.000 

3859.000 


20.000 

29.860 

39.980 

39.980 


32 

319 

501 

561 


15.0000 37 
1.0000 451 
8.0640 502 
1.0000 562 


2.8000 38 2.8000 
94.7000 452 95.0000 
8.2080 503 8.2080 
2.0000 563 3.0000 


643 1.0000 647 1.0000 641 1.0000 


31 36.0000 32 15.0000 


2.0000 319 1.0000 451 

96.9000 501 8.0640 502 

4.0000 561 1.0000 562 


94.7000 452 
8.2080 503 
2.0000 563 


95.0000 

8.2080 

3.0000 


1.0000 643 1.0000 647 2.0000 641 1.0000 
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11 

15 


31 


10 1.0000 
130.5254000 
14 3.2000 

160.6667000 
20 3.6000 

230.8000000 
27 1.0000 

319 1.0000 451 

501 8.0640 50 J 

561 1.0000 565 

5810.5692000 
5820.3152000 

5830.2277000 

5840.2277000 
611 300.0000 642 


4.0000 

1.0000 


36.0000 32 
94.7000 452 
8.2080 503 
2.0000 563 


1.0000 643 


15.0000 261 

95.0000 453 
8.2080 504 

3.0000 564 


1.0000 647 


1.0000 262 2 
95.0000 454 96 
12.6720 560 4 

4.0000 


3.0000 641 2. 
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TINE VELOCITY F.S.TE*. F.S.PRESS. 
LOO 3816. 94.70 .80640*001 

2.00 1021. 95.00 .82000*001 

3.00 3021. 95.00 .82060*001 

4.00 3859. 96.90 .12672*002 
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PSEMIN CODE LISTING 


oo ooooooooon 


^4 


PRtMIN CODE 

LANGLEY MINIVER PREPROCESSOR CODE 
REMTECH, INC. 1083 

BY: 

C. ENGEL 
C. SCHMITZ 

PH. 205-536-8581 

PROGRAM MAIN 


COMMON/ WARR AY/ W (700) 

COMMON/ W WAR AY/ WW ( 7 00 ) 

COMMON/ WWWARY/ WWW (700) 

COMMON/UNIT/ 1 IN, I OUT 
C0MM0N/TITLE/HTL1 

COMMON/POOEFF/TMCP ( 50 ) ,TCPM( 50 ) , NCPMT 
COMMON/M I SC/NC, JFKS, TFL AG, MFL AG, ANSI 
c COMMON/ I NP/N I , ND, INSERT, NCPMT S , TMCPS ( 50 ) , TCPMS ( 50 ) , JFL AG, (ADD 

CHARACTER *20 FNAM6 
CHARACTER*72 TITL1 
INTEGER ANS,TFLAG, ANSI, FLAG 
C I IN - INTERACTIVE INPUT UNIT NUMBER 
C I OUT - INTERACTIVE OUTPUT UNIT NUNBER 
IIN*1 
IOUT a 1 
C 

C MFL AG = 0 ENGLISH UNITS 

C MFL AG - 1 METRIC UNITS 

C FLAG ■ 1 CREATE NEW OUTPUT FILE 

C FLAG - 1 MODIFY AN EXISTING OUTPUT FILE 

C JFLAG = 1 DELETED LAST CASE 

C INSERT « 1 CURRENT CASE AN INSERTION CASE 

C TFL AG =1 TITLE EXISTS FOR CURRENT CASE 
C FNAM6 - NAME OF INPUT FILE TO BE MODIFIED 

C r I CAS£ R NUlvBER UT °° NTR0L PARAMETER for PREVIOUS CASE 

C ND - NUMBER OF CASES DELETED 

C Nl - NUNBER OF CASES INSERTED 

C TITL1 - TITLE FOR A SET OF CASES 

C NCPMT - NUMBER OF ENTRIES IN PRESSURE COEFFICIENT TABLE 

C NCPMTS - SAVED VALUE OF NCPMT 
C TMCP - MAGI NUNBERS 

C TMCPS - SAVED VALUES OF TMCP ARRAY 
C TCPM - CP ARRAY 

C TCPMS - SAVED VALUES OF TCPM ARRAY 
C W - WORKING ARRAY OF W VALUE TO BE MODIFIED 

C WW - W ARRAY MODIFIED TO DESIRED UNITS FOR OUTPUT 

c (INPUT ARRAY FOR LANMIN) 

INSERTS 

NC=0 



C-l 


\ 


Nl a Q 

NIX) 

MU Arx 1 
JFKSH) 

W(641 )-0.,J 
WRI1F (| OUT, 30) 

30 I ORMAH1H1,' INTERACT IVI INPlJi 
C OKI A 1 1 OK MOD I F Y (JUTIU'I FILE 


10 FARC MINIVER - LANMIN* ) 


35 

31, 


WRI1I (I OUT ,36) 

I GRMAT ( ///, 1 0X, * 01 1 ITU | I II I METHOD*//, 
t IX,' 1 . cm All A Nrw OUT HIT Mil*,/, 

1 IX, '2. MODIFY AN ITX I SI I NIT DONUT f 1 1 f * // 
i IX, 'OPT ION Ml FOTFD ?•) ' 

Nf AIK I IN, *, ERR* 35, 1, NO- 0000 ) 1 1 AO 
IF (FLAG. tip. 1 )G0 TO 20 
C NAME FILL FOR MODIFICATION 
10 WRITE (I OUT, 12) 

12 FORMAT!//, IX, 'WHAT IS THE NAME OF THE INPUT FILE ?») 
READ! I IN, 13,ERR S 1 0,END=9999)FNAM6 

13 FORMAT (A20) 

OPENfUNIT-^FILE^FNAMS, STATUS'- 'OLD' ,ERR a 8888) 
IF(MFLAG.EQ.1 ) WRITE! I0UT,25) 

25 FORMAT!/, IX, ‘INPUT FILE IS IN ENGLISH UNITS.', /.IX 
X 'ALL CHANGES MUST BE IN ENGLISH UNITS.'./) 

MFLAG=0 

33 IFOADD.EQ. 1 )G0 TO 37 

I F (FLAG.EQ.2)CALL INPUT! 1) 

IF! INSERT. EQ. 1 ) GO TO 37 
I F ( JFLAG.EQ. 1 )G0 TO 3 

34 IF(FLAG.EQ.2)G0 TO 200 
- CHOOSE UNITS TO BE USED 

WRITE! I0UT,21) 

FORMAT (IX, 'DO YOU WANT TO INPUT DATA IN ENGLISH OR ' 
$ 'METRIC ?') 

READ! I IN,45,ERR=20,END=9999)ANS 
MFLAG=0 

IF(ANS.EQ.1HM)MFLAG=1 
CONTINUE 
NONC+1 
W(647)=NC 

INPUT DATA FOR CASE 


20 

21 


37 


3 

38 

1 


39 

40 


WRITE! I OUT, 38 )NC 
FORMAT!/, 1X,» INTERACTIVE INPUT FOR CASE 
I F( JFKS.GT.OJGO TO (2,9900,1,1,2,1 >J*KS 
CALL TIMING 

I F ( JFKS.EQ.3.QR. JFKS.EO,g)G0 TO 2 
CALL TRAJ 
CALL ATMS 

(F(W(IO).EQ.I JCALL WNDTUN 
IF(W(10).Ep.2.)CALL ATMDTA 
WRITE (I OUT, 40) 

FORMAT !1H1,' DO YOU WANT TO RUN 


13) 


READ! I IN,45,ERR=39,END=9999)ANST 
FORMAT (AH 

IF (ANSI .EQ.1 HN)G0 TO 50 


A HEATING INDICATOR ?•) 


C-2 


N 

CALL HEAT IN 
GO TO 9000 
5>'J CALL HTRM7D 
WRITE (I OUT, 55) 

• r; FORMAT (1H1,» DO YOU WANT TO USE A HEAT TRANSFER ♦, 

$ 'MULTIPLICATION METHOD ?•) 

READ( I IN,45,ERR"54, END-9999) ANS 
IF(ANS.FQ.1HY)CALL HTMULT 
CALL TRANS 

69 WRITEd0UT,70) 

70 FORMAT ( 1H1 , 'DO YOU WANT CROSS FLOW ADJUSTMENT OPTION ?') 

Rf.AD( I IN,45,ERR s 69,END*9999)ANS 
I F ( ANS. EQ. 1 HY ) CAL L CROSS 
CALL FLOW(I) 

89 WRITEdOUT,90) 

90 FORMA7(1H1,'DO YOU WANT TO USE TIME DEPENDENT GEOMETRY ?') 
READd IN,45,ERR B 89,END s 9999)ANS 
I F (ANS. EQ. 1 HY )CALL TDGEOM 

9000 CONTINUE 
CALL OONTRL 
200 CONTINUE 

WR ITE ( I OUT, 1 00 ) NC 

100 F0RMAT(//,1X, ' ««*** INPUT COMPLETE FOR CASE M3. • ******* 

?10 WR I TE ( I OUT, 220 ) NC 

2‘^0 F0RMAT(/,1X, 'DO YOU WANT TO MAKE ANY MODIFICATIONS TO ', 

$ 'CASE ',13,' ?•) 

READd IN,45,ERR*210,END=9999)ANS 
IF(ANS.EQ.1HY)CALL MODIFY 
IF(ANS.EQ.1HV)G0 TO 200 
CALL UNITS 
CALL OUTPUT 

IF(W(641).EQ.2.)00 TO 9900 

C ZERO W ARRAY ACCORDING TO PROGRAM INPUT OONTROL PARAMETER 
C OF LAST CASE 
CALL INPUT(2) 

• F0NSERT.EQ.1 )G0 TO 33 
IF(JFLAG.E0.1)FLAG=1 
IF(FLAG.EQ.2)G0 TO 33 
GO TO 37 

9900 WRITEdOUT, 110) 

110 FORMAT <//, 1 X, ' **»#»#** OUTPUT FILE OOMPLETE ********* ) 

CALL EXIT 

8888 CONTINUE 
WRITE(I0UT,8889)FNAM6 

8889 FORMAT (/, IX, 'CANNOT OPEN FILE »,A20,/) 

CO TO 10 

9999 CONTINUE 
CALL EXIT 


o o 




SUBROUTINE INPUT(IB) 


C THIS ROUTINE HAS 2 FUNCTIONS 

£ !! ' 1 ?“!I INE R£A0S IN " ««AY DATA from old run file 

c ’ param^r^las/case 0 ^ 13 T0 PR06R * m INPUT °° ntrol 

c JFK - PROGRAM INPUT OQNTRQL PARAMETER 
C JFCP - PRESSURE COEFFICIENT INPUT OPTION 
COMMON/ WARRAY/W ( 7 00 ) 

COMMON/UNIT/ 1 IN, IOUT 
OOMMON/PCOEFF/TMCP ( 50 ) , TCPM ( 50 ) , NCPMT 
COMMON/T I TL E/T I TL 1 

COMMON/MI SC/NC, JFKS, TFL AG, MFLAG, ANSI 
c COMMON/ 1 NP/N I , ND, I N SERT, NCPMTS, TMCPS ( 50 ) , T CPMS ( 50 ) , J FL AG, IAOD 

CHARACTER*72 TITL1 
INTEGER AN$,TFLAG,ANS1 
DIMENSION L(5),X(5) 

CHECK FOR IB * 2 
IF( IB.EQ.2)G0 TO 5 

C CHECK IF LAST CASE WAS AN INSERTED CASE 
IF( INSERT.EQ. 1 )GO TO 400 
NC1*NC+1 

2 WRITE( IOUT,3)NC1 

3 Y ° U W,SH T ° ,NSERT A CASE BEF ° RE CASE 

READd IN,4,ERR=2,END s 9999)ANSS 
IF(ANSS.EQ.1HY)GO TO 300 
4 FORMAT (AD 
5 CONTINUE 

JFK*W(641 1+.0001 
I F ( IB.EQ.1 ) JFK® JFKS 
JFKS*W(641 >+.0001 
IF(JFK.EQ.1.AND. IB.EQ.D60 TO 160 
IF(JFK.EQ.3)G0 TO 130 
IF(JFK.EQ.4)G0 TO 100 
IF(JFK.EQ.5)G0 TO 120 
IF( JFK.EQ.6)G0 TO 110 
IF( IB.EQ.2)G0 TO 130 

100 CONTINUE ~ ” W ARRAY 

DO 1000 J=1 ,700 
1000 W(J)»0.0 

C no~~Do"noo'j~D8 2ER0 TIMING PARAMETER s and print control 

1100 W(J)=0.0 ' 

120 DO 1200 J s 9,49 ZER0 CASE DATA 

1200 W(J)=0.0 

DO 1300 J=201,210 
1300 W( J)=0.0 

DO 1400 J=261,650 
1400 W(J)=0.0 

IF( JFK.EQ.5.AND. IB.EQ.1 )G0 TO 160 


C-4 


130 


C — 

135 

C 

140 

C — 
150 
155 


C — 

160 
C — 

170 

171 
175 

180 


IF(JFK.EQ.6)G0 TO 130 
IF(IB.EQ.2)G0 TO 130 
JFK*0 
CONTINUE 

IF ( IB.EQ.2)W(641 )*FLQAT( JFK) 

IF( IB.EQ.2)RETURN 

TFLAG*1 


READ(9,135)TITL1 
FORMAT (A7 2) 


TITLE 


TIM,NG PARAMETER m °° ntro1 

FORMAT (3F 20. 6) 

I F( JFK. EQ. 3. OR. JFK. EQ. 6 )G0 TO 160 

REAdI9~55)W(W) R ° F TIME DEPENDENT T/ ' BLE ENTRIES FOR TRAJ. 

FORMAT (3F20. 6, 2F 10.6) 

N=W<50)+.0001 

NT*50+N 


READ(9,155) (W(K),W(K+50),W(K+100) 
W(641 )*0.0 


TRAJ DATA 

,W(K+160),W(K+600),K*51,NT) 


READ (9, 175) (L(J),X(J),J*1,5) 
DO 171 J=1,5 

IF(L(J).GT.O)W(L(J))=X(J) 
F0RMAT(5( 13, FI 0.6 ) ) 

I F (W(641 ) )170, 170, 1 80 
CONTINUE 


CASE DATA 


IF(W(641),NE.2.0)G0 TO 185 

181 WRITE (I OUT, 182) APPENDING CASES - 

182 $ F8 ^RENT 1 END D CASE U ?M NT T ° ADD ADD,TI0NAL CASES AFTER THE », 

READ( I IN»4,ERR=1 81 ,END=9999)ANS 
IF(ANS.EQ.IHY) I ADD® 1 
185 CONTINUE 

W ( 647 ) * W ( 647 ) -FLOAT ( ND ) +FLOAT ( N I ) 

JFK*W(641 )+.0001 
NC*W(647)+.0001 
JFCP»W( 649) +.0001 
IF(JFCP.LE.O)GO TO 210 
READ ( 9, 1 90 ) NCPMT 
190 FORMAT (13) 

RE AD ( 9, 200 ) (7MCP ( J ) , TCPM ( J ) , J= 1 , NCPMT ) 

200 FORMAT (2F 10.6) 

210 CONTINUE 


C — 
220 
230 


IF(W(641 ) .EQ. 2. )CL0SE(UNIT*9, STATUS* ’KEEP’ ) 

WRITE( I0UT,230)NC DELETED CASES 

FORMAT ( IX, *D0 YOU WISH TO DELETE CASE »,I3 * ?*) 
READ( I IN,4,ERR*220,END s 9999)ANS 
IF(ANS.EQ.1HN)G0 TO 290 
IF(W(641).EQ.2.)JFLAG=1 
IF(JFLAG.EQ.1)G0 TO 290 


240 WRITE(I0UT,250)NC,NC+1,NC 

250 FORMAT (IX, 'CASE DELETED. CASE ',13,' REDEFINED », 

$ 'AS CASE ',13) 

ND=ND+1 

IF(NC.NE,1 )60 TO 5 
IF (JFK.EQ.4)G0 TO 5 
260 WRITE (I OUT, 270) 

270 F0RMAT( IX, 'SINCE YOU WISHED TO DELETE CASE 1 DO YOU ALSO » 
S 'WANT TO DELETE THE *,/, IX, 'TITLE, TIMING, AND TRAJ. THAT • 
S 'CORRESPONDS TO CASE 1 ?•) 

READd IN,4,ERR=260,END B 9999)ANS 
I F(ANS.EQ. 1 HY)GO TO 5 
60 TO 160 
290 CONTINUE 
RETURN 


300 


CONTINUE 
ANSS* *N» 

INSERTS 

W(641 ) “FLOAT (JFKS) 
60 TO 290 


INSERTED CASES 


400 


CONTINUE 
INSERTS 
NI=NI+1 
NONC+1 

W ( 647 ) “FLOAT ( NC ) 

WRITEt I0UT,3)NC 
READ ( 1 1 N, 4 , ERR=430, END*9999) ANSS 
I F ( ANSS. EQ. 1 HY )60 TO 300 
60 TO 5 
9999 CONTINUE 
CALL EXIT 


LAST CASE INSERTED - 


430 



1 


SUBROUTINE ATMDTA 
C 

c R0 £il0N F 2 R ,NPUTING A ™° Sft,ERE °ATA FOR ATMOSPHERE 

OOMMON/WARRAY/W ( 700 ) 

COMMON/UN IT/ 1 IN, I OUT 

c COMMON/M I SC/NC, JFKS, TFLAG, MFL AG, ANSI 

INTEGER ANS, TFLAG, ANSI 

C MAY USE DATA FROM PREVIOUS CASE IF JFKS «1 OR 3 
IFfJFKS.NE.I.AND.JFKS.NE.J.Gj TO 9 3 

NUM1 «NOi 

WR I TE ( I OUT, 2 ) NC, NCM1 

* f ?“e’::? 3 “%;- iv atms - mta «« « •. 

READ ( 1 1 N, 1 0O, ERR*1 , END® 9999 ) ANS 
I F (ANS. EQ. 1 HY JRETURN 
CONTINUE 
WRITE (I OUT, 20) 

FORMAT(///^» INPUT ATMOSPHERIC DATA',//. 

j ? ’ 

WRITE(I0UT,40) 

FORMAT ( IX, * HOW MANY ALTITUDE VALUES WILL BP input ->i\ 

REAO(IIN,.,ERR=30,EN0.9999)W(400) ' NPUT 

NALT « NUNBER OF ALTITUDE ENTRIES 
NALT=W(400)+.001 
IF(MFLAG.EQ.0)WRITE(I0UT,50) 

IF(MFLAG.EQ.1 )WRITE( I0UT,51 ) 

FORMAT( IX, * ALT I TUDE (FT ) ,T- I NF (R ) , P- | NF ( LB/SFT ) n 

DO^TOOO ^ M ^NALT^ (M) # T " 1 NF (K } » p “ 1 NF (NEWTON/SQ. M) » ) 
WRITE (I OUT, 70) I 
FORMATdX, 12) 

1000 continue' * ,ERR ’ 60 ' ENO ' 9999,I,<40 0'* 1 >.««0*I).WSOO*I ) 
80 WRITE (IOUT, 90) 

D??U?Tf /,1X,,ARE ALL INPUTS CORRECT ?*) 

fo^tIai )° 0#ERRa80,ENDa9999)ANS 

IF(ANS. EQ. 1HN)G0 TO 10 
RETURN 

9999 CONTINUE 
CALL EXIT 
END 


10 

20 


30 

40 


50 

51 

60 

70 


90 

100 


C-7 


o o 


SUBROUTINE ATMS 


routine that chooses atmosphere option to be used 
COMMQN/WARRAY/W( 700) 

COMMON/ UN IT/ 1 IN, I OUT 

COMMON/M I SC/NC, JFKS , TFL AG, MFL AG, ANSI 


CHARACTERMO CHAR(5) 

INTEGER ANS, TFL AG, ANSI 
C TITLES FOR ATMOSPHERE OPTIONS 

DATA QHAR/'I 962 U.S. STANDARD ATMOSPHERE’, 

$ ’WIND TUNNEL OPTION’, 

$ ’INPUT ATMOSPHERIC DATA (ALT, T- I NF, P- I NF ) ' , 

$ *1963 PATRICK AIR FORCE BASE ATMOSPHERE’, 

$ '1971 VANDENBERG REFERENCE ATMOSPHERE'/ 

C USE ATMS, DATA FROM PREVIOUS CASE IF JFKS * 1 OR 3 
IF(JFKS.NE.1. AND. JFKS. NE.3)G0 TO 9 
NCM1=NC-1 


1 

2 


9 

10 

20 


30 

40 

50 

60 


WR I TE ( I OUT, 2 ) NC, NCM1 

S F< *F^ T CASE # • 1 ?3°^ S ? ’ > ^ l3# * ATM0SF * ERE data SAME AS ’, 

READ (I IN,60,ERR=1 ,END=9999)ANS 
I F (ANS. EQ. 1 HY )RETURN 
CONTINUE 


WRITE(I0UT,20) (l,CHAR( I ), 1=1,5) 

FORMAT (/// ,5X, * ATMOSPHERE DATA’ ,//, 1 X, ’ OPT I ONS 
S II,'. *,A40,/ ,4( 1 IX, II ,* . »,A40,/), 

J /, 1 X, » OPT I ON SELECTED ?») 


READ( 1 1 N, ”, ERR®1 0, END*9999) W( 1 0) 
J*IFIX(W( 10)) 

WRITE(I0UT,30)CHAR(J) 

FORMAT (/,1X,A40) 

WRITE( I0UT,50) 

FORMAT (/, IX,’ IS THIS OPTION CORRECT ?•) 
READ( 1 1 N, 60, ERR=40, END= 9999) ANS 
FORMAT (A1 ) 


IF(ANS.EQ.1HN)G0 TO 10 

IF(W(10).LT.4.)W(10)=W(10)-1. 

RETURN 

9999 CONTINUE 
CALL EXIT 
END 


no 


SUBROUTINE OONTRL 


ROUTINE FOR SETT I NG CONTROL FLAGS 
COMMON/ WARR AY/ W ( 700 ) 

COMMON/UN IT/ 1 IN, I OUT 
COMMON/MI SC/NC, JFKS, TFL AG, MFLAG, ANSI 
COMMON/ STRM/ I STRM, NNC 
C 

INTEGER ANS, TFL AG, ANSI 
C 

ISTRM-0 
NOW (647) 

10 WRITE( I0UT,20) 

20 FORMAT <1 HI, 1 0X, 'OONTROL FLAGS*,//) 

30 CONTINUE 
50 WRITE(IOUT,60)NC 

60 F0RMAT(1X,'Y0U ARE COMPLETING INPUT FOR CASE »,I3> 

65 FORMAT (A1 ) 

70 WRITE (I0UT,80)NC 

80 FORMAT (/, 1 X, * WHAT IS THE BODY POINT NUNBER FOR CASE », 

$ 13,* ?*) 

READ ( 1 1 N, », ERR*70, END* 9999 ) W ( 6 1 1 ) 

90 WR ITE ( I OUT, 1 00) NC 

100 FORMAT ( IX, • SHOULD LANMIN CREATE AN OUTPUT FILE FOR CASE » 

$ ,I3,» ?•) 

READd IN,65,ERR*90,END=9999)ANS 
IF(ANS.EQ.1HY)G0 TO 110 
W(643)*0.0 
GO TO 120 
110 W(643)*1.0 

C CHECK FOR HEATING INDICATOR OPTION 
120 IF(ANS1.EQ.1HY)G0 TO 140 /* HEATING INDICATOR W(642)*2.0 

C PRINT OPTION 

WRITE(IOUT,130) 

130 F0RMAT(/,1 OX, ‘PRINT CONTROL OPTIONS*,/, 

$ IX, *0. DETAILED "INTOUT* ,/, 

$ IX, *1. DETAILED h-.i SUMMARY PRINTOUT*,/, 

$ IX, *2. SUMMARY PRINTOUT*,//, 

$ IX, 'NOTE: IF AN OUTPUT FILE IS TO BE CREATED EITHER *,/, 

$ IX, 'OPTION 1. OR 2. MUST BE SELECTED*,//, 

S IX, 'OPTION SELECTED ?») 

READd IN, *, ERR* 1 20 , END* 9999 ) W ( 642 ) 

C OUTPUT UNITS OPTION 

140 WRITE(IOUT,150) 

150 F0RMAT(//,1 OX, 'OUTPUT UNITS OPTIONS',//, 

S 1X,»0. ENGLISH',/, 

$ 1X,*1. METRIC',//, 

$ IX, 'OPT ION SELECTED ?•) 

READ ( 1 1 N, *, ERR* 1 40 , END* 9999 ) W ( 648 ) 

C ■ STREAM I NE OPTION — 

160 WRITEOOUT, 170)NC+1 ,NC 

170 FORMAT (//, IX,' IS CASE *,I3,» ALONG THE ', 

$ 'SAME STREAM INE* ,/, IX, 'AS CASE »,I3,' ?') 

READd IN,65,ERR-160,END*9999)ANS 
IF(ANS.EQ.1HY)GQ TO 200 


C-9 


180 

190 


’”““**“771 " PROGRAM INPUT CONTROL PARAMETERS 

WRITE (I OUT, 190) 

FORMAT (//,1 OX,' INPUT CONTROL FLAG',// 

5 IX, '1. NEW CASE FOLLOWS USING TITLE, TIMING AND TRAJ • 

S / ,6X, *DATA FROM PREVIOUS CASE 1 ,/, * * 

{ IX, ‘2. END OF INPUT (LAST CASE)*,/, 

i !?icp\ / N f! £? U0WS USINS TRA1.DATA FROM PREVIOUS • 
$ * UNCHANGED? T TLE AND T,MING * ,N,T,AL CASE DATA *, 

$ IX, *4, NEW CASE FOLLOWS USING NEW TITLE. TIMING TRAl t 
J ’AND CASE DATA *,/,6X, '(INITIALLY ZERO W ARRAY )S/? AJ * * 

$ wmous wev!" EXCEPT ZER0 au CASE DATA FR0M •« 

$ ilw^CASEWAVA EXCEFT INIT,A1 - IZE ZER0 TIMING 
S IX, ’OPTION SELECTED ?*) 

READ( I IN,* 1 , ERR* 1 80, END S 9999 ) W ( 64 1 ) 

60 TO 300 


200 CONTINUE 
W ( 641 )*1 .0 
CALL STREAM 
300 CONTINUE 
RETURN 

9999 CONTINUE 
CALL EXIT 
END 


STREAM I NE CASES 




n o 


SUBROUTINE CROSS 


ROUTINE TO CHOOSE CROSS FLOW OPTIONS 
COMMON/ WARRAY/W ( 700 ) 

CQMMQN/UNIT/I IN, I OUT 
COMMON/MI SC/NC, JFKS,TFLAG,MFLAG, ANSI 

INTEGER ANS, TFL AG, ANSI 
C MAY USE CROSS FLOW DATA FROM PREVIOUS CASE IF JFKS - 1 OR 3 
IFUFKS.NE.1 .AND. JFKS.NE.3)G0 TO 9 
NCM1 -NC-1 

1 WR I TE ( I OUT, 2 ) NC, NCM1 

2 FORMAT <///,' IS CASE ',I3,» CROSS FLOW DATA THE SAME AS 
$ 'FOR CASE ',13,' ?•) 

READ( 1 I N, 160, ERR* 1 ,END*9999)ANS 
I F (ANS. EQ. 1 HY )RETURN 
9 OONTINUE 
10 WRITE(IOUT,20) 

20 FORMAT(/// , IX, 1 0X, * CROSS FLOW ADJUSTMENTS*,//, 

$ 10X,*1. CONSTANT WIDTH RECTANGLE ( IDEAL GAS)',/, 

S 10X,'2. CONSTANT WIDTH RECTANGLE (REAL GAS)',/, 

$ 10X, '3. SHARP EDGE DELTA WING ( IDEAL GAS)',/, 

$ 10X,'4, SHARP EDGE DELTA WING(REAL GAS)',//, 

$ IX, 'OPTION SELECTED ?•) 

READ( I IN, #,ERR*1 0, END® 9999) W( 201 ) 

IT=IFIX(W(20D) 

GO TO (100,200,300,400) IT 
GO TO 500 

C CROSS FLOW OPTION 1 

100 IF(MFLAG.EQ.O)WRITE( I0UT,110) 

I F (MFLAG. EQ. 1 ) WR I TE ( I OUT, 111) 

1 1 0 F0RMAT(//,1X, ' 1 . CONSTANT WIDTH RECTANGLE ( I DEAL GAS) ' , 
$ //, IX, 'RECTANGLE WIDTH(FT) ?') 

111 FORMAT (// ,1X,'1 . CONSTANT WIDTH RECTANGLE (I DEAL GAS)', 
$ //, IX, 'RECTANGLE WIDTH (M) ?«) 

READ( I IN, *, ERR®1 00,END®9999)W(202) 

120 IF(MFLAG.EQ.O)WRITE( IOUT,130) 

IF (MFLAG. EQ. 1 )WRITE( IOUT, 131 ) 

130 FORMAT (IX, 'EDGE RADIUS(FT) ?') 

131 FORMATdX, 'EDGE RADIUS(M) ?•) 

READ( I IN, *, ERR® 1 20, END® 9999) W( 205) 

140 WRITE(I0UT,1 50) 

150 FORMAT (IX, 'ANY CHANGES ?») 

READ ( I IN, 160, ERR® 1 40 , END-9999) ANS 
160 FORMAT(AI) 

IF(ANS.EQ.1HY)G0 TO 100 
GO TO 600 

C CROSS FLOW OPTION 2 

200 IF(MFLAG.EQ.0)WRITE(I0UT,210) 

I F (MFL AG . EQ. 1 ) WR I TE ( I OUT, 21 1 ) 

210 F0RMAT(//,1X, »2. CONSTANT WIDTH RECTANGLE (RE*! GAS)'/ 
i //, IX, 'RECTANGLE WIDTH(FT) ?') 

211 F0RMAT(//,1X,'2. CONSTANT WIDTH RECTANGLE (REAL GAS) ', 

$ //, IX, 'RECTANGLE WIDTH(M) ?') 

READ( 1 1 N, *, ERR-200, END* 9999) W( 202) 


220 

230 

240 

250 


C — 
300 
310 


320 


C — 
400 
410 


420 

430 


440 


C ~ 
500 

510 

600 

610 

620 


9999 


IF(MFLAG.EQ.Q)WRITE(|0UT,130) 
IF(MFLAG.EQ.1 )WRITE( IOUT, 131 ) 

READ(| IN, *, ERR*22Q, END=9999) W (205) 
WRITE (I OUT, 240) 

FORMAT (IX, ‘REAL GAS VELOCITY GRADIENT 
$ 'SURFACE, 1.0 SWEPT CYLINDER) ?') 
READ< I IN,*, ERR*230, END*9999) W( 204) 
WRITE( I OUT, 1 50) 

READ( | IN, 160,ERR*230,END*9999)ANS 
IF(ANS.E0.1HY)GO TO 200 
GO TO 600 


(0.31 FLAT ', 


CR0SS FLCW 0PT,0N 3 — 

FORMAT(//,1X, *3, SHARP EDGE DELTA WING( IDEAL GAS) 
$ //, IX, 'DELTA WING SWEEP ANGLE (DEG) ?») 

READ ( 1 1 N, ", ERR*? 00, END* 9999) W ( 203 ) 

WRITE( I0UT,1 50) 

READ (I I N, 1 60, ERR*320, END=9999) ANS 
IF(ANS.EQ.1HY)G0 TO 300 
GO TO 600 


woss nm opt,on 4 

FORMAT(//,iX,M. SHARP EDGE DELTA NINGIREAL GAS)' 
$ //, IX, 'DELTA WING SWEEP ANGLE (DEG) ?») 

READ( UN,*, ERR*400, END*9999) W( 203 ) 

WRITE(I0UT,430) 

FORMATdX, 'REAL GAS VELOCITY GRAD I ENT (SHARP EDGE ', 
S 'DELTA WING * 0.31, POINTED CONE ■ 1.0) ?') 

READ ( 1 1 N, *, ERR*420 , END* 9999 ) W ( 204 ) 

WRITE (I OUT, 150) 

READ ( 1 1 N, 1 60, ERR*440, END=9999) ANS 
IF(ANS.EQ.1HY)G0 TO 400 
GO TO 600 


CONTINUE 

WRITE(I0UT,510) 

FORMAT (//, IX, 'BAD OPTION',/) 

GO TO 10 
CONTINUE 
WRITE (I OUT, 620) 

F0RMAT(/,1X, ' IS THE OPTION OORRECT ?') 

READ ( 1 1 N, 1 60, ERR*6 1 0, END*9999) ANS 

IF(ANS.EQ.1HN)G0 TO 10 

RETURN 

CONTINUE 

CALL EXIT 


c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


c 


SUBROUTINE FLOW (ID) 

ROUTINE TO SET FLOWFIELD AND LOCAL PRESSURE OPTIONS 
ID - 1 CREATING NEW FLOWFIELD DATA 
ID - 2 MODIFYING OLD FLOWFIELD DATA 
COMMON/ WARR AY/ W ( 700 ) 

COMMON/UNIT/I IN, I OUT 

COMMQN/POQEF F/TMCP ( 50 ) ,TCPM( 50 ) , NCPMT 

COMMON/M I SC/NC, JFKS, TFL AG, MFLAG, ANSI 

FF(N) - FLOWFIELD OPTION 

FFA(N) - ANGLE CORRESPONDING TO FF (N) 

P(N) - LOCAL PRESSURE OPTION 

PA(N) - ANGLE CORRESPONDING TO P(N> 

FNAM5 - NAME FOR PRESSURE COEFFICIENT INPUT TABLE FILE 
DIMENSION FF(9),FFA(9),P(9),PA<9) 

INTEGER ANS, TFL AG, ANSI 
CHARACTERS 0 CHARF(4),CHARP<6) 

CHARACTER *20 FNAM5 


DATA CHARF/ ' SHRP-WEDGE ' , 

S ‘SHRP-CONE •, 

S »OBL IQ/NORM', 

$ ‘PARALLEL »/ 

DATA CHARP/‘CP-VS-MACH‘, 

$ ‘TAN- WEDGE ', 

$ ‘TAN-CONE ', 

S ‘OBLIQ-SURF*, 

$ ‘MOD-NEWT. », 

$ ‘PRANDT-MEY*/ 

C MAY USE FLOWFIELD DATA FROM PREVIOUS CASE IF JFKS = 1 OR 3 
IF(JFKS.NE.1 .AND. JFKS.NE.3)G0 TO 9 
NCM1 -NC-1 

1 WRITE(I0UT,2)NC,NCM1 

2 F0RMAT(1H1,‘ IS CASE *,13, • FLOWFIELD DATA THE SAME AS », 
$ ‘FOR CASE ',13* ?•) 

READd IN, 20, ERRS ,END=9999)ANS 
I F ( ANS. EQ. 1 HY JRETURN 

9 CONTINUE 
W(261)-1.0 
W(262)-2.0 
W(263)*0.0 
W(264)-0.0 
W(265)*0.0 
W(266)-0.0 
W(267)*0.0 
W(268)*0.0 
W(269)«0.0 
DO 5 1-1,9 
FF(I>— 1.0 
PCD— 1.0 
FFA( I )-0.0 
PA( I )-0.0 
5 CONTINUE 
C 


C- 13 


10 

15 

20 

c — 
c 

100 


110 


120 

130 

200 

220 

230 

300 

C — 
C 

500 

505 

510 

515 

520 

525 

600 

610 


WRITE ( IOUT, 1 5) 

FORMAT <///» 1 0X, * FLOWF IFFD AND FOCAL PRf S51JRF OPTION*!' //) 
FORMAT (A1) * 

IFdD.FQ.I )00 TO 1000 

CONVFRT W ARRAY TO FLOWF If LD/PRF SCURF lAOLL ~ 

CONTINUE 
1*1 
J«31 
K k 37 

CONTINUE 

IF(W( J) ,EQ. 0.0. OR, I.GT.9)GQ TO 500 
lf(W(J).LT.30.)G0 TO 200 
IF(W(J).GT.37.)G0 TO 120 
FF( I )°W(J)-34. 

FFA( I )=W(K) 

60 TO 130 
FF( I )=W( J)-35. 

FFAdJ^WOO 
J=J+1 
K=K-> 1 

IF( J.EQ.37) J=46 
I F(FF ( I ).EQ. 4. )GO TO 300 
IF(W( JJ.EQ.29. )G0 TO 220 
P(l )=W(J)-13. 

PA( I )=W(K) 

GO TO 230 
P( I )=6. 

PA( I )=W(K) 

J=J+1 

K=K+1 

IF(J.EQ.37)J=46 

1 = 1+1 

GO TO 110 


ECHO FLOWF I ELD/PRESSURE DATA 
CONTINUE 
1=1 


WRITE (I OUT, 505) 

F0RMAT(///,1X,25X, 'TOTAL EFFECTIVE ANGLE ») 
WRITE(I0UT,51 5) I 
FORMAT (IX, ' SET' ,11) 

i^i ,KGT ‘°* 0)WRITE(,0UT » 520)CHA RE(FF(l)),(FFA(J) 
F0RMAT(2X,'FF',6X,A10,5X,» ALPHA + »,9(F7.3s,' + »)) 
1F(P(|).GT.0.0)WRITE(I0UT,525)CHARP(P(|)),(PA(J),J« 
F0RMAT(2X, 'P* ,7X,A1 0,5X, 'ALPHA + ',9(F7.s!,' ♦ • ) 
1 = 1+1 ' 


> 1 , 1 ) 

, 1 ) 


GO^TO^IO^*"* • >«EQ.— 1 ,0)GO TO 600 

WRITE(IOUT,610) 

FORMAT (//, IX, ' ANY CHANGES ?») 

RE AD ( 1 1 N, 20, ERR*6 00, END=9999) ANS 
IF(ANS,EQ.1HY)G0 TO 1000 
GO TO 2000 


o o 


C INPUT FLOWFIELD/PRFSSURF DATA 
1000 CONTINUE 
1010 WRITE ( (OUT, 1 01 5) 

1015 FORMAT <///,1X,6(* -'),' FLOWFIELD* ,6( • -'),3X,6<» -'), 
t ' PRESSURE', 6{* IX, »-1. FLOWFIELD TYPE NOT NEEDED', 

S 6X,'-1. PRESSURE TYPE NOT NEEDED* ,/,2X, • 1 . SHARP WEDGE ', 
l 'SHOCK ANGLE', 1 IX, *1. INPUT CP VS MACH NO. TABLE', /,2X, 

S *2. SHARP CONE SHOCK ANGLE»,12X,'2. TANGENT WEDGE PRESSURE', 
$ /,2X,'3. OBLIQUE AND NORMAL SHOCK (90 DEG)', IX, 

$ '3. TANGENT COPE PRESSURE', /,2X, '4. PARALLEL SHOCK .PRES ' 

$ 'NOT NEEDED) ',2X,»4. OBLIQUE SURFACE PRESSURE' ,/,39X, 

S *5. MODIFIED NEWTONIAN', /,39X, '6. PRANDTL-MEYER EXP. ', 

S '(FF NOT NEEDED) ',//,! X,«? ( • *»),» INPUT OPTIONS IN PAIRS ', 

$ 'WITH ASSOCIATED DELTA ANGLES',4(* *'),/, IX, • (NOTEi TO ♦, 

S 'SIGNIFY END OF CASE USE -1.0 FOR BOTH FF AND PRESS. », 

$ 'OPTIONS.)',/, IX, 'FLOWFIELD, DEL ANGLE, PRESSURE, DEL ANGLE') 

DO 1020 1-1,9 
1030 WRITE( I OUT, 1035)1 
1035 FORMATOX, II ) 

READ( I IN, ", ERR-1 030,END-9999)FF( I ),FFA( I ),P( I ),PA( I ) 
IF(FF(I).EQ.-1..AND.P(I).EQ.-1.)G0 TO 5U0 
1020 CONTINUE 


OONVERT FLOWL I ELD/PRESSURE TABLE TO W ARRAY 

2000 CONTINUE 

DO 2001 1*1,9 

IF(FF(I ).GT.2.0)FF(I )«FF(l)+35.0 
I F (FF < I ).GE.1 .0.AND.FFd ) »LE.2.0)FF( I )®FF(I )+34.0 
IF(P(I ).EQ.6.0)P(I )*29.0 

2001 IF(P( I ) .GE. 1 .0.AND.P( I ).LE.5.0)P( I )=P( I )+13.0 
1-1 

J-31 

JP1-J+1 

K-37 

KP1-K+1 

2100 IF(FF(I).EQ.-1..AND.P(I).EQ.-1.)G0 TO 3000 
IF(P(I).EQ.-1.)G0 TO 2300 
IF(FF(I ).EQ.-1 .)G0 TO 2200 
W(J)-FF(I) 

W(K)-FFAd) 

W(JPD-Pd) 

W(KPD-PAd) 

J-J+2 
K-K+2 
2150 1-1+1 
JP1-J+1 
KP1-K+1 

IF( J.EQ.37) J*46 
IF( J.EQ.38) J-47 
IF( JP1 .EQ.37) J-46 
I F ( JP1 .EQ.38) J a 47 
GO TO 2100 
2200 W(J)-Pd) 

W(K)-PAd) 

J-J+1 
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K*K+1 

GO TO 2130 
2300 W(JWFC » 

W(K)*FFA( I ) 

J*J+1 
K*K+1 

GO TO 2150 



C CHECK FOR PRESSURE TYPE 

3000 CONTINUE 
DO 3001 L*1 ,9 

3001 IF(P(L).EQ.14.0)GO TO 3100 
GO TO 5000 

3100 CONTINUE 

WRITE( I OUT, 3110) 

3110 $ F »T^LE* 1 ) 1J0X#,,NPUT PRESSUR£ DEFICIENT VS. MACH NO. 
W(649)*1.0 

3120 WRITE(I0UT,3125) 

3125 F0RMAT(1X,»IS THERE AN INPUT FILE AVAILABLE ?») 

READ( I IN, 20, ERR*3 1 20, END*9999) ANS 
IF(ANS.EQ.1HN)G0 TO 4000 

C INPUT MACH TABLE VIA FILE 
3130 WRITEC I0UT,3135)NC 

3135 FORMAT (IX, 1 WHAT IS THE FILE NAME FOR CASE * . 13 • ?t) 
READ ( 1 1 N, 3 1 40, ERR*3 130, END® 9999 ) FNAM5 # 

3140 FORMAT (A20) 

OPEN( UNI T*8, F ILE*FNAM5 , STATUS* »OLD» , ERR*8888) 

READ (8,31 60 )NCPMT 
3160 F0RMAT(I3) 

WRITE(I0UT,3165)NCPMT 
31 ° 5 FORMAT (//, IX, 13, » MACH NUM3ERS' ) 

WRITE(I0UT,3166) 

3166 F0FMAT(//,5X, 'M-INF 1 ,T21, •CP’,/) 

DO 3333 l*1,NCPMT 
REA.M8,3170)TMCP( I ),TCPM( I ) 

3170 FORMAT(2F10.6) 

WRITE(,0UT,3175)TMCP(I ),TCPM(I ) 

3175 FORMAT(1X,2F10.6) 

3333 CONTINUE 

CL0SE(UNIT*8,ST/YTUS* •KEEP' ) 

3150 WRITE(I0UT,3185) 

3185 F0RMAT(/,1X, • IS THIS THE DATA YOU WANT ?») 

READ( 1 1 N, 20, ERR*3 1 80, END*9999) ANS 
iF(ANS.EQ. 1HY)G0 TO 5000 
CLOSE (UN I T*8, STATUS* ’KEEP 1 ) 

GO TO 3120 


1 


INPUT MACH TABLE VIA TERMINAL 
4000 WRITE (IOUT,4010) 

4010 FORMATOH1 ,** » « INPUT VIA TERMINAL * * «».//. IX 

NG.S WILL BE INPUT ? (MAXIMUM OF *0)*) 
READ( 1 1 N, «, ERR*4000, END=9999) NCPMT 
WRITE < I0UT,401 5) 
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4015 FORMAT ( IX, 'M- INF, CP' ) 

DO 4444 1*1, NCPMT 
4020 HR I TE ( I OUT, <>025 ) I 
4025 FORMATdX, !2) 

READ( I IN, *,ERR JB 4020,END=9999)TMCP( I ),TCPM( I ) 

4444 CONTINUE 
4030 WR I TE ( I0UT,4035) 

4035 FORMAT (/, 1 X, ' ARE THERE ANY CHANGES TO THE INPUT ?») 
READ( I IN,20,ERR=4030,FND=9999)AN5 
IF(ANS.EQ.1HY)G0 TO 4000 

5000 CONTINUE /» EXIT SUBROUTINE 

RETURN 

8880 WRITE (I OUT, 8889) 

8889 FORMAT (/, 1 X, ' UNABLE TO OPEN FILE.') 

GO TO 3120 
9999 CONTINUE 
CALL EXIT 
END 
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SUBROUTINE HEAT IN 
C 

C ROUTINE FOR SETTING W ARRAY VALUES FOR A HEATING 
C INDICATOR CASE 

COMMON/ WARRAY/W (700) 

COMMON/UN IT/ 1 IN, I OUT 
COMMON/T I TL E/T I TL 1 
C 

CHARACTER*72 TITL1 
INTEGER ANS 
C 


10 WRITE( IOUT,20) 

20 F0RMAT(///,1X,» HEATING INDICATOR*,//) 

30 FORMATtAI ) 


WRITEdOUT, 
40 FORMATdX, * 
$ IX,' 

$ IX, » 

$ IX, » 

$ IX,' 

W<12)*1 .0 
W<11)*1 .0 
W(24)*0.0 
W(31)*38.0 
W(37)*90.0 
W(32)»18.0 
W( 38) *90,0 
W(642)*2,0 
W( 644) *0.0 
W(646)*0.0 
W(315)*1 .0 
RETURN 


40) 


FAY AND RIDDEU',/, 

RADIUS * 1 FT SPHERE',/, 

WALL TEMP * 0 DEG. F',/, 

LEWIS NO. « 1.0',/, 

SUMMARY PRINT ONLY',/) 

/* RADIUS 

/* FAY _RIDDELL 

/* WALL TEMPERATURE (F) 

/* OBLIQUE SHOCK 
/« SHOCK ANGLE 

l* MODIFIED NEWTONIAN PRESSURE 
/* BODY ANGLE 
/* SUMMARY PRINT ONLY 


/* LEWIS NO. * 1.0 


9999 CONTINUE 
CALL EXIT 
END 
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SUBROUTINE HTMULT 


a 


ROUTINE TO SET HEAT TRANSFER MULT I PI I CAT I ON FACTORS 
CQMMQN/WARRAY/W(70Q) 

COMMON/UNIT/I IN, IOUT 
COMMON/M I SC/NC, JFKS, TFL AG, MFL AG, ANSI 
C 

INTEGER ANS,TFLAG,AN$1 

C MAY USE MULTIPLIER DATA FROM PREVIOUS CASE IF JFKS « 1 OR 3 
IFUFKS.NE.1. AND. JFKS. NE.3)GO TO 9 
NCM1-NC-1 

1 WRITE(IOUT,2)NC, NCM1 

2 FORMAT (///, 1 X, * I S CASE • , 13,* MULTIPLICATION FACTOR DATA ', 
$ 'THE SAME AS FOR CASE ',13, » ?') 

READd IN,30,ERR*1 ,END=9999)ANS 
I F (ANS.EQ. 1 HY )RETURN 
9 CONTINUE 
10 WRITE(I0UT,20) 

20 FORMAT (///, 1 X, 1 0X, * HEAT TRANSFER MULTIPLICATION FACTORS*, 

$ //,1X,*0PTI0N TYPES' ,/,10X, 

$ *1. CONSTANT VALUE*, /,1 OX, 

$ *2. MULTIPLIER A FUNCTION OF TIME»,/,10X, 

S *3. MULTIPLIER A FUNCTION OF FREESTREAM MACH NO.', 

$ //,1X,»N0TE: MULTIPLIERS ARE MULTIPLIEDd.E. », 

$ 'AMPLIFICATION * TYPE 1 *TYPE2*TYPE3 )»,//) 

30 FORMAT (A1 ) 

C MULTIPLIER TYPE 1 

100 WRITE (I OUT, 110) 

110 FORMAT (//, 1 X, ' 1 . CONSTANT VALUE', //,1X, 

S 'DO YOU WANT THIS TYPE ?') 

READd IN,30,ERR*100,END*9999)ANS 
IF(ANS.EQ.1HN)G0 TO 200 
120 WRITE (IOUT, 130) 

130 FORMAT( IX, 'ENTER LAMI NAR,TURBULENT VALUES') 

READd IN, »,ERR»120,END*9999)W(18),W(19) 

140 WRITE(IOUT,150) 

150 FORMAT ( IX, » ANY CHANGES ?') 

READd IN, 30, ERR- 140, END* 9999) ANS 
IF (ANS.EQ. 1HY) GO TO 100 

C MULTIPLIER TYPE 2 

200 WRITE( IOUT,210) 

210 F0RMAT(//,1X,'2* MULTIPLIER A FUNCTION OF TIME', 

$ //,1X, 'DO YOU WANT THIS TYPE ?») 

READ( 1 1 N, 30, ERR-200, END-9999) ANS 
IF(ANS.EQ.1HN)G0 TO 300 
220 WRITE( IOUT,230) 

230 FORMAT ( IX, * NUMBER OF TIMES (10 MAXIMUM) ?») 

READ( UN,*, ERR*220 , END® 9999 ) W ( 320 ) 

240 WRITE(I0UT,250) 

250 F0RMAT(1X,'TIME,LAM MULTIPL IER,TURB MULTIPLIER') 
IT«W(320)+.0001 
DO 1000 1-1, IT 
260 WRITE(I0UT,265) I 
265 F0RMAT(I2) 

READd IN, *,ERR=260,END»9999)W(320+I ),W(330+I ),W(340+I ) 
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1000 CONTINUE 
270 WR I TE ( I OUT, 1 50 ) 

READ ( II N, 30, ERR=27 0, END* 9999) ANS 
IF!ANS.EQ.1HY)G0 TO 200 

C MULTIPLIER TYPE 3 

300 WRITE! I0UT,305) 

305 FORMAT!//, IX, '3. MULTIPLIER A FUNCTION OF FREESTREAM', 
$ • MACH NO. »,/,5X,» LINEAR INTERPOLATION IN LOG(M-INF)', 

5 • VS LOG(MULTIPLIER) SPACES//, IX, *D0 YOU WANT THIS’, 

$ ' TYPE ?') 

READ ( 1 1 N, 30, ERR=300, END® 9999) ANS 
IF(ANS.EQ«1HN)G0 TO 400 
310 WRITE! IOUT,315) 

315 FORMAT! 1X,33HNUIBER OF MACH »S (10 MAXIMUM) ?) 

READ ( 1 1 N, *, ERR=3 1 0, END*9999) W(360) 

WRITE! I OUT, 320) 

320 FORMAT! IX, »M- INF, LAM MULT,TURB MULT') 

ITT«W(360)+.0001 
DO 2000 l®1 , ITT 
330 WRITE! I OUT, 265) I 

READ! I IN, *,ERR®330,END=9999)W(360+I ),W(370H ),W(380+I ) 

C TEST FOR VALUE LESS THAN OR EQUAL TO 0.0 
IF(W(360+I).LE.0.0)G0 TO 333 
IF(W(370+I ).LE. 0.0)60 TO 333 
IF(W(380+I ).LE. 0.0)60 TO 333 
GO TO 2000 
333 WRITE! I OUT, 335) 

335 FORMAT! IX, 'VALUE CANNOT BE LESS THAN OR EQUAL TO 0.0',/) 
WRITE! IOUT,320) 

60 TO 330 
2000 CONTINUE 
340 WRITE! I OUT, 150) 

READ! I IN,30,ERR=340,END®9999)ANS 
IFtANS.EQ. 1HY)60 TO 300 
C CHECK NUNBER OF TABLE ENTRIES (MIN®2) 

IF(ITT.LT.2.0)GO TO 500 
DO 3000 1=1 , ITT 
C CONVERT TO LOG VALUE 

W(360+l )=AL0G1 0(W(360+I )) 

W(370+I)=AL0610(W(370+I)) 

W(380+l )=AL061 0! W( 380+ 1 ) ) 

3000 CONTINUE 
400 CONTINUE 
RETURN 

500 WRITE! I0UT,510) 

510 FORMAT! IX, 'MIN I MUM NUI^ER MACH NO. - 2') 

60 TO 310 
9999 CONTINUE 
CALL EXIT 
END 
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SUBROUTINE HTRMTD 


ROUTINE FOR SELECTING HEAT TRANSFER METHOD 

XLEWNO - LEWIS NUteER 
C0MM0N/WARRAY/W(700) 

COMMON/UN IT/ 1 IN, I OUT 
COMMON/MI SC/NC, JFKS, TFL AG, MFLAG, ANSI 

CHARACTER*50 CHAR(II) 

INTEGER ANS, TFL AG, ANSI 
C TITLES FOR HEAT TRANSFER OPTIONS 

OATA CHAR/ 'HEMISPHERE STAGNATION POINT*, 

$ »CATO/ JOHNSON SWEPT CYLINDER*, 

$ ‘ECKERT REF. ENTHALPY FLAT PLATE METHOO', 

S * ECKERT/SPAULD I NG-CH I FLAT PLATE METHOD', 

S 'BOEING RHO-MU FLAT PLATE METHOD*, 

S 'BECKWITH/GALLAGHER SWEPT CYLINDER METHOO*, 

$ 'BOEING RHO-MU SWEPT CYLINDER METHOD*, 

* 'LEES/DETRA-HIDALGO HEMISPHERE DISTRIBUTION*, 

$ » LEES IDE ORB ITER HEATING', 

$ »FLAP REATTACHMENT HEATING', 

$ 'FIN-PLATE PEAK INTERFERENCE HEATING*/ 

C MAY USE DATA FROM PREVIOUS CASE IF JFKS * 1 OR 3 
IF( JFKS.NE.l .AND. JFKS. NE.3)G0 TO 9 
NCM1-NC-1 

1 WRITE! I0UT,2)NC, NCM1 

2 FORMAT! 1H1,» IS CASE »,I3,» HEAT TRANSFER DATA SAME AS », 

$ »FOR CASE *, 13,* ?») 

READ! 1 1 N, 70, ERR* 1 ,END*9999)ANS 
I F ! ANS. EQ. I HY )RETURN 
9 CONTINUE 

10 WRITE! IOUT,20)(I,CHAR( I ), 1*1,11) 

20 FORMAT!///,* HEAT TRANSFER METHOD*,//, IX, ’OPTIONS*, 

$ 3X, 12,*. *,A50,/,1 0! 1 IX, 12,' . *,A50,/>) 

C CHOOSE A HEAT TRANSFER METHOD 

30 WRITE!I0UT,40) 

40 FORMAT! IX, 'OPTION SELECTED ?•) 

READ! JIN,*, ERR-30, END-9999) W( 1 1 ) 

50 WRITE! I0UT,60) 

60 FORMAT! IX, 'SHOULD RAREFIED FLOW HEATING BE INCLUDED ?') 

READ! I IN,70,ERR*50,END*9999)ANS 
70 FORMAT !A1 ) 

W( 646) -0.0 

IF!ANS.EQ.1HY)W(646)-1 .0 
80 WRITE! IOUT,90) 

90 FORMAT ! IX, ' I S THE HEAT TRANSFER OPTION CORRECT ?•> 

READ! 1 1 N,70,ERR*80, END-9999) ANS 
IF(ANS.EQ.1HN)G0 TO 10 

GO TO <100,200,300, 400,500,600,700,800, 900,1000,1 100), W(11) 
GO TO 1200 

c heat transfer METHOD 1 

100 I F (MFLAG. EQ. 0) WRITE! I OUT, 1 10)CHAR( 1 ) 

IF(MFLAG*EQ»1 ) WR I TE < I OUT, 1 11 )CHAR(1) 

110 FORMAT! 1H1,»1. »,A50,//,1X, 'BODY RADIUS (FT) ?•) 
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Ill FORMATdHI, ‘1. *,A50,//,1X, 'BODY RADIUS (M) ?*) 

READ ( 1 1 N, *, ERR*1 00, END*9999) W( 1 2) 

120 WRITE! (OUT, 130) 

130 FORMAT! IX, • LEW IS NO. * 1 .0 OR 1 .4 ?•) 

READ (I IN, *, ERR-1 20, END*9999)XLEWN0 
IF!XLEWN0.EQ.1 .0)W(315)«1 .0 
IF(XLEWN0.EQ.1 ,4)W(315)»0.0 
I F ! W< 646 ) .EQ.1 .0) WRITE! I OUT, 140) 

160 WrTteT I OUT^I 70?"^ * op, T I ON HAS BEEN SELECTED.*) 

170 FORMAT! IX, • ANY CHANGES ?') 

READ! I IN,70,ERR=161 JD*9999)AN$ 

IF!ANS.EQ.1HY)60 TO 100 
GO TO 1300 


200 

210 

211 

220 

230 

240 


260 


300 

310 

311 

320 

330 

340 

350 

360 

370 


375 

376 


heat transfer method 2 — 

IF!MFLAG.EQ.O)WRITE! I0UT,21 0) CHAR (2) 

IFtMFLAG.EQ.I )WRITE! I0UT,21 1 )CHAR!2) 

FORMAT! 1 HI , *2. »,A50,//,1X, ‘BODY RADIUS (FT) ?») 

FORMATdHI ,*2. *,A50,//,1X, ‘BODY RADIUS (M) ?») 

READ < I IN,*, ERR=200, END*9999) W( 1 2) 

WRITE! IOUT,230) 

FORMAT ( IX, * SWEEP ANGLE ?») 

READ! I IN, *,ERR S 220,END=9999)W(17) 

WRITE! (OUT, 130) 

READ! I IN, *,ERR=240,END=9999)XLEWN0 
IF(XLEWN0.EQ.1 .0)W(315)*1 .0 
IF(XLEWN0.EQ.1 ,4)W(315)*0.0 
WRITE! I OUT, 170) 

READ! I IN,70,ERR*260,END*9999)ANS 
IF(ANS.EQ.1HY)G0 TO 200 
GO TO 1300 

HEAT TRANSFER METHOD 3 

IF(MFLAG.EQ.O)WRITE( I0UT,310)CHAR(3) 

IFtMFLAG.EQ.I ) WR I TE ( I OUT, 311 )CHAR(3) 

FORMATdHI, *3. *,A50,//, IX, ‘RUNNING LENGTH (FT) ?‘) 

FORMATdHI, *3. *,A50,//,1X, ‘ RUNNING LENGTH (M) ?») 

READ! I IN,*,ERR=300,END=9999)W( 13) 

WRITE! I OUT, 330) 

FORMAT (IX, ‘TURBULENT MANGLER FACTOR ?») 

READ! I IN,*,ERR«320,END=9999)W(16) 

WRITE! I0UT,350) 

FORMAT !1X,‘ LAMINAR MANGLER FACTOR ?») 

READ ( I IN,*, ERR“340, END=9999) W ( 1 5) 

WRITE! IOUT,370) 

FORMATdX, ‘SURFACE DISTANCE TO START OF TURBULENT B.L.*. 

$ /# 1 X, * TH I S DISTANCE IS SUBTRACTED FROM THE RUNNING *, 

$ ‘LENGTH*,/, IX, ‘FOR TURBULENT HEATING CALCULATIONS. *,/, IX, 

$ ‘DESIRED LENGTH ?•) * * 

READ! I IN,*,ERR«360,END S 9999)W(354) 

WRITE! I OUT, 376) 

FORMATdX, 1 IS AN AUTOMATIC VIRTUAL ORIGIN CORRECTION •, 

S ’DESIRED ?*) 


READ! 1 1 N, 70, ERR=37 5, END® 9999) ANS 
IF(ANS.EQ. 1HN)W(29)*0.0 
IF(ANS.EQ. 1HY)W(29)®1 .0 


C-22 


360 WRITE* I0UT, 170) 

READ* I IN,70, ERR-360 , END- 9999 ) AN S 
IF*ANS.EQ.1HY)GQ TO 300 
60 TO 1300 

C HEAT TRANSFER METHOD 4 

400 IF*MFLAG.EQ.0)WRIT£(IOUT,405)CHAR*4) 

IF*MFLAG.EQ*1 )WRITE(I0UT,406)CHAR(4) 

405 FORMAT* 1H1,»4. *,A50,//, IX, 'RUNNING LENGTH *FT) 7') 

406 FORMAT* 1H1, *4. »,A50,//,1X, 'RUNNING LENGTH <M> ?») 

READ* II N, *, ERR-400, END-9999 ) W * 1 3 ) 

410 WRITE* I OUT, 370) 

READ* I IN,*,ERR-41 0,END-9999)W<354) 

415 WRITE* I0UT,416) 

416 FORMAT* IX,' IS AN AUTOMATIC VIRTUAL ORIGIN CORRECTION ', 
$ 'DESIRED 7') 

READ* 1 1 N,70,ERR*41 5, END-9999) ANS 
IF(ANS«EQ*1HN)W(29)*0.0 
IF(ANS.EQ.1HY)W*29)-1 .0 
420 WRITE* I0UT,330) 

READ* 1 1 N, «, ERR-420, END-9999) W( 1 6 ) 

430 WRITE* I0UT,350) 

READ ( 1 1 N, », ERR-43 0, END-9999 ) W * 1 5 ) 

440 WRITE(I0UT,445) 

445 FORMAT*/, 1 X, • REYNOLDS-ANALOGY FACTOR' ,//,1 IX, 

$ »0. COLBURN' ,/,1 IX, »1. VON KARMAN»,//,1X, 

$ 'DESIRED FACTOR ?•) 

READ* I IN,*,ERR-440,END-9999)W(319) 

IF*W*646).NE.1.0)GO TO 470 
450 WRITE* I0UT,455) 

455 FORMAT*/, IX, 'RAREFIED FLOW OPTIONS' ,/,1 IX, 

$•0. C0NE»,/,11X,'1 • FLAT PL ATE',/, IX, 

$ 'DESIRED OPTION ?') 

READ* I IN,*,ERR»450,END“9999)W(314) 

IF(W*314).NE.0.0)G0 TO 470 
460 WRITE* I0UT,465) 

465 FORMAT*/, IX,' RAREFIED CONE OPTIONS' ,/,1 IX, 

$•0. SHARP CONE' ,/,1 IX, *1. BLUNT OONE',/, IX, 

$ 'DESIRED OPTION ?•) 

READ* II N, *, ERR-460, END-9999) W* 650 ) 

470 WRITE* IOUT,170) 

READ* 1 1 N,70,ERR»470, END-9999) ANS 
IF*ANS.EQ.1HY)G0 TO 400 
GO TO 1300 

C HEAT TRANSFER MTHOD 5 

500 IF*MFLA6.EQ.0)WRITE* I0UT,510)CHAR(5) 

IF(MFLAG*EQ*1 ) WRITE* I OUT, 51 1 )CHAR(5) 

510 FORMAT* 1H1,»5. »,A50,//, IX, 'RUNNING LENGTH *FT) ?') 

511 F0RMAT*1H1 ,'5. *,A50,//, IX, 'RUNNING LENGTH *M) ?') 

READ* 1 1 N, *, ERR-500, END-9999) W* 1 3 ) 

520 WRITE* IOUT,370) 

READ * II N, •, ERR 520, END-9999) W( 354) 

525 WRITE* I OUT, 526) 

526 FORMAT* IX,' IS AN AUTOMATIC VIRTUAL ORIGIN CORRECTION 
$ 'DESIRED 7') 

READ* 1 1 N, 70, ERR-525, END- 9999) ANS 
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530 


600 

610 

611 


620 


630 

640 


700 

710 

711 

720 

730 


800 

810 

811 

820 

825 

826 

830 

835 


840 


IF(ANS,EQ.1HN)W(29)*-0.o 
IF(ANS.EQ. 1HY)W(29)*1 .0 
WRITE( IOUT, 170) 

READ( I IN,70,ERR»530,END-9999)ANS 
IF(ANS.EQ.1HY)GO TO 500 
GO TO 1300 

heat transfer ncthod 6 

IF(MFLAG.EQ.0)WRITE(IOUT,610)CHAR(6) 

IF(MFLr\G.EQ.1 ) WR I TE ( I OUT, 611 )CHAR(6) 

™MAT(1H1,’6. »,A50,//,1X, ’RADIUS (FT) ?’) 

J2SJT. <1H1 * ,# * *»A50,//, IX, ’RADIUS (M) ?•) 

READ( I IN,*, ERR*600,END B 9999)W( 12) 

WRITE (I OUT, 620) 

IF(W(646).EQ.1 ,0)WRITE(IOUT,630) 

$ F »SELECTED X ‘) N0TEl THE RAREF,ED CYl,n DER OPTION HAS BEEN 
WRITE (I OUT, 170) 

READ( I IN,70,ERR*640,END a 9999)ANS 
IF(ANS.EQ. 1HY)G0 TO 600 
GO TO 1300 

HEAT TRANSFER ACTHOD 7 

IF(MFLAG.EQ.O)WRITE( I0UT,710)CHAR(7) 

1F(MFLAG.EQ.1 )WRITE( I OUT, 711 )CHAR(7) 

F0RMAT(1H1,*7. ’,A50,//, IX, ’RADIUS (FT) ?’) 

0RMAT(1H1,’7. »,A50,//, IX, ’RADIUS (M) ?») 

READ ( 1 1 N, *, ERR®700, END®9999) W( 1 2) 

WRITE (I OUT, 620) 

WRITE (I OUT, 170) 

READ ( I | N, 70, ERR=730, END® 9999) ANS 
IF(ANS.EQ.1HY)G0 TO 700 
GO TO 1300 

heat transfer method 8 

IF(MFLAG.EQ.0)WRITE(IOUT,810)CHAR(8) 

IF(MFLAG.EQ.1 ) WR I TE ( t OUT, 8 1 1 )CHAR(8) 

F0RMAT(1H1 , »8. »,A50,//, IX, ’RADIUS (FT) ?») 

RFAn?MM H i'cDD aA ^ A50 »//#1 x # , ^ RADIUS (M) ?») 

READ( 1 1 N, *, ERR-800, END-9999) W( 1 2) 

IF(MFLAG.EQ.0)WRITE(I0UT,825) 

I F ( MFL AG . EQ. 1 ) WR I TE ( I OUT, 826 ) 

FORMAT( IX, ’RUNNING LENGTH (FT) ?’) 

FORMAT( IX, ’RUNNING LENGTH (M) ?») 

READ ( 1 1 N, *, ERR=820, END“9999) W( 1 3 ) 

WRITE( I0UT,835) 

FORMAT( IX, ’LOCAL BODY SLOPE (DEG.) ?•) 

READ(| IN, *, ERR»830,END=9999)W( 17) 

WRITE (I OUT, 840) 

FORMAT (/ , IX, ’NOTE: THIS OPTION IS USED WITH’./. 

t 7V ter ncrrir\Ai * ,M,n * 


850 


-»w . _ w wjtu ni in,/, 

< 0P I!°? /* AND NORMAL SHOCK, SHOCK 

5 .SJSf “ 90 » P OPTION 5. MODIFIED ’, 

S ’NEWTONIAN, ANGLE ■ BODY SLOPE •) 

WR|TE(IOUT,170) 
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READ( I IN,70,ERR-850,END-9999)ANS 
IF(ANS.EQ.1HY)G0 TO 800 
60 TO 1300 

c heat TRANSFER METHOD 9 

900 IF(MFLA6.EQ.0)WRITE(I0UT,910)CHAR(9) 

IF(MFLAG.EQ.1 >WRITE( I0UT,91 1 )CHAR(9) 

910 F0RMAT(1H1,»9. »,A50,//,1X, 'FULL SCALE VEHICLE », 

$ 'USE RAOIUS • 1.0 FT. ',/, IX, 'RADIUS 'HD ?•> 

911 F0RMAT(1H1,»9. *,A50,//,1X,'FULL SCt,.. '. VEHICLE ', 

$ 'USE RADIUS ■ .3048 M. »,/, IX, 'RADIUS (M) ?') 

READ( 1 1 N, *, ERR- 900, END-9999) W< 1 2) 

920 !F(MFLAG.EQ.O)WRITE(IOUT,925) 

I F (MFL AG. EQ. 1 ) WR I TE ( I OUT, 926 ) 

925 FORMAT( IX, 'WINDWARD WALL ENTHALPY',/, IX, 

$ '(480 BTU/LEM - 2000 R ASSUMED IF ZERO IS INPUT) ?') 

926 FORMAT( IX, 'WINDWARD WALL ENTHALPY',/, IX, 

$ 'C.12E6 JOULES/KG » 1111 K ASSUMED IF ZERO IS INPUT) ?») 
READ ( 1 1 N, ", ERR-920, END-9999) W ( 21 ) 

WRITE(I0UT,930) 

930 FORMATt/, 1 X, ' NOTE t THIS OPTION IS USED WITH',/, 

$ 7X,»FF OPTION 3. OBLIQUE AND NORMAL SHOCK, SHOCK ', 

S 'ANGLE » 90.«,/,7X,» P OPTION 5. MODIFIED NEWTONIAN,', 

$ ' ANGLE • 90.') 

940 WRITE(I0UT,170) 

READ( 1 1 N, 7 0, ERR- 940, END-9999) ANS 
IF(ANS.EQ.1HY)G0 TO 900 
GO TO 1300 

c HEAT TRANSFER NETHOD 10 

1000 IF (MFLAG.EQ. 0 ) WR I TE ( I OUT, 1 005) CHAR (10) 

IF(MFLAG.EQ.1 )WRITE( (OUT, 1006)CHAR(10) 

1005 F0RMAT(1H1 ,»10. »,A50,//,1X, 'RUNNING LENGTH ', 

$ 'TO HINGE LINE (FT) ?') 

1006 F0RMAT(1H1 ,'10. ',A50,//, IX, 'RUNNING LENGTH *, 

$ 'TO HINGE LINE (M) ?') 

READ ( 1 1 N, *, ERR- 1000, END-9999 ) W ( 1 3 ) 

1010 WRITE(I0UT,370) 

READ( II N, *, ERR-1 01 0, END-9999) W( 354) 

1015 WRITE(IOUT,1016) 

1016 FORMAT ( IX, ' I S M AUTOMATIC VIRTUAL ORIGIN CORRECTION ', 

$ 'DESIRED 7') 

READ( I IN,70,ERR»1015,END»9999)ANS 
I F (ANS. EQ. 1 HN) W( 29) -0 .0 
IF(ANS.EQ.1HY)W(29)«1 .0 
1020 WRITE(IOUT,1025) 

1025 FORMATdX, 'TURBULENT MANGLER FACTOR ?») 

READ( I IN, *, ERR-1 020, END-9999) W( 16) 

1030 WRITE (IOUT, 1035) 

1035 FORMATdX, 'LAMINAR MANGLER FACTOR ?') 

READ( II N, », ERR- 1 030, END-9999) W( 1 5) 

1040 IF(MFLAG.EQ.O)WRITE( IOUT, 1045) 

I F (MFL AG . EQ. 1 ) WR I TE ( I OUT, 1 046 ) 

1045 FORMATdX, 'FLAP LENGTH (FT) ?») 

1046 FORMATdX, 'FLAP LENGTH (M) ?•) 

READ( 1 1 N, ", ERR- 1 040, END-9999) W( 22) 

1050 WRITE( IOUT, 1055) 


1055 FORMAT ( IX, • REYNOLDS ANALOGY FACTOR* f A OX. 

$ *0. COLBURN* ,/ ,1 OX, * 1 . VON K ARMAN*,/, 

$ IX, ‘DESIRED FACTOR ?•) 

RE AD ( II N, «, ERR- 1050, END- 9999 ) W ( 3 1 9) 

WRITE( IOUT, 1 060) 

1060 FORMAT( IX, * NOTEs FLAP ANGLE IS INPUT THROUGH *. 

J ‘LAST PRESS. OPTION ANGLE.*) 

1070 WRITE (IOUT, 170) 

READ( 1 1 N, 7 0, ERR- 1 070, END-9999) ANS 
IF(ANS.EQ.1HY)G0 TO 1000 
GO TO 1300 

C HEAT TRANSFER METHOD 11 

1100 IF(MFLAG.EQ.0)WR)TE(I0UT,1105)CHAR(11) 

I F (MFL AG. EQ. 1 ) WR I TE ( I OUT, 1 1 06 )CHAR (11) 

1 1 05 FORMAT ( 1H1 , • 1 1 . *,A50,//, IX, 'RUNNING LENGTH », 

$ 'TO FIN LEADING EDGE (FT) ?') 

1 106 FORMAT ( 1H1 , ' 1 1 . »,A50,//, IX, 'RUNNING LENGTH *, 

% *T0 FIN LEADING EDGE (M) ?») 

READ ( 1 1 N, *, ERR- 1100, END-9999) W (13) 

1110 WRITE( IOUT.370) 

READ ( 1 1 N, **, ERR- 1110, END-9999) W( 354) 

1115 WR I TE ( IOUT, 1 1 16) 

1116 F0RMAT(1X,» IS AN AUTOMATIC VIRTUAL ORIGIN CORRECTION ». 

$ 'DESIRED ?') 

READ ( 1 1 N, 70, ERR- 1115, END-9999) ANS 
IF(ANS.EQ.1HN)W(29)=0.0 
IF(ANS.EQ.1HY)W(29)-1 .0 
1120 IF (MFLAG.EQ.O ) WR ITE ( I OUT, 1125) 

IF(MFLAG.EQ. 1 )WRITE( IOUT, 1 1 26) 

1125^F0RMAT( IX, » D| STANCE ALONG FIN TO POINT OF INTEREST », 

1126 FORMAT( IX, 'DISTANCE ALONG FIN TO POINT OF INTEREST ', 

I (M) ?*) 

READ( I IN, »,ERR=1 1 20, END-9999) W( 25) 

1130 WRITE(I0UT,1135) 

1135 F0RMAT(1X,'FIN ANGLE AT ALPHA - BETA * 0 (DEG) ?') 

READ( I IN,*,ERR®1 1 30, END-9999) W( 26) 

1140 WRITE( IOUT, 1 145) 

1145 FORMATdX, 'FIN ANGLE 0PTI0N',/,10X, 

$ *0. FIN ANGLE - INPUT FIN ANGLE',/, 1 0X, 

$ '1. FIN ANGLE - INPUT FIN ANGLE + ALPHA»,/,10X, 

S »2. FIN ANGLE - INPUT FIN AI^LE ♦ BETA',//,1X, 

$ ‘OPTION DESIRED ?•) "" # 

READd IN,*,ERR-1 140,END«9999)W(30) 

1150 WRITE( IOUT,330) 

RE AD ( 1 1 N, ", ERR- 1 1 50 , END- 9999 ) W ( 1 6 ) 

1160 WRITE( I OUT, 350) 

READ( I IN, *,ERR-1 1 60,END=9999) W( 1 5) 

1170 WRITE (I OUT, 445) 

READd IN. *,ERR-1170,END-9999)W(319) 

1180 WRITE( IOUT, 170) 

READd IN,70,£RR-1 180, END=9999)ANS 
IF(ANS.EQ.1HY)G0 TO 1100 
GO TO 1300 

C B AD OPTION — 
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1200 CONTINUE 

WRITE(I0UT,1210) 

1210 FORMAT!//, IX, 'BAD OPTION*,//) 

GO TO 10 

C SET WALL TEMPERATURE 

1300 CONTINUE 

1310 IF(MFLAG.EQ.0)WRITE(I0UT,1320) 

IF(MFLAG.EQ.1 ) WRITE! I0UT,1321 ) 

1320 FORMAT!///, IX,** * * HALL CONDITIONS * • *»,//, IX, 

$ »HALL TEMPERATURE (DEG F) ?•) 

1321 FORMAT!///, IX,'* • * HALL CONDITIONS * * *»,//, IX, 

S «HALL TEMPERATURE (DEG K) ?') 

READ( 1 1 N, », ERR- 1310, END- 9999) H( 24) 

C SET HALL EMISSIVITY — 

1330 HRITE(I0UT,1340) 

1340 FORMATdX, *HALL EMISSIVITY ?•) 

READ(IIN,*,ERR-1330,END-9999)W(23) 

C CONTINUATION OPTION 

1400 IF(MFLAG.EQ.O)WRITE( IOUT, 1 41 0) 

I F ( MFL AG . EQ . 1 ) HR I TE ( I OUT, 1411) 

1410 F0RMAT(///,1X,»* * * CONTINUATION OPTION * * *»,//,1X, 

$ *D0 YOU HANT TO PROVIDE AN INITIAL HEATING LOAD GT *, 

$ *0.0 (BTU/SQ.FT) ?•) 

1411 F0RMAT(///,1X, ** * * CONTINUATION OPTION * * «',//, IX, 

$ *D0 YOU HANT TO PROVIDE AN INITIAL HEATING LOAD GT *, 

$ *0.0 (JOULED/SQ.M) ?•) 

READ( 1 1 N,70,ERR»1 400, END-9999) ANS 
H( 209) -0.0 

IF(ANS.EQ.1HY)H(209)»1 .0 
IF(ANS.EQ. 1HN)G0 TO 1300 
1420 IF(MFLAG.EQ.O)HRITE( I OUT, 1430) 

I F(MFLAG.EQ. 1 )HRITE( IOUT, 1 431 ) 

1430 FORMATdX,* INITIAL LOAD (BTU/SFT) ?•) 

1431 FORMATdX,* INITIAL LOAD (JOULED/SQ.M) ?*) 
READ(IIN,*,ERR=1420,END»9999)H(316) 

C TEST RAREFT.D FLOW FLAG 

1500 CONTINUE 

IF(H(646).NE.1.0)G0 TO 1600 

GO TO (1600,1510,1510,1600,1510,1600,1510,1600,1510,1510 
$ ,1510), Hdl) 

1510 HR I TE( I OUT, 1520) 

1520 FORMAT (/, IX,* THE RAREFIED FLOH OPTION CANNOT BE USED HITH », 
$ *HEAT TRANSFER* ,/,1X, ‘OPTION SELECTED.*,//, IX, 

$ * * RAREFIED FLOW OPTION IS DISENABLED «»,//,1X, 

$ »D0 YOU WISH TO CHANGE HEAT TRANSFER OPTIONS ?») 

W(646)-0.0 

READ( I IN,70,ERR»1 51 0, END-9999) ANS 
IF(ANS.EQ.1HY)G0 TO 10 
1600 CONTINUE 
RETURN 

9999 CONTINUE 
CALL EXIT 
END 
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SUBROUTINE MODIFY 


C 

C THIS ROUTINE ENABLES USER TO PICK A SECTION OR ROUTINE 
C TO BE MODIFIED OR REDEFINED 
COMMON/ WARR AY/ W (700) 

COMMON/UNIT/I IN, I OUT 

COMMON/MI SC/NC, JFKS,TFLAG,MFLAG, ANSI 

C 

INTEGER ANS,TFLAG, ANSI 
C 

17 WRITE(I0UT,16) 

10 FORMAT (///,1 OX, ’MODIFICATION SECTIONS’,//, 

$ IX,' 1. TIMING PARAMETERS',/, 

$ IX, » 2. TRAJECTORY DATA',/, 

$ IX,' 3. ATMOSPHERE DATA',/, 

$ IX, » 4. FLOWFIELD AND PRESSURE DATA',/, 

$ IX, » 5. CROSSFLOW DATA',/, 

$ 1X,» 6. TRANSITION CRITERIA',/, 

$ IX, » 7. HEAT TRANSFER OPTION’,/, 

$ IX,’ 8. HEATING MULTIPL IERS* ,/, 

$ IX,' 9. GEOMETRY DATA',/, 

$ IX, '10. CONTROL PARAMETERS',/, 

$ IX, *11. HEATING INDICATOR',/, . 

$ IX, '12. OR CHANGE A SPECIFIC VARIABLE IN W ARRAY’,// 

$ IX, 'SECTION TO BE MODIFIED ?') 

READd IN,«,ERR=17,END=9999)SECT 

GO TO (1,2,3, 4, 5, 6,7, 8,9, 10,11, 12), SECT 

GO TO 17 

1 CALL TIMING 
GO TO 30 

2 CALL TRAJ 
GO TO 30 

3 CALL ATMS 
IF(W(10).EQ.1.)CALL WNDTUN 
IF(W(10).EQ.2.)CALL ATMDTA 
GO TO 30 

4 CALL FL0W(2) 

GO TO 30 

5 CALL CROSS 
GO TO 30 

6 CALL TRANS 
GO TO 30 

7 CALL HTRMTD 
GO TO 30 

8 CALL HTMULT 
GO TO 30 

9 CALL TDGEOM 
GO TO 30 

10 CALL CONTRL 
GO TO 30 

11 CALL HEAT IN 
GO TO 30 

12 GO TO 200 

30 WRITE( I0UT,40)NC . 

40 F0RMAT(/,1X, 'DO YOU WISH TO MODIFY ANY OTHER SECTIONS , 
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$ 'FOR CASE ',13,' ?•) 

READ< 1 1 N, 50, ERR-30, END«9999) ANS 
50 FORMAT (At) 

IF(ANS.EQ.1HY)G0 TO 17 
RETURN 

C MODIFY A SPECIFIC W NUMBER 

200 CONTINUE 
210 WRITE (I OUT, 220) 

220 FORMAT </, IX, ‘WHAT IS THE W NUMBER YOU WISH TO CHANGE ?•) 
REAO( II N, », ERR-21 0, END-9999) J 
230 WRITE(I0UT,240)J 

240 FORMAT ( IX, ' WHAT IS THE VALUE FOR W(»,I3,') ?•) 

READ< 1 1 N, *, ERR* 230, ENO-9999) W( J ) 

250 WRITEdOUT, 260) 

260 FORMATdX, 'IS THERE ANOTHER W VARIABLE YOU WISH TO ', 

$ 'CHANGE ?•) 

READ( I IN, 50, ERR-250, END-9999) ANS 
IF(ANS.EQ.1HY)G0 TO 200 
GO TO 30 
9999 CONTINUE 
CALL EXIT 


oo 


SUBROUTINE STREAM 


ORIGINAL PAGE is 
OF POOR QUALITY 


c 

C ROUTINE FOR CREATING STREAML INE CASES 
C 

c X - X(FT) OR X(M) 

C SA - SHOCK ANGLE (OEG) 

C BA - BODY ANGLE (OEG) 

C BP - BODY POINT NUMBER 
C CASE - CASE NUMBER 

C XPRT - PROGRAM INPUT CONTROL PARAMETER FOR CASE 
C J1 - W INDEX FOR CASE 
C J2 - W index FOR X 
C J3 - W INDEX FOR SA 
C J4 - W INDEX FOR BA 
C J5 - W INDEX FOR BP 
C J6 - W INDEX FOR XPRT 
C JK - FLAG FOR SPECIAL CASES OF STREAMLINE 
C ~ (PARALLEL SHOCK OR PRANDT-MEYER OPTIONS) 
COMMON/ WARRAY/W( 700 ) 

COMMON/UN IT/ 1 IN, I OUT 
COMMON/ STRM/ I STRM, NNC 
„ COMMON/M I SC/NC, JFKS, TFL AG, MFLAG, ANSI 

INTEGER TFLAG, ANSI 


CHECK FOR TIME DEPENDENT GEOMETRY OPTION 
IF(W(560).GT.0.0)G0 TO 200 

C CRE ^„f„^r? R ^. r /! L£ F0R ST0RING STREAM INE CASE OATA 
OPEN( UN I T*4,F I L E* • STREAM. TMP* , STATUS* 1 UNKNOWN' ) 
NWONC+1 
ISTRM*1 
J 1*647 
J2=13 

DETERMINE SHOCK ANGLE AT PREVIOUS CASE - 

18,46,-1 

IF(W( I ).GT.34.AND.W( I ).LT.40)G0 TO 27 
DO 26 1*36,31,-1 
1 1 1 * 1+6 

IF(W( I ).GT.34.AND.W( I ),LT.40)G0 TO 27 


25 

DO 25 1* 

111*1-3 

IF(W(I). 

26 

DO 26 i*: 
1 1 1*1+6 
IF(W( 1 ).( 

27 

J3*l 1 1 


J3F*I 


30 


40 

50 


DO 30 |*48 E 4 ERM,NF B0DY AN6LE AT PREV,0LIS CASE 

11*1-3 

IF(W( I ).GT.13.AND.W( I ).LT.30)G0 TO 50 
DO 40 1*36,31,-1 
1 1 * 1+6 

IF(W(I).GT.13.AND.W(| ).LT.30)G0 TO 50 

J4* 1 1 

J4F*I 

JK*0 

CHECK FOR PARALLEL SHOCK OR PRANDT-MEYER CASES 

IF(W( J4FKEQ.29) JK*1 
•F( J3F.GT. J4F) JK*2 

C-30 


f 

k 




t 







IF( JK.EQ.1 ) J4“0 ORIGINAL PAGE IS 

IF( JK.EQ.2) J3-0 OF POOR QUALITY 

J5-611 

J6-641 

XPRT-1 .0 

HRITE(IOUTJ) 

IF(MFLAG.EQ.O.AND.JK.EQ.O)WRITE(IOUT,2> 

IF(MFLAG.EQ.O. AND. JK.EQ.1 )WRITE< I0UT,3) 

IF(MFLA6.EQ.0.AND. JK.EQ.2)WRITE( I0UT,4) 

IF(MFLA6.EQ.1 .AND. JK.EQ.O)WRITE< I0UT,5) 

I FtMFLAG.EQ.I. AND. JK.EQ.1 )WRITE( I0UT,6) 

IF(MFLAG.EQ.1 .AND. JK.EQ.2)WRITE( I0UT,7) 

1 FORMAT ( 1H1 , 1 0X, 1 STREAM I NE CASES’ ,//*1 IX, •LIMITATIONS* j^*'^**^ 

$ »2X, ’ALL LOCATIONS ON STREAM. I NE HAVE NO INTERVENING SHOCKS OR’, 
$ • EXPANSION FANS’, /,5X,’*»,2X,»MJLTIPLIERS SA1C FOR ALL X », 

$ ’LOCATIONS’ ,/,5X,»*»,2X, ’FIXED RUNNING LENGTHS WITH ALPHA’,/, 

$ 5X,’**,2X,*N0TE: SHOCK ANGLE • SHOCK ANGLE OR SHOCK GENERATOR ’, 
$ ’ANGLE’ ,//,1X,»( INPUT DATA FOR EACH CASE - - USE NEG X ’, 

$ ’TO INDICATE END OF DATA)’) 

2 FORMAT(/, IX, 'X(FT), SHOCK ANGLE, BODY ANGLE,B.P.NO. ') 

3 F0RMAT(/,1X,’X(FT), SHOCK ANGLE,B.P.NO.’> 

4 F0RMAT1 /, 1 X, ’ X ( FT ) ,BODY ANGLE, B.P. NO.’ ) 

5 FORMATt/, IX, ’X(M), SHOCK ANGLE, BODY ANGLE, B.P. NO.* ) 

6 F0RMAT(/,1X,’X(M), SHOCK ANGLE, B.P. NO. • ) 

7 F0RMAT(/,1X,*X(M),B0DY ANGLE,B.P.NO.») 

60 WRITE(I0UT,70)NNC 

70 FORMATdX, 13) 

^(MFLAG^EQ^AND. JK.EQ.O)READ( I IN,*,ERR=60,END=9999)X,SA,BA,BP 
I F (MFL AG. EQ.O . AND. JK. EQ. 1 ) READ( 1 1 N, «, ERR»60,END»99' , 9)X, SA, BP 
IF(MFLAG.EQ.O.AND. JK.EQ.2)READ( I IN,*,ERR=60,END»9999)X,BA,BP 
IF(MFLAG.EQ.1 .AND. JK.EQ.O)READ( I IN,*,ERR s 60,END s 9999)X,SA,BA,BP 
I F(MFLAG.EQ. 1 .AND. JK. EQ. 1 )READ( I IN,*, ERR=60,END=9999)X, SA,BP 
IF(MFLAG.EQ.1 .AND. JK.EQ.2)READ( I IN,*,ERR-60,END-9999)X,BA,BP 
IF(MFLAG.EQ.1 )X-X*3.28084 

WRITE TO TEMPORARY FILE 


WRITE(4,80)J1,CASE, J2,X, J3,SA, J4,BA, J5,BP 
WRITE(4,80) J6,XPRT 
80 F0RMAT(1X,5(I3,F9.4,1X)) 

IF(X.LT.0.0)G0 TO 100 
NNONNC+1 
GO TO 60 
100 CONTINUE 

REWIND (UNIT-4) 

RETURN 

200 CONTINUE 

WRIVE(I0UT,210) 

210 FORMAT(/, IX, 'CANNOT RUN STREAM. I NE CASES AND USE », 
) »TI1C DEPENDENT GE01CTRY • ' ) 


RETURN 

9999 CONTINUE 
CALL EXIT 
END 
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SUBROUTINE TDGEOM 

ROUTINE FOR SETTING TIME DEPENDENT GEOMETRY DATA 
COMMQN/WARRAY/WC700) 

COMMON/UN IT/ 1 IN, I OUT 
COMMON/MI SC/NC, JFKS, TFL AG, MFLAG, ANSI 
C 

INTEGER ANS, TFL AG, ANSI 

C MAY USE GEOfrCTRY DATA FROM PREVIOUS CASE IF JFKS « 1 OR 3 
IF< JFKS. NE.1. AND. JFKS. NE.3)GO TO 9 
NCM1 *NC-1 

1 WR I TE ( I OUT, 2 ) NC, NCM1 

2 FORMAT (///, IX,* IS CASE ',I3,» TIME DEPENDENT GEOMETRY », 

$ ‘DATA THE SAME AS FOR CASE *,13, • ?•) 

READd IN,30,ERR*1,END*9999)ANS 
I F (ANS. EQ. 1 HY )RETURN 
9 CONTINUE 
10 WRITEdOUT,20) 

20 FORMAT (///, 1 0X, * TIME DEPENDENT GEOMETRY*,//) 

30 FORMAT (A1 ) 

C CHOOSE NUNBER OF TIME DEPENDENT ENTRIES — 

40 WRITE(I0UT,50) 

50 FORMAT ( IX, » NJMBER OF TIME DEPENDENT ENTRIES ? (MAX*10)») 

READ (I IN,#,ERR*40,END*9999) I I 
W( 560) “FLOAT d I ) 

60 IF(MFLAG.EQ.O)WRITE( I0UT,70) 

I F (MFL AG. EQ. 1 ) WR I TE ( I OUT, 7 1 ) 

70 FORMATdX, *TIIC(SEC), RADI I (FT), LENGTH (FT), SLOPE OR SWEEP*, 
$ *(DEG) ») 

71 FORMATdX, »TII€(SEC), RADI KM), LENGTH(M), SLOPE OR SWEEP*, 

$ »(DEG) ») 

c input values for each time 

DO 1000 1=1,11 

READC I IN,*, ERR=60, END=9999) W( 560+ 1 ) , W( 570+ 1 ) , W( 580+ 1 ) , 

$ W(590+l ) 

1000 CONTINUE 
80 WRITE(I0UT,90) 

90 FORMATdX, 'ANY CHANGES ?») 

READd IN,30,ERR»80 END-9999)ANS 
IF(ANS.EQ.1HY)G0 TO 10 
100 CONTINUE 
RETURN 

9999 CONTINUE 
CALL EXIT 
END 
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SUBROUTINE TIMING 


ROUTINE THAT SETS TIMING PARAMETERS AND PRINT CONTROL VALUES 

SUM • TOTAL NUMBER OF PRINT TIMES 
COMMON/ WARRAY/W ( 700 ) 

COMMON/UN IT/ 1 IN, I OUT 

INTEGER ANS 

SUM1-0.0 
SUM2-0.0 
SUM3-0.0 
WRITE(I0UT,20) 

FORMAT ( 1 HI, 5X,» SPECIFY PRINTOUT INTERVALS*,/) 
WRITE(I0UT,40) ' 

FORMAT( IX, » INITIAL TIME (SEC)*) 

READd IN, »,ERR*30,END*9999)W(1) 

WRITE (I OUT, 60) 

FORMAT < IX, * PR INTOUT INTERVAL 1 (SEC) »,1 OX, ‘DELTA TIME*) 
READd IN,«,ERR*50,END*9999)W(2) 

WRITE(I0UT,80) 

FORMAT ( IX, ‘SECOND TIME (SEC)*) 

READd IN,*,ERR*70,END*9999)W(3) 
IF(W(2).GT.0.0)SUM1*((W(3)-W(1))/W(2)) 
write(iout,ioo) 

FORMAT (IX, » PR INTOUT INTERVAL 2 (SEC) *,1 0X,» DELTA TIME*) 
READd IN, «,ERR*90,END*9999)W(4) 

WRITEdOUT, 1 20) 

FORMATdX, ‘THIRD TIME (SEC)*) 

READd IN, *,ERR*110,END*9999)W(5) 

IF(W(4).GT.0,0)SUM2*( (W(5)-W(3) )/W(4) ) 

WRITEdOUT, 140) 

FORMATdX, ‘PRINTOUT INTERVAL 3 (SEC) *,1 OX,* DELTA TIME*) 
READd IN, »,ERR*130,END«9999)W(6) 

WRITEdOUT, 160) 

FORMATdX, FOURTH TIME (SEC)*) 

READ (UN.*, ERR* 150, END* 9999) W ( 7 ) 
IF(W(6).GT.0.0)SUM3»((W(7)-W(5))/W(6)) 

W(8)*1 .0 
WRITEdOUT, 180) 

F0RMAT(/,1X, ‘ARE THE PRINTOUT TIMES CORRECT ?') 

READd IN,190,ERR*170,END*9999)ANS 
FORMAT (A1 ) 

IF(ANS.EQ.1HN)G0 TO 10 
SUM*SUM1 +SUM2+SUM3 

C TEST FOR GREATER THAN 100 PRINT TIMES 
IF(SUM, LT. 1O0.)GO TO 1000 
WRITEdOUT, 200) 

200 FORMAT(/, IX, 'MAXIMUM NUIBER OF TIMES(IOO) EXCEEDED.*,/) 

GO TO 30 
1000 CONTINUE 
RETURN 

9999 CONTINUE 
CALL EXIT 


10 

20 

30 

40 

50 

60 

70 

80 


90 
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160 
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180 

190 


END 
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SUBROUTINE TRAJ 
C 

C ROUTINE THAT COPIES TRAJ DATA TO W ARRAY 
C 

C LINE - TRAJ TABLE LINE NUMBER 
C FNAM1 - NA« OF TRAJ FILE 
C TITLE - TITLE OF TRAJ FILE 
C BFLAG - 0 YAW ANGLE NOT INCLUDED 
C BFLAG » 1 YAW ANGLE INCLUDED 

COMMON/WARRAY/W(700) 

COMMON/ UN IT/ 1 IN, I OUT 
COMMON/M I SC/NC, JFKS, TFL AG, MFL AG, ANSI 
C 

CHARACTERS FNAM1 
CHARACTERS TITLE 
INTEGER ANS, BFLAG, TFL AG, ANSI 
C 


ORIGINAI 

POOF 


10 WRITE( I0UT,20) 

20 K ORMAT (1 HI, 1 0X, ‘TRAJECTORY INPUT 1 ,/) 

30 WRITE( I0UT,40) 

40 FORMATdX, ‘DO YOU HAVE A TRAJECTORY INPUT FILE ?») 

READd IN,50,ERR=30,END*9999)ANS 
50 FORMAT (A1) 

(F(ANS.EQ.IHN)GO TO 130 

C INPUT TRAJECTORY VIA FILE 
60 WRITE(I0UT,70) 

70 FORMAT (IX, ‘WHAT IS THE FILE NAIC ?') 

READd IN,80,ERR*60,END*9999)FNAM1 
80 FORMAT (A20) 

OPEN ( UN I T*7 , F I L E-FN AMI , STATUS* ‘OLD*, ERR* 8888) 
READ(7,85)TITLE 
85 FORMAT (A80) 

C BFLAG * 1 IF YAW ANGLE INCLUDED 
C W(50) - NUMBER OF TIME DEPENDENT TABLE ENTRIES 
READ(7,90)BFLAG,W(50) 

90 FORMAT(I2,F10.5) 

NFTS*IFIX(W(50)) 

C READ TIME, ALT, VEL, ANGLE OF ATTACK, AND YAW ANGLE (OPTIONAL) 

IF(BFLAG.EQ.0)READ(7,99,END*98) (W( 50+ I ),W( 100+1 ),W( 150+1 ), 

$ W(210+l ), 1-1 ,NPTS) 

IF(BFLAG.EQ.1 )READ( 7, 100, END-98) (W(50+l ),W(100+I ),W(150+I ), 
$ W(210+l ),W(650+I ), 1*1 ,NPTS) 

98 CONTINUE 

99 FORMATUE15.4) 

100 FORMAT (5E1 5.4) 

IF(MFLAG.EQ.O)GO TO 103 
DO 101 12*101 ,200 

101 W( IZ)*W( IZ)/3. 28084 

103 CONTINUE 

CLOSE ( UNIT-7 , STATUS* ‘KEEP’ ) 

C 

C DISPLAY TRAJECTORY 

104 CONTINUE 


1*0 


PAGE IS 

QUALITY 
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— \r 


105 


110 


112 


111 


113 


115 


120 


123 

125 


CONTINUE 


IF(BFLAG.EQ.O.AND. 
I F (BFLAG. EQ.O. AND. 
IF (BFLAG. EQ.1 .AND. 
IF (BFL AG. EQ.1 .AND. 
FORMATS, IX,' 

$ 'ATTACK' ,/,' 

S ' (DEG)*,/) 

F0RMAT(/,1X,» 

$ ’ATTACK* ,/,« 

$ » (DEG)*,/) 

F0RMAT(/,1X,» 

$ 'ATTACK BETA*, 
$ • (DEG) 

F0RMAT(/,1X, ' 

$ 'ATTACK BETA*, 
S ' (DEG) 

l»l+1 


MFLAG.EQ.0)WRITE(I0UT,1 
MFL AG. EQ.1 )WRITE( I0UT,1 
MFLAG.EQ.O)WRITE( I OUT, 1 
MFLAG.EQ.1 )WRITE< IOUT, 1 
TI1C ALTITUDE 
(SEC) (FT) 


TIME 

(SEC) 


ALTITUDE 

(M) 


TIME ALTITUDE 

/,' (SEC) 

(DEG)*,/) 

Tlte ALTITUDE 

/,' (SEC) 

(DEG)*,/) 


10 ) 

12 ) 

11 ) 

13) 

VELOCITY ANGLE ', 
(FT/SEC)*, 

VELOCITY ANGLE », 
(M/SEC)*, 

VELOCITY ANGLE ', 
(FT) (FT/SEC)*, 

VELOCITY ANGLE ', 
(M) (M/SEC)', 


IF(BFL AG.EQ.0)WRITE( IOUT, 120) I, W(50+I),W( 100+1), 
$ W ( 1 50+ I ),W(210+I ) 


I F (BFLAG. EQ. 1 ) WRITE< IOUT, 1 20) I , W< 50+ 1 ),W(100+I ), 

S W(150+I),W(210+I),W(650+I) 

FORMATOX, I2,5E12.4) 

IF(I.EQ.NPTS)GO TO 210 

IF( I .EQ.20.0R. I .EQ.40)READ( I IN,1 25, ERR-1 23, END-9999) 

IFU.EQ.20.0R. I.EQ.40)GO TO 105 

FORMAT(IX) 

GO TO 115 


INPUT TRAJECTORY VIA TERMINAL 
130 CONTINUE 
140 WRITE( IOUT, 150) 

150 FORMATOX, 'WHAT IS THE NIMJER OF TRAJECTORY POINTS ? », 

S '(50 TRAJ.PTS. MAXIMUM)') 

READ( 1 1 N, *, ERR-1 40, END-9999) NPTS 
W ( 50 ) -FLOAT ( NPTS ) 

160 WRITE( I0UT,170) 

170 FORMATOX, 'WILL BETA VALUES BE INPUT ?') 

READO IN,50,ERR-160,END-9999)ANS 
BFLAG-0 

I F ( ANS. EQ. 1 HY )BFL AG-1 

I F (BFLAG. EQ. 0 . AND. MFL AG. EQ. 0 ) WR I TE ( I OUT, 1 80 ) 

I F (BFL AG. EQ. 0 . AND . MFL AG . EQ . 1 ) WR I TE ( I OUT, 1 82 ) 

I F (BFL AG. EQ. 1 . AND. MFL AG. EQ. 0 ) WR I TE ( I OUT, 1 81 ) 

I F (BFLAG. EQ. 1 . AND. MFLAG. EQ. 1 ) WR I TE ( I OUT, 1 83 ) 

180 FORMATOX, 'TYPE IN THE FOLLOWING TRAJECTORY VARIABLES ', 
$ »SEPERATED BY COMMAS*, /,1 1 ( » -»),» 50 TIMES MAXIMUM', 

$ 11(* -'),/,' Tl>€( SEC), ALTITUDE (FT), VELOCITY (FT/SEC),', 
$ 'ANGLE OF ATTACK(DEG) •) 

182 FORMATOX, 'TYPE IN THE FOLLOWING TRAJECTORY VARIABLES *, 
$ 'SEPERATED BY COMMAS', /,1 1 ( » -'),» 50 TIMES MAXIMUM', 

S IK' -'),/,' T I >€ (SEC), ALT ITUDE(M), VELOCITY (M/SEC),', 

S 'ANGLE OF ATTACK(DEG) ') 

181 FORMATOX, 'TYPE IN THE FOLLOWING TRAJECTORY VARI/fcLES ', 
) 'SEPERATED BY COMMAS', /,1 1 ( • -•),» 50 TIMES MAXIMUM', 
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$ IK* -*),/,* T I ME (SEC), ALTITUDE (FT), VELOCITY (FT/SEC), ', 

$ * ANGLE OF ATTACK(DEG),BETA(DEG) ') 

183 FORMAT (IX,* TYPE IN THE FOLLOWING TRAJECTORY VARIABLES ', 

$ ‘SEPERATED BY COMMAS' ,/,1 1 ( • -'),' 50 TIMES MAXIMUM*, 

S IK* -»),/,' TIME(SEC),ALTITUDE(M), VELOCITY (M/SEC),*, fVV V 1*5 

$ * ANGLE OF ATTACK(DEG),BETA(DEG) •) OF POOR G*. : AUYY 

DO 1000 1-1 ,NPTS 
190 WRITE( IOUT,200) I 
200 F0RMAT(1X, 12) 

I F (BFL AG. EQ. 0 ) READ ( 1 1 N, «, ERR- 1 90, END- 9999) W( 50+ 1 ) , 

$ W(100+I),W(150+I),W(210+I) 

I F(BFLAG. EQ. 1 )READ( 1 1 N, *, ERR-1 90,END«9999) W( 50+ 1 ), 

$ W( 100+1 ),W( 150+1 ),W(210+I),W(650+I) 

1000 CONTINUE 
GO TO 104 


CHANGE TRAJECTORY 
210 WRITE( I0UT,220) 

220 F0RMAT(/,1X,»D0 YOU WISH TO CHANGE ANY OF THE ♦, 

$ ‘TRAJECTORY INPUT ?») 

READO I N, 50, ERR-210, END-9999) ANS 
(F(ANS.EQ.IHN)GO TO 280 

260 CONTINUE 
WRITE( I OUT, 261 ) 

261 FORMAT (/, IX, 'DO YOU WISH TO CHANGE AN EXISTING LINE ?») 
READd I N, 50, ERR-260, END-9999) ANS 
IF(ANS.EQ.1HY)G0 TO 400 

262 WRITE (I OUT, 263) 

263 FORMATOX, *D0 YOU WISH TO ADD A NEW LINE TO THE TRAJ. ?») 
READ( I IN,50,ERR=262,END=9999)ANS 

IF(ANS.EQ.1HY)G0 TO 500 

264 WRITE (I OUT, 265) 

265 FORMATdX, *D0 YOU WISH TO REMOVE A LINE FROM THE TRAJ. ?») 
READ ( 1 1 N. 50, ERR-264, END-9999) ANS 

IF(ANS.EQ.1HY)G0 TO 700 

266 WRITE (I OUT, 267) 

267 FORMATdX,* DO YOU WANT TO CHANGE THE WHOLE TRAJ. ?') 

READ( I IN, 50, ERR-266, END-9999) ANS 

I F (ANS. EQ. 1 HY )G0 TO 130 

— 

C OOPY TRAJECTORY TO FILE 
280 CONTINUE 
290 WRITE(IOUT,300) 

300 F0RMAT(1X,*D0 YOU WISH TO WRITE THIS INPUT TO A FILE ?*) 
READ ( 1 1 N, 50, ERR-290, END-9999) ANS 
IF(ANS.EQ.1HN)G0 TO 330 
310 WRITE(IOUT,320) 

320 FORMATdX,* WHAT IS THE NEW FILE NA)€ ?') 

READ ( 1 1 N, 80, ERR-3 1 0, END- 9999) FNAM1 
OPEN ( UN I T-7 , F I L E-FNAM1 , STATUS- » NEW ' ) 

325 WRITE(I0UT,326) 

326 FORMATdX,* WHAT IS THE TITLE OF THE TRAJECTORY ?•) 

READd IN,85,ERR*325,END=9999)TITLE 
WRITE(7,85)TITLE 

WRITE(7,90)BFLAG, W(50) 




1 
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IF(BFLAL.EQ.0)WRITE(7,99) (W(50+l ),W(100+I ),W( 150+1 ), 

S W(210+I),I«1,NPTS) 

IF(BFLAG.EQ.1 )WRITE(7,100) (W(50+l ),W(100+I ),W(150+I ), 
S W(210+l ),W(650+I ), |«1 ,NPT$) 

CLOSE(UNIT-7, STATUS* 'KEEP* ) 

330 WRITE(I0UT,340) 

340 F0RMAT(/,1X, 'TRAJECTORY INPUT IS COMPLETE') 

RETURN 


C 

C 


400 

410 

420 

430 

440 

441 

442 

443 


CHANGE A LINE OF TRAJECTORY 
CONTINUE 
WRITE(I0UT,410) 

FORMAT (IX, 'WHICH LINE DO YOU WISH TO CHANGE ?•) 

READC I IN, »,ERR-400, END-9999) LINE 
IF(BFLAG.EO.O)WRITE(IOUT,120)LINE,W(50+LINE),W(100+LINE), 
$ W(150+LINE),W(210+LINE) 

IF(BFLAG.EQ.1 ) WRITE ( (OUT, 120) LINE, W( 50+L INE), WOOO+L INE), 
$ W( 150+LINE),W(210+LINE),W(650+LINE) 
IF(BFLAG.EQ.0.AND.MFLAG.EQ.0)WRITE(I0UT,440)LINE 
I F (BFLAG.EQ. 0 . AND.MFLAG. EQ. 1 ) WR I TE ( I OUT, 441 ) L I NE 
I F (BFL AG. EQ. 1 . AND. MFLAG. EQ. 0 ) WR I TE ( I OUT, 442) L I NE 
I FCBFLAG.E0.1 .AND.MFLAG. EQ.1 )WRITE( I0UT,443)LINE 
F0RMAT(/,1X, 'TIME (SEC), ALTITUDE (FT), VELOCITY(FT/SEC),*, 


$ 'ANGLE OF ATTACK(DEG) »,/,1X, 12) 

FORMAT (/, 1 X, ♦ T I ME ( SEC) , ALT I TUDE ( M ) , VELOC I TY (M/SEC) , » , 
$ 'ANGLE OF ATTACK(DEG) ',/,1X, 12) 

FORMAT (/ , 1 X, • T I ME ( SEC) , ALTITUDE ( FT ) , VELOCITY (FT/SEC) , « 
$ 'ANGLE OF ATTACK (DEG), BETA (DEG) *,/,1X, 12) 

FORMAT (/, 1 X, ' T I ME ( SEC) , ALT I TUDE (M ) , VELOC I TY (M/SEC) , • , 
$ 'ANGLE OF ATTACK (DEG), BETA (DEG) »,/,1X, 12) 

I F (BFL AG. EQ. 0 ) READ ( 1 1 N, *, ERR-43 0, END«9999) W( 50+L I NE ) , 
$ W(100+LINE),W(150+L!NE),W(210+LINE) 

I F (BFLAG.EQ. 1 )READ( 1 1 N, *, ERR-430, END* 9999) W( 50+L I NE) , 
$ W(100+LINE),W(150+LINE),W(210+LINE),W(650+LINE) 

GO TO 210 


C ADD A LINE TO TRAJECTORY 
500 CONTINUE 

WRITE(I0UT,510) 

510 F0RMAT(1X,'WHAT IS THE NUNBER OF THE LINE YOU WISH TO ADD ?') 
READ( I IN, *, ERR-500, END-9999) L I NE 
DO 600 K*NPTS,LINE,-1 
J*K+1 


W(50+J)-W(50+K) 

W(100+J)*W(100+K) 

W( 150+J)*W( 150+K) 

W(210+J)*W(210+K) 

W( 650+J )*W( 650+K ) 

600 CONTINUE 

520 IF(BFLAG.EQ.O. AND.MFLAG. EQ.O) WRI TE ( I0UT,440)L I NE 
IF(BFLAG.EQ.0.AND.MFLAG.EQ.1 )WRITE( I0UT,441 >L INE 
IF(BFLAG.EQ.1 .AND.MFLAG. EQ.O)WRITE( I0UT,442)LINE 
I F(BFLAG. EQ. 1 .AND.MFLAG. EQ. 1 ) WR I TE ( I0UT,443 ) L I NE 
I F (BFL AG. EQ. 0 ) READ ( 1 1 N, *, ERR-520, END-9999) W( 50+L INE), 
S W( 100+L INE), W( 1 50+L INE), W(210+L INE) 
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IF(BFLAG.EQ.1 )READ ( I IN#*, ERR-520# END-9999) W( 50+L I NE ) , 
$ W ( 1 00+L I NE ) # W ( 1 50+L I NE J , W ( 21 0+L I NE ) # W( 650+L I NE ) 
NPTS-NPTS+1 
W ( 50) -FLOAT ( NPTS ) 

60 TO 210 


DELETE A LINE OF TRAJECTORY 
700 CONTINUE 

WRITE( IOUT#710) 

710 FORMAT (IX,' WHAT IS THE NUI€ER OF THE LINE YOU WISH TO DELETE ?') 
READ ( 1 1 N# «# ERR-7 00# END- 9999 ) L I NE 
DO 800 J-LINE# NPTS-1 
K-J+1 

W(50+J)«W(50+K) 

W(100+J)=W(100+K) 

W ( 1 50+ J ) ®W ( 1 50+K ) 

W(210+J)*W(210+K) 

W ( 650+ J ) ®W ( 650+K ) 

800 CONTINUE 
NPTS® NPTS-1 
W ( 507*FL0AT ( NPTS ) 

60 TO 210 
8888 CONTINUE 

WRITE( IOUT,350)FNAM1 
350 FORMAT ( IX, • CANNOT OPEN *,A20) 

60 TO 60 
9999 CONTINUE 
CALL EXIT 
END 


OF POOR QUALITY 
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SUBROUTINE TRANS 




RUTINE TO CHOOSE TRANSITION OPTION 
COMMON/ HARRAY/HC 700) 

COMMON/ UN IT/ 1 IN, I OUT 

COMMON/MI SC/NC, JFKS,TFLA6,MFLAG, ANSI 


CWGifrJAt. PAGE |9 
OF POOR QUALITY 


INTEGER ANS,TFLAG,ANS1 

C MAY USE TRANSITION DATA FROM PREVIOUS CASE IF JFKS ■ 1 OR 3 
IF( JFKS.NE.1 .AND. JFKS.NE.3)G0 TO 9 
NCM1-NC-1 

1 HRITE( I0UT,2)NC, NCM1 

2 F0RMATC1H1 ,’ IS CASE • # I3 # * TRANSITION DATA THE SAME », 

$ ’FOR CASE »,l 3,» ?*> 

READCI IN, 130, ERR-1, END-9999) ANS 
I F (ANS.EQ. 1 HY )RETURN 

9 CONTINUE 

10 NR I TEC I OUT, 20) 

20 F0RMAT(///,10X,» TRANSITION OPTIONS*,//, IX, 

$ ’OPTIONS 1, TIME DEPENDENT} LAM TO TURB',/, 

$ 1 IX, *2. TIME DEPENDENT} TURB TO LAM’,/, 

$ 11X,»3. REYNOLDS NO. DEPENDENT’,/, 

$ 1 IX, *4. RE-THETA’,/, 

S 11X,»5. MDAC-E TRANSITION’,/, 

S 11X,»6. MDAC-E TABLE LOOK-UP*,/, 

$ 11X, *7. NAR RE VS ME TABLE LOOK-UP’,/, 

S 11X,»8. RE-THETA/ME’,//, 

$ IX, ’OPTION SELECTED ?•) 

READCI IN, ”, ERR- 1 0, END-9999) H ( 27 ) 

IT- IFIXCHC27) ) 

GO TO (100,200,300,400,500,600,700,800), IT 
GO TO 900 

TRANSITION OPTION 1 

100 HRITEC ICUT, 105) 

105 F0RMAT(//,1X,'1. TIME DEPENDENT} LAM TO TURB*, 

S //, IX, ’TRANSITION BEGINS AT TI)C(SEC) ?») 

READC 1 1 N, ”, ERR-1 00, END-9999) W( 14) 

110 WR I TE ( I OUT, 115) 

115 FORMATC IX, ’FULLY TURBULENT AT TIME(SEC) ?•) 

READCI IN, *,ERR-110,END-9999)H(20) 

IFCWC29) .EQ.O.O)GO TO 118 

116 WRITECIOUT, 1 17) 

117 FORMATC IX, ’A VIRTUAL ORIGIN CORRECTION CANNOT BE MADE *, 
S’NITH THIS OPT ION.',/, IX, ’RESET VIRTUAL ORIGIN OPTION ?») 

READCI IN, 130,ERR»1 16, END-9999) ANS 
I F CANS. EQ. 1 HY )H( 29) -0 .0 

118 CONTINUE 

120 HR I TEC I OUT, 125) 

125 FORMAT (IX, ’ANY CHANGES ?•) 

READC I IP, 130, ERR-1 20, END-9999) ANS 

130 FORMAT (A1) 

IF(ANS.EQ.1HY)G0 TO 10 
GO TO 1000 

C TRANSITION OPTION 2 

200 HR I TEC I OUT, 205) 


V 
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205 


210 

215 


216 

217 


218 

220 


F0RMAT<//,1X,'2. 

$ 'FULLY TURBULENT 


TIME DEPENDENTi TURB TO LAM',//, IX, 
rtT TIME (SEC) ?•) 


READ( I IN,»,ERR«200,END«9999)W< 14) 
WRITE(I0UT,215) 

FORMAT ( 1 X, ' L AM I N AR FLOW BEGINS AT TIME(SEC) ?•) 
REAO(| IN,*,ERR-210,END«9999)W(20) 

IF(W(29) «EQ.O.O)GO TO 218 
WRITE<I0UT,217) 


FORMAT < IX, » A VIRTUAL ORIGIN CORRECTION CANNOT BE MADE ' 
$*WITH THIS OPTION. »,/,1X, 'RESET VIRTUAL ORIGIN OPTION ? 
READd IN,130,ERR-216,END-9999)ANS 
IF(ANS.EQ.1HV)W(29)-0.0 
CONTINUE 
WRITE (I OUT, 125) 

REAO( 1 1 N, 1 30, ERR-220, END-9999) ANS 
IF(ANS.EQ.1HY)G0 TO 10 
GO TO 1000 


' TRANSITION OPTION 3 

300 WRITE(I0UT,305) 

305 F0RMAT(//,1X,»3. REYNOLDS NO. DEPENDENT*, //,1 OX, 

! * RE- LOCAL REYNOLDS NO. FOR FLAT PLATE OPTIONS' ,/,1 OX, 
$ 'RE“(RE-INF)D FOR SWEPT CYLINDER OPTIONS* ,/, IX, 

$ 'HIGHEST LAMINAR RE ?») 

READ( 1 1 N, «, ERR-300, END-9999) W( 14) 

310 WRITE(I0UT,315) 

315 FORMAT ( IX, ♦ LOWEST TURBULENT RE ?•) 

READ (I I N, *, ERR-3 1 0, END- 9999) W ( 20 ) 

320 WRITE (I OUT, 125) 

READd IN, 130,ERR»320,END«9999)ANS 
(F(ANS.EQ.IHY)GO TO 10 
GO TO 1000 


’ TRANSITION OPTION 4 

400 WRITE(I0UT,405) 

405 F0RMAT(//,1X, *4. RE-THETA DEPENDENT*,//, IX, 

$ 'HIGHEST LAMINAR VALUE ?») 

READ< 1 1 N, «, ERR-400, END-9999) W( 14) 

410 WRITE(I0UT,415) 

415 FORMAT( IX, 'LOWEST TURBULENT VALUE ?') 

READ (I I N, *, ERR-4 1 0, END-9999) W ( 20 ) 

420 WRITE(I0UT,125) 

READd IN, 130, ERR-420, END-9999) ANS 
IF(ANS.E(\ HY)GO TO 10 
GO TO 1000 

; TRANSITION OPTION 5 

500 WRITE (I OUT, 505) 

505 F0RMAT(//,1X,'5. MDAC-E TRANSITION: RE-THETA/', 

$ '<ME«(RH0«V/MU)**.2)',//,1X, 

$ 'HIGHEST LAMINAR VALUE ?') 

READC 1 1 N, ", ERR-500, END-9999) W( 14) 

510 WRITE(I0UT,515) 

515 FORMAT( IX, 'LOWEST TURBULENT VALUE ?') 

READ( I IN,*,ERR-51 0,END=9999)W(20) 

520 WRITEdOUT, 125) 

READ( 1 1 N, 1 30, ERR-520, END-9999) ANS 
IF(ANS.EQ.1HY)G0 TO 10 


C-40 


ORIGINAL PAGE 19 
OF POOR QUALITY 


GO TO 1000 


TRANSITION OPTION 6 


C 

600 WRITEdOU*. . 05) 

605 FORMAT <//,1X,»6. MDAC-E TABLE LOOK-UP' ,//, IX, 

S 'OPTIONAL INPUT* (L)FULLY TURB/(L)TRAN ONSET', 

$ /,1X, INPUT THE RATIO OR 0.0 **»•'",/, IX, 
$ 'INPUT DESIRED ?•) 

READd IN,*,ERR-60?,END-9999)W(28) 

610 WRITE( I OUT, 125) 

READ (I I N, 1 30, ERR-6 1 0, END-9999) ANS 
IF(ANS.EQ.1HY)G0 TO 10 
GO TO 1000 

C TRANSITION OPTION 7 

700 WRITE (I OUT, 7 05) 

705 F0RMAT(//,1X, *7. NAR RE VS. ME TABLE LOOK-UP', 
$ //,1X, 'NO INPUT REQUIRED') 

710 WRITE(I0UT,125) 

READd I N, 1 30, ERR-7 1 0, END-9999) ANS 
IF(ANS.EQ.1HY)G0 TO 10 
GO TO 1000 

C TRANSITION OPTION 8 

800 WRITE (I OUT, 805) 

805 FORMAT (//, IX, *8. RE-THETA/ME DEPENDENT',//, IX, 

$ 'HIGHEST LAMINAR VALUE ?*> 

READd I N,«,ERR»800, END-9999) W( 14) 

810 WRITE( I OUT, 81 5) 

815 FORMAT < IX, 'LOWEST TURBULENT VALUE ?•) 

READ( I IN, ", ERR-81 0, END-9999) W( 20) 

820 WRITE( IOUT, 1 25) 

RE AD (I I N, 1 30, ERR-820, END-9999) ANS 
IF(ANS.EQ.1HY)G0 TO 10 
GO TO 1000 

C BAD OPTION 

900 CONTINUE 

WRITE (I OUT, 910) 

910 FORMAT (//, IX, 'BAD OPTION NU)©ER«,//) 

GO TO 10 
1000 CONTINUE 
RETURN 

9999 CONTINUE 
CALL EXIT 


SUBROUTINE WNDTUN 
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C ROUTINE FOR INPUTING ATMOSPHERIC DATA FOR WIND TUNNEL OPTION 
COMMON/ WARRAY/W( 700) 

COMMON/UN IT/ 1 IN, I OUT 

COMMON/MI SC/NC, JFKS, TFL AG, MFL AG, ANSI 

c 

INTEGER ANS, TFL AG, ANSI 

C NPTS - NUMBER OF TIME DEPENDENT TABLE ENTRIES 
C MAY USE DATA FROM PREVIOUS CASE IF JFKS - 1 OR 3 
I F( JFKS. NE.1. AND. JFKS. NE. 3)60 TO 9 
NCM1*NC-1 

1 WR I TE ( I OUT, 2 ) NC, NCM1 

2 F0RMAT(1H1,*IS CASE ',13,' WIND TUNNEL DATA SAI€ AS • 

S 'FOR CASE *,I3,» ?•) 

READ< 1 1 N, 80, ERR* 1 ,END-9999)AN$ 

I F(ANS.EQ.IHY) RETURN 
9 CONTINUE 

NPTS*W(50)+.0001 
10 WRITE (I OUT, 20) NPTS 

20 FORMAT(///,» WIND TUNNEL OPTION',//, 

$ * INPUT STATIC TEMPERATURE AND PRESSURE AS A FUNCTION ', 
$ 'OF TIME.',/,' Tire AND FREESTREAM VEL. ARE INPUT ♦, 

$ 'IN TRAJ. DATA WITH ALT. SET * 0.0', /,1X 
S 6( ' -'), 12,' VALUES REQUIRED', 6( ' -•),/) 
IF(MFLAG.EQ.O)WRITE( I0UT,30) 

I F(MFLAG.EQ.1 ) WR I TE < i OUT, 31) 

30 FORMAT d X, » T- 1 NF ( R ) , P- 1 NF < LB/ SFT ) • ) 

3 1 FORMAT ( IX, 'T- I NF (K ) , P- I NF ( NEWTON/ SO. M) ' ) 

DO 1000 1*1, NPTS 

40 WRITE( I0UT,50) I 
50 FORMATdX, 12) 

READ( I IN, *,ERR=40,END=9999)W(450+I ),W(500+I ) 

1000 CONTINUE 
60 WRITE(I0UT,70) 

70 FORMATdX, 'ARE ALL INPUTS CORRECT ?») 

READ( I IN,80,ERR*60,END a 9999)ANS 
80 FORMAT (A1 ) 

IF(ANS.EQ.1HN)G0 TO 10 
RETURN 

9999 CONTINUE 
CALL EXIT 


SUBROUTINE UNITS 


OF PQt'tft 

C ROUTINE FOR CONVERTING TO METRIC UNITS 

C COPIES WORKING ARRAY W INTO OUTPUT ARRAY WW 
COMMON/ WARRAY/W ( 700 ) 

COMMON/WWARAY/WW( 700) 

COMMON/ WWWARY/WWW( 700) 

COMMON/UN IT/ 1 IN, I OUT 

COMMON/MI SC/NC,JFKS,TFLAG,MFLAG, ANSI 

c 

INTEGER TFLAG,ANS1 

C 

00 1000 1*1 ,700 
WWW( I )*WW( I ) 

WW( I )«W( I ) 

1000 CONTINUE 

C MFLAG - 0 ENGLISH UNITS 
IF(MFLAG.EQ.O)GO TO 2000 
WW(12)*W( 12) *3.28084 
WW(13)*W(13)*3. 28084 
WW( 22) *W( 22) *3.28084 
WW( 25) *W( 25) *3.28084 
WW(202)*W(202) *3.28084 
WW(205)-W< 205) *3.28084 
WW( 354) -W( 354) *3.28084 
DO 1 1*571,580 

1 WW( I )*W( I )*3. 28084 
DO 2 1*101,200 

2 WW( I )*W( I ) *3 .28084 
DO 3 1*401,450 

3 WW( I )-W( I )*3. 28084 
DO 4 1-451,500 

4 WW( I )*W( I )*1 .8 
DO 5 1*501,550 

5 WW< I )-W( I ) *0.0208855 
WW(316)*W(316)*8.81 143E-5 
WW(21)*W( 21) *4.301 89E-4 
WW(24)*(W(24)»1 .8-459.7) 

2000 CONTINUE 
RETURN 
END 


oi 'au rv 
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,4 


C 

C 

c 

c 


c 

c 


SUBROUTINE OUTPUT 

ROUTINE TO CREATE INPUT FILE FOR LANMIN 

FNAM4 - OUTPUT FILENAME (INPUT FILE FOR LANMIN) 
COMMON/WARRAY/W (700) 

COMMON/ W WAR AY/ WW( 700) 

COMMON/ WWWARY/ WWW ( 700 ) 

COMMON/UN IT/ 1 IN, I OUT 

COMMON/ POCEF F/TMCP ( 50 ) , TCPM ( 50 ) , NCPMT 

OOMMON/T I TLE/T I TL 1 

COMMON/MI SC/NC, JFK$,TFLAG,MFLAG, ANSI 
COMMON/ STRM/ I STRM, NNC 

INTEGER ANS,TFLAG,ANS1 
CHARACTER*20 FNAM4 
CHARACTER«72 TITL1 
DIMENSION IS(5) 


IF(NC.GT.1 .0)G0 TO 26 
OBTAIN OUTPUT FILE NAME AND OPEN 
10 WRITE( I0UT,20) 

20 F0RMAT(1H1,1 OX, ♦ CREATE OUTPUT FILE',//,1X, 

S 'WHAT IS THE NAME OF THE OUTPUT FILE TO BE CREATED ?') 
READ ( 1 1 N, 25, ERR= 1 0, END=9999) FNAM4 
OPEN ( UN I T*8, F I L E*FN AM4 , STATUS® » NEW • , ERR*8^88) 

25 F0RMAT(A20) 

c Zm : T v T0 ^™;s F i« t '' oooR ° 1 ng to pro6ram ,nput °° ntro1 

26 CONTINUE 

IF(JFKS.NE.O)GO TO ( 100,9999,30,30,1 00,30), JFKS 
IF(TFLA6.NE.t )G0 TO 33 ' 

WRITE(I0UT,31)NC 

F0RMAT(1X,»D0 YOU WISH TO CHANGE THE TITLE FOR CASE M3) 
READ(IIN.32,ERR«30,END=9999)ANS ' 

FORMAT (At ) 

I F(ANS.EQ.1HN)G0 TO 52 
WRITE( IOUT,40)NC 

F0RMAT(/,1X,»WHAT IS THE TITLE FOR CASE »,I3,» ? »,5X. 

$ '(NOTE: 72 CHAR. LIMIT)') 

READO IN,50,ERR*33,ENP*9999)TITL1 
FORMAT (A72) 

WRITE(8,55)TITL1 ~ f,TLE 

FORMAT (A72) 

C TIMING PARAMETERS AND PRINT CONTROL 

WRITE(8,60) (WW( I ), 1*1 ,8) 

FORMAT (3F 20. 3) 

IF( JFKS.EQ.3.0R. JFKS.EQ.6)G0 TO 100 

W TINt DEPENDENT traj table entries - 
WRITE (8,70) WW( 50) 

F0RMAT(3F20.3,2F 1 0.3 ) 

N«WW(50)+.001 
NT®50+N 

DATA 
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31 
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33 
40 
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52 
55 


60 
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Oi '--jk; ; ? 

OF f-Coi? ‘ 


WRITE (8,70) (WW(K),WW(K+50),WW(K+100),WW(K+160) ,WW(K+600), 
$K«51,NT) 

CASE DATA 

IF IZFLAG • 1 THEN THE ZERO VALUE FOR THE W MUST BE PRINTED 
100 >9 
110 JS*1 

120 IF(JS.GT.5)G0 TO 130 
IZFLAG-0 

IF( J.6T.49.AND. J.LT.201 ) J*201 
IF(J.GT.210.AND. J.LT.261 )J*261 
IF(J.GT.700)GO TO 200 
IF(J.EQ.641)J-J+1 

IF(WW(J).EQ.0.0.AND.WWW(J).NE.0.0)IZFLAG«1 

IF(WW(J).NE.0.0.0R.IZFLAG.EQ.1)IS(JS)«J 

IF(J.EQ.649.AND.WW( J).GT.O.O)GO TO 300 

IF(WW(J).GT.1000.0.OR.WW(J).LT.1.0.AND.WW(J).6T.0.0)GO TO 400 

IF(WW(J).NE.O.O.OR. IZFLAG.EQ.1 )JS»JS+1 

J“J+1 

GO TO 120 

C WRITE CASE DATA TO OUTPUT FILE 

130 WRITE(8,80)(IS(K),WW(IS(K)),K«1,5) 

80 FORMAT<5( I3,F9.4,1X)) 

81 FORMAT(5( I3,F9.7,1X)) 

82 F0RMAT(5(I3,F9.1,1X)) 

GO TO 110 

C WRITE PROGRAM INPUT CONTROL PARAMETERS FOR CASE 

200 IS(JS)«641 

WRITE(8,80) ( IS(K),WW( IS(K) ),K«1 , JS) 

C (MECK FOR STREAMLINE CASES 

I F ( I STRM. EQ. 1 )G0 TO 500 
250 CONTINUE 

I F ( I NSERT. EQ.O) JFKS s WW(641 )+.0001 
I F ( JFKS . EQ. 2 ) CLOSE ( UN I T*8 ) 

RETURN 

C PRESSURE COEFFICIENT INPUT TABLES 

300 CONTINUE 

WRITE (8,80) ( IS(K),WW( IS(K) ),K S 1, JS) 

WRITE(8,90)NCPMT 
90 F0RMATO3) 

WR ITE (8,95) (TMCP (K ) ,TCPM(K ) ,K»1 , NCPWT) 

95 F0RMAT(2F10.6) 

J-J+1 
GO TO 110 

C - OPTIONAL FORMAT FOR WRITING CASE DATA TO OUTPUT FILE 

400 CONTINUE 

IF(JS-1 .GE.1 ) WRITE (8, 80) ( IS(K),WW( IS(K) ),K*1 , JS-1 ) 

IF(WW( J) »LT.1000»)WRITE(6,81 ) J,WW( J) 

I F (WW( J). GE.1 000.) WRITE (8,82) J,WW(J) 

J*J+1 
GO TO 110 

C — — - COPY STREAM INE CASES FROM TEMPORARY STREAM INE FILE — 

500 CONTINUE 

READ(4,80)J1,W(J1 ), J2,W( J2), J3,W( J3), J4,W( J4),J5,W( J5) 
READ(4,80) J6,W( J6) 

WRITE(8,80)J1,W(J1),J2,W(J2),J3,W(J3),J4,W(J4),J5,W(J5) 
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IF(W( J1 ).EQ.NNC-1 )G0 TO 510 
WRITE(8,80)J6,W(J6) 

GO TO 500 

C CHOOSE PROGRAM INPUT CONTROL PARAMETER FOR LAST STREAM I NE CASE 
510 CONTINUE 

WRITE(I0UT,520) 

520 F0RMAT(//, 1 0X, » INPUT CONTROL FLAG*,//, 

$ IX, M. NEW CASE FOLLOWS USING TITLE,TIMING AND TRAJ. ', 

S /,6X,'DATA FROM PREVIOUS CASE',/, 

$ IX, *2. END OF INPUT (LAST CASE)',/, 

$ IX, '3. NEW CASE FOLLOWS USING TRAJ. DATA FROM PREVIOUS ', 

S 'CASE. »,/,6X,' NEW TITLE AND TIMING. INITIAL CASE DATA », 

$ 'UNCHANGED.',/, 

$ 1X,»4. NEW CASE FOLLOWS USING NEW TITLE, TIMING, TRAJ. ', 

$ 'AND CASE DATA. *,/,6X,'( INITIALLY ZERO W ARRAY)',/, 

S IX, »5. SAME AS (1) EXCEPT ZERO ALL CASE DATA FROM ', 

$ 'PREVIOUS CASE*,/, 

$ IX, *6. SAME AS (3) EXCEPT INITIALIZE ZERO ALL TIMING ', 

S 'AND CASE DATA',//, 

$ IX, 'OPTION SELECTED ?•) 

READ (I IN,*, ERR*5 1 0, END*9999) W ( J6 ) 

WRITE C8,80)J6,W(J6) 

CLOSE (UN I T=4, STATUS* 'DELETE' ) 

NC*W( 647) +.0001 lffl 

GO TO 250 ORIGINAL PAGE Fa 

8888 CONTINUE OF POOR QUALITY 

WR ITE ( I OUT, 8889 ) FNAM4 

8889 FORMAT (/, IX, 'CANNOT OPEN »,A20,/) 

GO TO 10 

9999 CONTINUE 
CALL EXIT 
END 
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